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Chapter 

Introduction 



A new field of discrete differential geometry is presently emerging on the 
border between differential and discrete geometry. Whereas classical differ- 
ential geometry investigates smooth geometric shapes (such as surfaces) , and 
discrete geometry studies geometric shapes with finite number of elements 
(such as polyhedra), the discrete differential geometry aims at the devel- 
opment of discrete equivalents of notions and methods of smooth surface 
theory. Current interest in this field derives not only from its importance 
in pure mathematics but also from its relevance for other fields like com- 
puter graphics. Recent progress in discrete differential geometry has lead, 
somewhat unexpectedly, to a better understanding of some fundamental 
structures lying in the basis of the classical differential geometry and of the 
theory of integrable systems (this is schematically presented on Fig. ^) . The 
goal of this book is to give a systematic presentation of current achievements 
in this field. 

The classical period of development of surface theory resulted in the 
beginning of the 20-th century in an enormous amount of knowledge on 
numerous special classes of surfaces, coordinate systems and their transfor- 
mations, which is summarized in extensive volumes by Darboux (DaltaDa2j . 
Bianchi [BJ etc. One can say that the local differential geometry of spe- 
cial classes of surfaces and coordinate systems has been completed during 
this heroic period. Mathematicians of that era have found most (if not all) 
geometries of interest and knew nearly everything about their properties. 
It was observed that such special geometries as minimal surfaces, surfaces 
with constant curvature, isothermic surfaces, orthogonal and conjugate coor- 
dinate systems, Ribaucour sphere congruences, Weingarten line congruences 
etc. have many similar features. Among others we mention here Backlund 
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Figure 1: Fundamental consistency principle of the discrete differential ge- 
ometry as a conceptual basis of the differential geometry of special surfaces 
and of the integrability. 
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and Darboux type transformations, with remarkable permutability proper- 
ties investigated mainly by Bianchi, and existence of special deformations 
within the class (associated family). Geometers realized that there should 
exist a unifying fundamental structure behind all these common proper- 
ties of quite different geometries. And they were definitely looking for this 
structure |jHE2j. 

Much later, after advent of the solitons theory in the the last quarter of 
the 20-th century, these common similar features were recognized to be asso- 
ciated to integrability of the underlying differential equations. However, the 
current status of the notion of integrability remains unsatisfactory from the 
point of view of a mathematician. There exists no commonly accepted defi- 
nition of the integrability (as the title of the volume "What is integrability?" 
jZj clearly demonstrates). Different scientists suggest different properties as 
defining ones. Usually, one just refers to some common features as Darboux 
transformations etc., exactly like the differential geometers of the classical 
period did. 

A progress in understanding of the unifying fundamental structure the 
classical differential geometers were looking for, and simultaneously in un- 
derstanding of the very nature of integrability, came from the efforts to 
discretize all these theories. It turns out that many sophisticated prop- 
erties of differential- geometric objects find their simple explanation within 
the discrete differential geometry. The early period of its development is 
documented in the works of Sauer and Wunderlich 'Sa' 'W|. The modern 
period began with the work by Bobenko and Pinkall BobPl|lBobP2j and by 
Doliwa and Santini |DoSll ICDSj . A closely related development of the spec- 
tral theory of difference operators on graphs was initiated by Novikov with 
collaborators |NoDl INoll INo2j . see also 'DyN' for a further development 
of a discrete complex analysis on triangulated manifolds. Discrete differ- 
ential geometry deals with multidimensional discrete nets (i.e., maps from 
the regular cubic lattice into or some other suitable space) speci- 
fied by certain geometric properties. In this setting, discrete surfaces ap- 
pear as two-dimensional layers of multidimensional discrete nets, and their 
transformations correspond to shifts in the transversal lattice directions. 
A characteristic feature of the theory is that all lattice directions are on 
equal footing with respect to the defining geometric properties. Due to this 
symmetry, discrete surfaces and their transformations become indistinguish- 
able. We associate such a situation with the multidimensional consistency 
(of geometric properties, resp. of equations which serve for their analytic 
description). The multidimensional consistency, and therefore the existence 
and construction of multidimensional discrete nets, relies just on certain 



8 



CHAPTER 0. INTRODUCTION 



incidence theorems of elementary geometry. 

Conceptually, one can think of passing to a continuum limit by refin- 
ing mesh size in some of the lattice directions. In these directions the net 
converges to smooth surfaces whereas those directions that remain discrete 
correspond to transformations of the surfaces (see Fig. [JJ . Differential ge- 
ometric properties of special classes of surfaces and their transformations 
follow in this way from (and find their simple explanation in) the elemen- 
tary geometric properties of the original multidimensional discrete nets. In 
particular, difficult classical theorems about permutability of the Backlund- 
Darboux type transformations (Bianchi permutability) for various geome- 
tries follow directly from the symmetry of the underlying discrete nets, and 
are therefore built in the very core of the theory. Thus the pass from differ- 
ential geometry to elementary geometry via discretization (or, in an opposite 
direction, the derivation of the differential geometry from the discrete dif- 
ferential geometry) leads to enormous conceptual simplifications, and the 
true roots of the classical theory of special classes of surfaces are found in 
various incidence theorems of elementary geometry. However, these elemen- 
tary roots become deeply hidden in the classical differential geometry, since 
the continuum limit from the discrete master theory to the classical one 
is inevitably accompanied by a break of the symmetry among the lattice 
directions, which always yields essential structural complications. 



Figure 2: From the discrete master theory to the classical theory: surfaces 
and their transformations appear by refining two of three net directions. 

To give a rigorous justification of this philosophy, one needs to prove con- 
vergence of the procedure just described. Theorems of this kind were lacking 
in the literature until recently. The first results of this kind have been proven 
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by the authors in a common work with Matthes |BobMaS 11 lBobMaS2j . This 
makes the general philosophy of the discrete differential geometry to a firmly 
established mathematical truth for several important classes of surfaces and 
coordinate systems, like conjugate nets, orthogonal nets, including general 
surfaces parametrized along curvature lines, surfaces with constant negative 
Gaussian curvature, and general surfaces parametrized along asymptotic 
lines. For some other classes, like isothermic surfaces, the convergence re- 
sults still wait to be rigorously established. 

But finding simple discrete explanations for complicated differential ge- 
ometric theories is not the only outcome of this development. It is well 
known that differential equations which analytically describe interesting spe- 
cial classes of surfaces are integrable (in the sense of the theory of integrable 
systems), and, conversely, many of interesting integrable systems admit a 
differential-geometric interpretation. Having identified the roots of the inte- 
grable differential geometry in the multidimensional consistency of discrete 
nets, we are led to a new (geometric) understanding of the integrability it- 
self. First, we adhere to the viewpoint that the central role in this theory 
belongs to discrete integrable systems. In particular, all the great variety 
of integrable differential equations can be derived from several fundamen- 
tal discrete systems by performing different continuous limits. Further, and 
more important, we come to a constructive and almost algorithmic defini- 
tion of integrability of discrete equations as their multidimensional consis- 
tency, introduced by the authors in |BobSulj (and independently in |Nij| ). 
It turns out that this definition captures enough structure to yield such 
traditional attributes of integrable equations as zero curvature representa- 
tions and Backlund-Darboux transformations (which, in turn, serve as the 
basis for applying analytic methods like inverse scattering, finite gap inte- 
gration, Riemann-Hilbert problems, etc.). A continuous counterpart (and 
consequence) of the multidimensional consistency is the well-known fact that 
integrable systems never appear alone but are organized into hierarchies of 
commuting flows. 

The geometric way of thinking about the discrete integrability has also 
led to introducing novel concepts into the latter. One of the reasons to con- 
sider discrete integrable systems on the regular square lattice is the desire 
to have a proper model for parametrized surfaces. However, an immanent 
and important feature of various parametrizations of surfaces is the exis- 
tence of distinguished points, where the combinatorics of coordinate lines 
change (like umbilic points, where the combinatorics of the curvature lines 
is special). This compels us to introduce quad-graphs, which are cell de- 
compositions of topological two-manifolds with quadrilateral faces. Their 
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elementary building blocks are still quadrilaterals but are attached to one 
another in a manner which can be more complicated than in 7?. This no- 
tion has been introduced into the context of discrete differential geometry in 
|BobP3j , and a systematic development of the theory of integrable systems 
on quad-graphs has been undertaken by the authors in [BobSulj . 

The structure of the book follows the logic of this introduction. We 
start in Chapter ^ with an overview of some classical results of the surface 
theory, focusing on transformations of surfaces. The geometries considered 
here include general conjugate and orthogonal nets in spaces of arbitrary 
dimension, asymptotic nets on general surfaces, as well as special classes of 
surfaces, like isothermic ones and surfaces with constant negative Gaussian 
curvature. There are no proofs in this chapter: on one hand, these tedious 
analytic proofs can be found in the original literature, and on the other 
hand, the discrete approach which we develop in the subsequent chapters 
will lead to conceptually transparent and technically much simpler proofs. 

In Chapter 121 we define and investigate discrete analogs of the classical 
geometries discussed in the previous chapter, focusing on the idea of multi- 
dimensional consistency of discrete nets. It turns out that all these discrete 
analogs are reductions of discrete conjugate nets, which are multidimensional 
nets with the combinatorics of Z"^, consisting of planar quadrilaterals. Im- 
posing additional constraints on the geometry of elementary quadrilaterals, 
one comes to discrete orthogonal nets, discrete asymptotic nets, discrete 
isothermic surfaces, discrete surfaces with constant negative Gaussian cur- 
vature etc. 

Then in Chapter 13] we develop an approximation theory for hyperbolic 
difference systems, which is applied to derive the classical theory of smooth 
surfaces as a continuum limit of the discrete theory. We prove that dis- 
crete nets of Chapter 121 approximate the corresponding smooth geometries of 
Chapter n simultaneously with their transformations. Bianchi's and Eisen- 
hart's permutability theorems for transformations appear in this approach 
as simple corollaries. 

In Chapter |1] we formulate the concept of multi-dimensional consis- 
tency as a defining principle of integrability. We derive basic features of 
integrable systems, such as the zero curvature representation, Backlund- 
Darboux transformations, from the consistency principle. At this point, the 
theory makes an interesting conceptual turn. First, we generalize the under- 
lying combinatorial structure of the two-dimensional theory from the regular 
square lattice to arbitrary quad-graphs, i.e., cell decomposition of surfaces 
with all quadrilateral faces. Introducing these generalized combinatorics has 
been partly motivated by the desire to have proper discrete models for ge- 
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ometrically significant special points of parametrized surfaces, like umbilic 
points in the curvature line parametrization. But then the multidimen- 
sional consistency allows us to regard integrable systems on quad-graphs as 
systems on regular square lattices restricted to quadrilateral surfaces. 
At this point, interesting interrelations with the discrete geometry and the 
combinatorial analysis come onto the scene. 

Finally, in Chapters |S1 these ideas are applied to discrete complex 
analysis. We study Laplace operators on graphs, discrete harmonic and 
holomorphic functions. The linear discrete complex analysis appears here 
as a linearization of the theory of circle patterns. The consistency principle 
allows us to single out distinguished cases where we obtain more detailed 
analytic results (like Green's function and isomonodromic special functions). 

Essential parts of this book are based on results obtained jointly with 
Vsevolod Adler, Tim Hoffmann, Daniel Matthes, Christian Mercat and Ul- 
rich Pinkall. We warmly thank them for inspiring collaboration. 
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Chapter 1 

Classical differential 
geometry 

In this chapter we discuss some classical results of the differential geome- 
try of nets (parametrized surfaces and coordinate systems) in M^, mainly 
concentrated around the topics of transformations of nets and of their per- 
mutability properties. This classical area was very popular in the differential 
geometry of the 19th and of the first quarter of the 20th century, and is well 
documented in the fundamental treatises |Bi[ IDall IDa2l IE1[ IE21 ITzij and 
others. Our presentation mainly follows these classical treatments, of course 
with modifications which reflect our present points of view. We do not trace 
back the exact origin of the concrete classical results: often enough this 
turns out to be a complicated task in the history of mathematics which still 
waits for its competent investigation. 

For the classes of nets described by essentially two-dimensional systems 
(special classes of surfaces such as surfaces with a constant negative Gaussian 
curvature or isothermic surfaces), the permutability theorems, mainly due 
to Bianchi, are dealing with a quadruple of surfaces (depicted as vertices of 
a so-called Bianchi quadrilateral). Given three surfaces of such a quadruple, 
the fourth one is uniquely deflned, see Theorems I1.2UI and 11.241 

For the classes of nets described by essentially three-dimensional systems 
(conjugate nets; orthogonal nets, including general surfaces parametrized by 
curvature lines; Moutard nets; general surfaces parametrized by asymptotic 
lines), the situation is somewhat different. The corresponding permutability 
theorems (Theorems 11.31 H"^ II. Ill and II. lb')) consist of two parts. The first 
parts present the traditional view and are dealing with Bianchi quadrilater- 
als. In our opinion, this is not the proper setting in the three-dimensional 
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context, and the non- uniqueness of the fourth net in these theorems re- 
flects this. The natural setting for permutability is given in the second 
parts, where the permutability is associated with an octuple of nets, de- 
picted as vertices of a combinatorial cube, so that the eighth net is uniquely 
determined by other seven ones. We found the first instance of this kind 
of statement in |E2j . §24 ("extended theorem of permutability" for conju- 
gate nets), and propose therefore to term such an octuple of nets as an 
Eisenhart cube; in the modern literature on integrable systems, this kind of 
permutability theorems in classical differential geometry was (re-) discovered 
in |GTj (for orthogonal nets). Our discrete philosophy makes the origin of 
such permutability theorems quite transparent. 

A few remarks on notations: we denote independent variables of a net 
/ : ^ by u = (ui, ...,Um)e M™, and we set di = d/dui. We write 
"^h-./is = {u ^ IK™' '■ Ui = 0, i ^ ii,...,is} for s-dimensional coordinate 
planes (coordinate axes, if s = 1). We always suppose that the dimension 
of the ambient space > 3. 

1.1 Conjugate nets and their transformations 

Conjugate nets. This classical notion can be defined as follows. 

Definition 1.1 A map / : M™ ^ is called an m-dimensional conjugate 
net in M.^ , if didjf € span{dif,djf) at any u G and for all pairs 1 < 
i ^ j < m. 

A parametrized surface in the three-space (m = 2, = 3) is a conjugate 
net, if its second fundamental form is diagonal in this parametrization. Such 
a parametrization exists for a general surface in the three-space. It is impor- 
tant to note that Definition ll.il as well as Definition 11.21 below, are dealing 
with projectively invariant notions only, and thus belong to the projective 
differential geometry. In this setting the space where conjugate net lives 
should be interpreted as an affine part of P(M^"'~-'^). 

^From Definition 11.11 there follows that conjugate nets are described by 
the following (linear) differential equations: 

didjf = Cjidif + Cijdjf, i 7^ j, (1.1) 

with some functions Cij : M™ — > M. Compatibility of these equations is 
expressed by the following system of (nonlinear) differential equations: 



diCji^ — CijCjif -\- CjiCiP; Cji^Cii^, i ^ j ^ k ^ i, 



(1.2) 
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which thus spht off from eqs. for /. System (|1.1() . (|1.2|) is hyperbohc 

(see Sect. the fohowing data define a well-posed Goursat problem for 
this system and determine a conjugate net / uniquely: 

(Qi) values of / on the coordinate axes 23 j for 1 < i < m, i.e., m smooth 
curves /fs. with a common intersection point /(O); 

(Q2) values of Cij, Cji on the coordinate planes Tnj for all 1 < z < j < m, 
i.e., m{m — 1) smooth real- valued functions CjjI"Sjj of two variables. 

Alternative analytic description of conjugate nets. Given the func- 
tions Cij , define functions hi : M™ — > M as solutions of the system of differ- 
ential equations 

dihj = djhj , j- (1.3) 

Compatibility of this system is a consequence of H1.2|l . Define vectors Vi = 
h~^dif. It follows from (|l.lj) and (|1.H|) that these vectors satisfy the following 
differential equations: 

diVj = ^CjiVi, i^j. (1.4) 
Thus, defining the rotation coefficients as 

rij 

we end up with the following system: 

dif = hiVi , (1.6) 
diVj = PjiVi , 3, (1.7) 
dihj = hiPij , i^j. (1.8) 

Rotation coefficients satisfy a closed system of differential equations, which 
follow from eqs. (|1.2|) upon substitution (|1.5j) : 

diPkj = PkiPij , i^j^k^i. (1.9) 

Eqs. (jl.yj) . known as the Darboux system, can be regarded as compatibility 
conditions of the linear differential equations (jl.Tj) . 

Observe an important difference between two descriptions of conjugate 
nets: while the functions Cij describe the local geometry of a net, this is not 
the case for the rotation coefficients l3ij. Indeed, to define the latter, one 
needs to find hi as solutions of differential equations (|1.3() . 

Transformations of conjugate nets. The most general class of trans- 
formations of conjugate nets was introduced by Jonas and Eisenhart. 
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Definition 1.2 A pair of m- dimensional conjugate nets f,f^ : M 
is called a Jonas pair, if three vectors dif, dif^ and 6f = f'^ — f are co- 
planar at any point u G M"^ of the definition domain and for any 1 < i < m. 
The net f^ is called a Jonas transform of the net f . 

This definition yields that Jonas transformations are described by the fol- 
lowing (linear) differential equations: 

dif^ = aAf + bi{f+-f). (1.10) 

Of course, functions ai,bi : M"^ M have to satisfy (nonlinear) differential 
equations, which express the compatibility of eqs. (|1.1U|) with ()1.1|) : 

diUj = (oi - aj)cij + bi{aj - I), (1.11) 
dibj = c±bj + cp, - bjbi, (1.12) 
^j^tj = aiCij + bi{aj - 1). (1-13) 

Following data determine a Jonas transform f~^ of a given conjugate net / 
uniquely: 

(Ji) a point /+(0); 

(J2) values of Oj, 6j on the coordinate axes for 1 < i < m, i.e., 2m 
smooth real- valued functions fflifs. , feiTsj of one variable. 

Observe a remarkable conceptual similarity between Definitions 11.11 and 
11.21 Indeed, one can interpret the condition of Definition II . II as planarity of 
infinitesimal quadrilaterals {f{u), f{u + eiCi), f{u + e^Ci + ejCj), f{u + ejej)), 
while the condition of Definition 11.21 can be interpreted as planarity of in- 
finitesimally narrow quadrilaterals {f{u),f{u + eiei),f^{u + eiCi), f~^{u)). 



Classical formulation of the Jonas transformation. Our formula- 
tion of Jonas transformations is rather different from the classical one, which 
can be found, e.g., in |E2j . The latter is based on the formula 

whose data are: an additional solution (p : — > M of eq. (|1.1() . a parallel 
to / net g : M*" R^, and the function ip : R, associated to (f> in the 

same way as g is related to /. We now demonstrate how to identify these 
ingredients within our approach and how they are specified by the initial 
data (Ji,2)- 
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There follows from eqs. Hl.ll() - ()1.13() : 

a.(^)=cA + c..^-^- (1-15) 

The symmetry of the right-hand sides of eqs. (|1.15j) . (|1.12|) yields the exis- 
tence of the functions (j), 0"^ : M™" — > M such that 

did hj djd)^ 

^ = -, ^ = h, l<i<m. (1.16) 

These equations define (j), (j)^ uniquely up to respective constant factors, 
which can be fixed by requiring 0(0) = (j)~^{0) = 1. An easy computation 
based on eqs. ()1.15() . H1.12() shows that the functions cj},(j)~^ satisfy the fol- 
lowing equations: 

didj(j) = Cijdj(j) + Cjidi(j), (1-17) 
didj(j)+ = c±dj(t>-^ + cpi(l>+, (1.18) 

for all 1 < i 7^ J < m. Thus, a Jonas transformation yields some additional 
scalar solutions (f) and of the equations describing the nets / and 
respectively. It is clear that the solution cp is directly specified by the initial 
data (J2)- Indeed, these data yield the values (p along the coordinate axes, 
through integrating the first equations in l|1.16|) : these values determine the 
solution of eq. (jLlZf) uniquely. 

Further, introduce the quantities 

Then a direct computation based on eqs. (|1.1UI) . (|l.llj) - (|1.13|) . and (|1.16|1 

shows that the following equations hold: 

dig = aidif, (1.20) 

diip = aidi(p, (1.21) 

where 

a,- — 1 , , 

= ^ . (1.22) 

Thus, g is a parallel net to /, and V is an associated to (p function, in 
Eisenhart's terminology. Another computation leads to the relation 



diOj = Cij{ai - aj). 



(1.23) 
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The same argument as above shows that the data (J2) yield the values of 
(/>"*", and thus the values of ai, on the coordinate axes 23 j. This uniquely 
specifies the solutions of the compatible linear system ()1.23|) . This, in 
turn, allows for a unique determination of the solutions g, ip of eqs. (|1.2fl|) . 
(n^Tl) with the initial data g{0) = /+(0) - /(O) and ^(0) = 1 (here the 
data (Ji) enter into the construction). Thus, the classical formula (|1.14|) is 
recovered. 

One can iterate Jonas transformations and obtain a sequence /, f~^, 
(/+)+, ... , of conjugate nets. We will see that this can be interpreted as 
generating a conjugate net of dimension M = m + 1, with m continuous 
directions and one discrete direction. The most remarkable property of 
Jonas transformations is the following permutability theorem. 

Theorem 1.3 (Permutability of Jonas transformations) 

1) Let f be an m- dimensional conjugate net, and let f^^^ and /^^^ be its 
two Jonas transforms. Then there exists a two-parameter family of conjugate 
nets /^"^^^ that are Jonas transforms of both f^^'^ and /*-^-*. Corresponding 
points of the four conjugate nets f, f^^\ f^"^^ and f^^'^^ are co-planar. 

2) Let f be an m- dimensional conjugate net. Let f^^\ /^-^^ and f^^^ be its 
three Jonas transforms, and let three further conjugate nets f^^'^\ /^^^^ and 
/(13) 5e given such that is a simultaneous Jonas transform of and 
f^^\ Then there exists generically a unique conjugate net 

/(123) that IS a 

Jonas transform of f^^'^\ f^"^^^ and f^^^^. The net is uniquely defined 

by the condition that any its point is co-planar with the corresponding points 
of f^^\ Z^*"'^ (''iT'd f^''^'^ for any permutation (ijk) o/(123). 

The situations described in this theorem can be interpreted as conjugate 
nets of dimension M = m + 2, resp. M = m + 3, with m continuous and 
two (resp. three) discrete directions. 

The theory of discrete conjugate nets allows one to put all directions 
on an equal footing and to unify the theories of smooth nets and of their 
transformations. Moreover, we will see that both these theories may be seen 
as a continuum limit (in some precise sense) of the fully discrete theory, if 
the mesh sizes of all or some of the directions becomes infinitely small (see 
Fig. [21) . This way of thinking is the guiding idea and the philosophy of the 
discrete differential geometry. 

1.2 Orthogonal nets and their transformations 

Orthogonal nets. An important subclass of conjugate nets is fixed in the 
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following definition. 

Definition 1.4 A conjugate net f : MJ^ is called an m- dimensional 

0-net (orthogonal net) in M^, if there holds dif _L djf at any u € and 
for all 1 < i ^ j < m. Such a net is called an orthogonal coordinate system 
if m = N . 

Two-dimensional (m = 2) orthogonal nets are called 0-surfaces. An O- 
surface in is nothing but a surface parametrized along its curvature 
lines, or, otherwise said, parametrized so that both the first and the sec- 
ond fundamental forms are diagonal. Such a parametrization exists for a 
general surface in in the neighborhood of a non-umbilic point. Note 
that Definition 11.41 is dealing only with notions which are invariant under 
Mobius transformations. Thus orthogonal nets (as well as their Ribaucour 
transformations from Definition 11.51 below) belong to the Mobius differen- 
tial geometry. It will be important to preserve this symmetry group under 
discretization. 

For an analytic description of an orthogonal net / : M*" M^, introduce 
metric coefficients hi = \dif\ and (pairwise orthogonal) unit vectors Vi = 
h~^dif. Then there hold eqs. ()1.6() ~ ()1.9() . supplemented by the orthogonality 
constraint 

diPij + dj(3ji = -{diVi, djVj) , i^ j. (1.24) 

Indeed, eq. (|1.7j) holds since / is a conjugate net and Vj are orthonormal, 
and serves as a definition of rotation coefficients (3ji. Eq. (|1.8() is a direct 
consequence of (|1.6|) . (|1.7|) . To derive eq. (|1.24j) . one considers the identity 
didj{vi,Vj) = 0. So, a distinctive feature of orthogonal nets among general 
conjugate ones is that the rotation coefficients /3ji reflect the local geometry. 
In the same spirit, a solution to the system (\l.'d^ is given by the locally 
defined metric coefficients hi. 

Eq. (|1.24|) is an admissible constraint for the system (|l.Hj) ~ (|1.9j) . This 
has the following meaning: eq. ()1.24|) involves two independent variables i, j 
only, and it is therefore sensible to ask for it to be fulfilled on the coordinate 
plane 53 jj. One can easily check that if a solution to the system H1.6|) - H1.9|) 
fulfills eq. (|1.24|) on all coordinate planes for 1 < i < j < m, then it is 
fulfilled everywhere on M™. 

System (|1.6)) - (|1.9j) . (|1.24|) is not hyperbolic, therefore it is less clear what 
data form a well-posed problem for it. It can be shown (see Sect. 13. 4() that 
the following data can be used to determine an orthogonal net / uniquely: 

(Oi) values of / on the coordinate axes Hi for 1 < i < m, i.e., m smooth 
curves fl's^, intersecting pairwise orthogonally at /(O); 
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(O2) m{m — l)/2 smooth functions jij : M for all 1 < i < j < m, 

which have the meaning of jij = ^{diPij — djl3ji)\'Bij- 

Ribaucour transformations of orthogonal nets. 

Definition 1.5 A pair of m- dimensional orthogonal nets f, : M"* 
is called a Ribaucour pair, if the corresponding coordinate curves of f and 
envelope one-parameter families of circles, i.e. if at any u € and 
for any 1 < i < m the straight lines spanned by the vectors dif , dif^ at the 
corresponding points f , are interchanged by the reflection in the affine 
hyperplane orthogonal to 6f = /"*" — /, which interchanges f and f^. The 
net f^ is called a Ribaucour transform of f . 

To describe a Ribaucour transformation analyticahy, we write: 

5./+ = r,(a./-2^/^5/), (1.25) 

with some functions rj : — > M which obviously coincide (up to a sign) 
with the quotients of the corresponding metric coefficients, rf = (hf /hi)'^. 
Further, denote i = \5f\ and introduce the unit vector y = £~^6f, so that 
= f + £y. Then, in the case rj > 0, we find: 

vl = Vi-2{vi,y)y, diy = \ei{vf + Vi), (1.26) 

with the functions 6'i : ^ R defined as Oi = {ri - l)hi/l. Eqs. (|T3?)|l 
imply equations for the metric coefficients: 

h+ = hi + e^£, d^£ = -{vi,y){ht + hi). (1.27) 

(In the case rj < one has to change the sign of the quantities vf, hf in 
eqs. (|1.26p . ()1.27|) .) Compatibility of the system (|1.26|) yields that 6i have 
to satisfy certain differential equations: 

/3+ = Pij - 2{vi,y)0j, diOj = \ei{[3± + Pi,). (1.28) 

Following data determine a Ribaucour transform of a given orthogonal 
net / uniquely: 

(Ri) a point /+(0); 

(R2) values of 9i on the coordinate axes for 1 < i < m, i.e., m smooth 
functions ^jf-Bi of one variable. 
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According to the general philosophy, iterating Ribaucour transforma- 
tions can be interpreted as adding an additional (discrete) dimension to 
an orthogonal net. The situation arising by adding two or three discrete 
dimensions is described in the following fundamental theorem. 

Theorem 1.6 (Permutability of Ribaucour transformations) 

1) Let f be an m- dimensional orthogonal net, and let f^^^ and f^"^^ be 
its two Ribaucour transforms. Then there exists a one-parameter family of 
orthogonal nets Z^^^-* that are Ribaucour transforms of both f^^^ and f^'^\ 
Corresponding points of the four orthogonal nets f , f^^\ /^^^ and f^^'^^ are 
concircular. 

2) Let f be an m- dimensional orthogonal net. Let f^^\ /^^^ and f^^^ be 
its three Ribaucour transforms, and let three further orthogonal nets f^^'^\ 
/(23) and be given such that f^^^^ is a simultaneous Ribaucour transform 
of f^^^ and /'--'■*. Then there exists generically a unique orthogonal net f^^'^^^ 
that is a Ribaucour transform of f^^'^\ Z*-^^^ and Z^^^-*. The net /(^^3) is 
uniquely defined by the condition that the corresponding points of f^'^\ f^^^\ 
j(ik) g^^^ y(i23) ^j,g concircular for any permutation (ijk) o/(123). 

The theory of discrete orthogonal nets will unify the theories of smooth 
orthogonal nets and of their transformations. 

Mobius-geometric description of orthogonal nets. Since orthog- 
onal nets belong to the Mobius differential geometry, it is useful to describe 
them with the help of the corresponding apparatus (a sketch of which is 
given in Appendix to Chapter^. This has major conceptual and technical 
advantages. First, this description linearizes the invariance group of orthog- 
onal nets, i.e., the Mobius group of the sphere (which can be considered 
as a compactification of by a point at infinity). Further, using the Clif- 
ford algebra model of the Mobius differential geometry enables us to give a 
frame description of orthogonal nets, which turns out to be a key technical 
device. 

In this formalismus, the ambient space for points and hyperspheres of 
the conformal A^-sphere is the projectivized Minkowski space 

P(M^+i'i). The standard basis of the Minkowski space M^+^'i is de- 
noted by {ei, . . . ,e7v+2}- We denote also eo = ^(e7v+2 + ^N+i) and eoo = 
^(e7v+2 — ^N+i)- The points of the conformal A-sphere are elements of the 
projectivized light cone P(L^"'"-^), i.e., straight line generators of L^"*"^. The 
Euclidean space is identified, via 

vro : 9 / ^ / = / + eo + \f\^e^ G , (1.29) 
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with the section Qq of the cone by the affine hyperplane {.^o = 1}) 

where is the eo-component of ^ € in the basis {ei, . . . ,eN, eo, eoo}- 

Orientation preserving Euchdean motions of are represented as conju- 
gations by elements of ^Koo, the isotropy subgroup of in Spin^(A^ + 1, 1). 

It is not difficult to derive the following nice characterization of orthog- 
onal nets, due to Darboux (its second half follows directly from eq. H1.29() l. 



Theorem 1.7 A conjugate net f : — > M.'^ is orthogonal, if and only 
if l/l^ satisfies the same equation as f does, in other words, if the 

corresponding / = vro o / : is a conjugate net in W'^^^'^. 

As easily seen, metric coefficients hi = \dif \ satisfy also hi = \dif\. Hence, 
vectors Vi = h~^dif = Vi + 2{f,Vi)eao have the (Lorentz) length 1. Since 
(/) /) = 0) one readily finds that (/, Vi) = and hi = —{diVi, f). 

Theorem 1.8 (Spinor frame of an O-net) For an orthogonal net f : 
M^, i.e., for the corresponding conjugate net f : M"^ — > Qq^, there 
exists a function ^ : M™" "K^o ( called a frame of / ), such that 

f = ip-^eoi;, (1.30) 
Vi = 'il>~'^eiil>, 1 < i < m, (1-31) 

and satisfying the system of differential equations: 

diip = —GiipSi, Si = \diVi, 1 < i < m. (1.32) 

Note that for an orthogonal coordinate system (m = N) the frame ■i/' is 
uniquely determined at any point by the requirements (|1.3()|) and (|1.31|) . 

It is readily seen that the unit tangent vectors Vi satisfy eq. (|1.7|) with 
the same rotation coefficients (3ji = {diVj,Vi) = —{diVi,Vj). With the help of 
the frame ip we extend the set of vectors {vi : 1 < i < m} to an orthonormal 
basis {-Ofc : 1 < A: < iV} of T^-Q^ : 

vj^ = ip-'^e^ip, l<k <N. (1.33) 

Correspondingly, we extend the set of rotation coefficients according to the 
formula 



Pki = {diVk,Vi) = 



-{diVi,Vk) = -{diVi,tp ^e^ip), l<i<m, l<k<N. 
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Recall that we also have: 

hi = -{diVi,f) = -{diVi,'ip~^eoip), 1 < i < m. 

Thus, introducing vectors Si = 'ipSit(^~^, we have the following expansion 
with respect to the vectors e^: 

5, = il)Siip~^ = \%lj{diVi)tp'^ = ~^ l^ki^k + ^i^oo- (1-34) 

kjl^i 

It is easy to see that eq. ()1.9() still holds, if the range of the indices is 
extended to all pairwise distinct with 1 < i, j < m and 1 < A; < A^, and that 
the orthogonality constraint H1.24() can be now put as 

diPij + djiSji = - ^ f^kiPkj ■ (1-35) 
The system consisting of (|1.9() . (|1.35p carries the name of the Lame system. 

1.3 Moutard nets and their transformations 

We introduce now Moutard nets |Mouj without a geometric motivation, 
but they will play an extremely important role in the subsequent geometric 
considerations. 

Definition 1.9 A map / : — > is called an M-net (Moutard net), if 
it satisfies the Moutard differential equation 

did2f = qi2f (1.36) 

with some qi2 : R^ ^ R. 

On the first sight, the notion of M-net is not related to that of a conjugate 
net. In particular, there do not exist M-dimensional M-nets with M > 3. 
However, the relation is easily established: if : R^ — > R is any solution 
of the same Moutard equation H1.36|) (for instance, any component of the 
vector /), then y = f : R^ R^ is a special conjugate net in R^: 

did2y = -{d2\ogv)diy - {di\ogv)d2y. 

Such nets y are called conjugate nets with equal invariants., and they were in- 
tensively studied in the classical projective differential geometry |Tzij . Thus, 
in a projective space the class of M-nets coincides with the class of conjugate 
nets with equal invariants. 

Following data determine an M-net / uniquely: 
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(Ml) values of / on the coordinate axes 23i, ^2, i.e., two smooth curves 
fl'Bi with a common intersection point /(O); 

(M2) a smooth function qi2 : ^ M, having the meaning of the coefficient 
of the Moutard equation. 



Definition 1.10 Two M-nets f,f+ : R'" are called Moutard trans- 

forms of one another, if they satisfy (linear) differential equations 

dif^ + dif = pi{f+-f), (1.37) 
d2f^-d2f = P2(/+ + /), (1.38) 

with some functions pi,p2 : R'^ ^ R (or similar equations with all plus and 
minus signs interchanged). 

The functions pi, p2, specifying the Moutard transform, have to satisfy 

(nonhnear) differential equations, which express compatibility of eqs. ()1.37|) . 
(fT3H|) with eq. 

diP2 = d2Pi = -qi2+PiP2, (1.39) 

qt2 = -qi2 + 2pip2. (1.40) 

Following data determine a Moutard transform f"^ of a given M-net /: 

(MTi) a point /+(0) G M^; 

(MT2) values of the functions pi on the coordinate axes for i = 1,2, i.e., 
two smooth functions Pif^i of one variable. 

Classical formulation of the Moutard transformation. Due to the 

first equation in (|1.39|) . for any Moutard transformation there exists a func- 
tion : ^ M, unique up to a constant factor, such that 

Pi = P2 = ^- (1-41) 

The last equation in H1.39() implies that 6 satisfies eq. H1.36|) . This scalar 
solution of eq. (|1.36|) can be specified by its values on the coordinate axes 23j 
(i = 1,2), which are readily obtained from the data (MT2) by integrating 
the corresponding eqs. (|1.41|) . This establishes a bridge to the classical 
formulation of the Moutard transformation (see, e.g., |Mou[ ITzij ) . They 
used to specify a Moutard transform /"*" of the solution / of the Moutard 
equation ()1.36|) by an additional scalar solution 6 of this equation, via eqs. 
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()1.37|) . H1.38() with (|1.41|) . From these equations one can conclude that 
solves the Moutard equation with the transformed potential 

q+=q-2d^d2loge=^^^, e+ = i (1.42) 

In our formulation, the origin of the function 6 becomes clear: it comes from 
Pi: P2 by integrating the system (|1.41j) . Eq. (|1.42j) is then nothing but an 
equivalent form of eq. ()1.4U|) . 

Theorem 1.11 (Permutability of Moutard transformations) 

1 ) Let f be an M-net, and let f^^^ and f^'^^ be its two Moutard transforms. 
Then there exists a one-parameter family of M-nets /^^^^ that are Moutard 
transforms of both f^^^ and f^'^\ 

2) Let f be an M-net. Let f^^\ /^^^ and f^'^^ be its three Moutard trans- 
forms, and let three further M-nets f^^'^\ /^^^^ cLi^-d f^^^^ be given such that 

is a simultaneous Moutard transform of f^"^^ and f^^\ Then there exists 
generically a unique M-net that is a Moutard transform of f^^'^\ /*^^^^ 

and 

1.4 Asymptotic nets and their transformations 

Definition 1.12 A map f : ^ is called an A-surface (a surface 
parametrized along asymptotic lines), if at any point the vectors dff, 9|/ 
lie in the tangent plane to the surface f , spanned by dif, 82 f. 

Thus, the second fundamental form of an A-surface in is off-diagonal. 
Such a paramerization exists for a general surface with a negative Gaussian 
curvature. Definition 11.121 like the definition of conjugate nets, contains 
projectively invariant notions only. Therefore A-surfaces belong actually to 
the geometry of the three-dimensional projective space. In our presentation, 
however, we will use for convenience additional structures of (Euclidean 
structure and the cross-product). For the projective interpretation of these 
constructions, see |KoPj . A convenient description of A-surfaces is provided 
by the Lelieuvre representation which states: there exists a unique (up to 
sign) normal field n : — s- M'^ to the surface / such that 

dif = din X n, ^2/ = n x d2n. (1-43) 

Cross-differentiation of eq. (|1.43j) reveals that did2n x n = 0, that is, the 
Lelieuvre normal field satisfies the Moutard equation 



did2n = qi2n 



(1.44) 
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with some qi2 : — > M. This reasoning can be reversed: integration of eqs. 
(|1.4I-{|) with any solution n : — > of the Moutard equation generates an 
A-surface / : ^ M^. 

Theorem 1.13 A- surf aces in are in a one-to-one correspondence, via 
the Lelieuvre representation \1.4^^ , with M-nets in R^. 

An A-surface / is reconstructed uniquely (up to a translation) from its 
Lelieuvre normal field n. In its turn, an M-net n is uniquely determined by 
the initial data (Mi^2)j which we denote in this context by (Ai^2): 

(Ai) values of the Lelieuvre normal field on the coordinate axes Si, 232, 
i.e., two smooth curves nl-^. with a common intersection point n(0); 

(A2) a smooth function qi2 : ^ M, having the meaning of the coefficient 
of the Moutard equation for n. 

Definition 1.14 A pair of A- surfaces f,f^ : — > is called a Wein- 
garten pair, if, for any u € M'^, the line [f{u), f'^{u)\ is tangent to both 
surfaces f and f'^ at the corresponding points. The surface is called a 
Weingarten transform of the surface f . 

It can be demonstrated that the Lelieuvre normal fields of a Weingarten 
pair /, of A-surfaces satisfy (with the suitable choice of their signs) the 
following relation: 

/+ - / = n+ X n. (1.45) 

Differentiating the last equation and using the Lelieuvre formulas (|1.43|) for 
/ and for /+, one easily sees that the normal fields of a Weingarten pair are 
related by (linear) differential equations: 

5in^ + 5in = pi(n^ — n), (1-46) 
d2n~^ — = p2{n^-\-n), (1-47) 

with some functions pi,P2 '-^"^ IR- Thus: 

Theorem 1.15 The Lelieuvre normal fields n, n"*" of a Weingarten pair f , 
/■^ of A-surfaces are Moutard transforms of one another. 

A Weingarten transform f~^ of a given A-surface / is reconstructed from a 
Moutard transform of the Lelieuvre normal field n. The data necessary 
for this are the data (MTi^2) for n: 



(Wi) a point n+(0) G M^. 
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(W2) values of the functions pi on the coordinate axes for i = 1,2, i.e., 
two smooth functions pi l"^- of one variable. 

Following statement is a direct consequence of Theorem 11.111 

Theorem 1.16 (Permutability of Weingarten transformations) 

1) Let f be an A-surface, and let f^^^ and f^"^^ be its two Weingarten 
transforms. Then there exists a one-parameter family of A-surfaces /^^^^ 
that are Weingarten transforms of both f^^') and f^'^\ 

2) Let f be an A-surface. Let f^^\ /^^^ and f^^^ be its three Weingarten 
transforms, and let three further A-surfaces f^^'^\ /^^^^ cL^d f^^^^ be given 
such that Z^*-'^ is a simultaneous Weingarten transform of /^*) and f^^\ 
Then there exists generically a unique A-surface that is a Weingarten 
transform of f^^'^\ /^^^^ and f^^^^. The net /(^^^) is uniquely defined by 
the condition that any its point lies in the tangent planes to f^^'^\ f^'^^^ and 
/^^^•* at the corresponding points. 

1.5 Surfaces with constant negative Gaussian 
curvature and their transformations 

Up to now, we discussed special classes of coordinate systems in space, or 
special parametrizations of a general surface. Now, we turn to the discussion 
of several special classes of surfaces. The distinctive feature of these classes 
is the existence of transformations with certain permutability properties. 

Definition 1.17 An A-surface / : — > is called a K-surface (or a 
pseudospheric surface^, if its Gaussian curvature K is constant, i.e., does 
not depend on u . 

K-surfaces constitute one of the most prominent examples of integrabil- 
ity in the differential geometry. One of the approaches to their analytical 
study is based on the investigation of the angle between asymptotic lines 
which is governed by the famous sine-Gordon equation. This approach was 
transferred to the discrete setting in iBobPlj . see also a presentation in 
|BobMaST] based on the notion of consistency. In the present paper, we take 
an alternative route, based on the study of the Gauss map of K-surfaces. 
Following are the classical characterization results. 

Theorem 1.18 

1) An A-surface / : ^ is a K-surface, if and only if the functions 
\dif\ = Oi depend on Ui only (i = 1,2). 
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2) The Lelieuvre normal field n : M? ^ M.^ of a K-surface with K = —1 
takes values in the sphere C M^, thus coinciding with the Gauss map. 
Conversely, any M-net in the unit sphere is the Gauss map and the 
Lelieuvre normal field of a K-surface with K = —1. There holds: \din\ = at 
for i = 1,2, with the same functions ai = C(i{ui) as in 1). 

Thus, K-surfaces are in a one-to-one correspondence with M-nets in 
(otherwise said, with Lorentz-harmonic S^-valued functions). It is important 
to observe that the coefficient qi2 of the Moutard equation satisfied by a a 
Lorentz-harmonic S^-valued function completely determined 

by n, more precisely, by its first order derivatives: 



Therefore, following data determine the Gauss map n of a K-surface /: 

(K) values of the Gauss map on the coordinate axes 23i, 232, i-e., two smooth 
curves ?^|''Bj in intersecting at a point n(0). 

The K-surface / is reconstructed from n uniquely, up to a translation, via 
formulas (|1.43j) . 

Historically the first class of surface transformations with remarkable 
permutability properties was introduced by Backlund. 

Definition 1.19 (Backlund transformation) A Weingarten pair of K- 
surfaces /, /"^ : ^ forms a Backlund pair, if the distance \f~^ — f \ is 
constant, i.e., does not depend on n S M^. 

The Gauss maps n, n"*" of a Backlund pair of K-surfaces /, are related 
by the Moutard transformation H1.46() . ()1.47|) . From these equations there 
follows easily that for a Backlund pair the quantity (n,ra^) is constant; thus, 
the intersection angle of the tangent planes at the corresponding points of 
a Backlund pair is constant. Moreover, eq. (|1.45j) yields that this constant 
angle is related to the constant distance between / and /+ via 



The fact that n,n+ € allows one to express the coefficients pi, p2 in eqs. 
(|1.46|) . (|1.47|) in terms of the solutions themselves: 



qi2 = {did2n,n) 



{din, d2n). 



(1.48) 



|/+-/|2 = l-(n,n+)2. 



{din,n'^) — {n,din'^) 



{din, n^) 



(1.49) 



Pi = 



2 - 2{n,n+) 
{n,d2n~^) — {d2n,n~^) 
2 + 2{n,n+) 



1 — (n, n+) 
-((92^, n+) 
1 + (n, n+) 



P2 = 



(1.50) 
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With these expressions, eqs. (|1.46p . H1.47() become a compatible system of 
first order differential equations for n"*", therefore the following data deter- 
mine a Backlund transform /"*" of the given K-surface / uniquely: 

(B) a point n+(0) e S^. 

Permutability of Backlund transformations is due to Bianchi: 

Theorem 1.20 (Permutability of Backlund transformations) 

Let f be a K-surface, and let f^^^ and /^^^ be its two Backlund trans- 
forms. Then there exists a unique K-surface /^^^^ which is simultaneously a 
Backlund transform of f^^^ and of f^'^\ The fourth surface f^^'^^ is uniquely 



defined by the properties = |/(2)_/[ and |/(i2)_y(2)| ^ 



We will see how the theory of discrete K-surfaces unifies the theories of 
smooth K-surfaces and of their Backlund transformations. 

1.6 Isothermic surfaces and their transformations 

Classically, theory of isothermic surfaces and their transformations was con- 
sidered as one of the highest achievements of the local differential geometry. 

Definition 1.21 An 0-surface / : — > is called an I-surface (isother- 
mic surface), if its first fundamental form is conformal, possibly upon a 
re-parametrization Ui — > (pi{ui) {i = 1,2) of the dependent variables, i.e., 
if at any point u G of the definition domain there holds |9i/p/|02/P = 
ai{ui)/a2{u2). 

In other words, isothermic surfaces are characterized by the relations 8182 f 
span(5i/, ^2/) and 



with some s : — > M+ and with the functions Ui depending on Ui only 
(i = 1,2). These conditions may be equivalently represented as 



did2f = {82 log s)dif + (5i log 5)^2/, {dif, d2f) = 0. (1.52) 



The following property of isothermic surfaces can actually serve as their 
another characterization. 




{dif, d2f) = 0, \dif\^ = ais\ \d2f\^ = a2s\ 



(1.51) 
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Theorem 1.22 (Dual I-surface) Let f -.M."^ —>■ be an isothermic sur- 
face. Then the M.^ -valued one-form df* defined by 



is closed. The surface f*:M.'^—>- M^, defined (up to a translation) by the 
integration of this one-form, is isothermic, with 

{dif*,d2f*) = 0, \dif*\^ = ai.s'^ \d2f*\^ = a2S-\ (1.54) 

The surface f* is called dual to the surface f, or the Christoffel transform 
of the surface f. 

Another important class of transformations of isothermic surfaces build the 
Darboux transformations. 

Definition 1.23 (Darboux transformation) A Ribaucour transform 
: — of a given isothermic surface / : — > is called a Dar- 
boux transform, if its first fundamental form is likewise conformal, possibly 
upon a reparametrization of the dependent variables, i.e., if at any point u £ 
of the definition domain there holds |5i/^p/|52/'^P = ai{ui) / a2{u2) ■ 

Introduce the corresponding function : ^ for the surface 
/+, and denote r = s+/s : ^ M+. Thus, |9i/+|V|5i/P = (s^/s? = 

for i = 1,2. Comparing this with the definition (|1.25|) of Ribaucour 
transformations we see that one of the two possibilities holds: 

(i) ri = r2 = r, or (ii) ri = — r2 = —r. 

It can be demonstrated that in the case (i) the surface f~^ is with necessity 
a Mobius transformation of /; we will not consider this trivial case further. 
In the case (ii) one gets proper Darboux transformations. An important 
property of the Darboux transformations is the following: the quantity 

c= — = — (1.55) 

ss^ ss^ 

is constant, i.e., does not depend on n G M^. It is called a parameter of the 
Darboux transformation. Following data determine a Darboux transform 
of a given isothermic surface / uniquely: 

(Di) a point /+(0); 
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(D2) a real number c, designated to be the constant H1.55() . 

As usual, we regard iterating a Darboux transformation as adding a 
third (discrete) dimension to a two-dimensional isothermic net. The main 
classical result on Darboux transformations is the following theorem, which 
assures that one can add several discrete dimensions in a consistent way. 

Theorem 1.24 (Permutability of Darboux transformations) 

Let f be an isothermic surface, and let f^^^ and f^'^^ be its two Darboux 
transforms, with parameters Ci and C2, respectively. Then there exists a 
unique isothermic surface Z*-^^^ which is simultaneously a Darboux transform 
of f^^^ with the parameter C2 and a Darboux transform of f^'^^ with the 
parameter ci. The surface Z*-^^^ is uniquely defined by the condition that 
the corresponding points of the four isothermic surfaces are concircular, and 
have a constant cross-ratio 

q{fj^'\f^''\f^'^) = "-^. 

Remark. The real cross-ratio of four concircular points a, b,c,d G 
may be defined as 

q{a, b, c, d) = {a- b){b - c)"^(c - d){d - a)"\ (1.56) 

where the points are interpreted as elements of the Clifford algebra C£(M^). 
In more down-to-earth terms, since the four points are co-planar, we may 
identify the plane where they lie with the complex plane, and interpret in 
the above formula the symbols a, b, c, d as complex numbers. 

The theory of discrete isothermic surfaces unifies the theories of smooth 
isothermic surfaces and of their Darboux transformations. 

Mobius-geometric characterization of isothermic surfaces and 
their Darboux transformations. It is easily checked that conditions 
(|1.5H) are invariant with respect to affine transformations of M'^ , as well as 
with respect to the inversion /—>//(/,/). In other words, the notion of 
isothermic surfaces belongs to the Mobius differential geometry. The same 
holds for their Darboux transformations. Therefore, it is useful to charac- 
terize these notions within the Mobius-geometric formalism. (However, the 
notion of the dual surface, or Christoffel transformation, is essentially based 
on the Euclidean structure of the ambient space M^.) 

To find such a characterization, note first of all that eqs. p.5'2j) are 
equivalent to 

did2f = {d2 log s)dif + ((9i log 5)^2/ 
for the image / : ^ of / in the quadric C L^+i. 
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Theorem 1.25 The lift s = s ^/ : — s- L^"^^ of an isothermic surface 
/ : — to the light cone ofM.^^^'^ satisfies the Moutard equation 

did2S = qi2S, (1.57) 

with qi2 = sdid2{s~^) . 

Conversely, given an M-net s : — > L^"*"^ in the light cone, define 
s : ^ M+ and f ■.R'^ ^ by 

s = s-^{f + eo + 1/peoo) 

(so that s~^ is the eo-component, and s~^ f is the R.^ -part of s in the basis 
ei, . . . , e/vT, eo, Goo). Then f is an isothermic surface. 

Note that the functions ai = {diS,dis) {i = 1,2) depend on Ui only and 
coincide with the namesake functions from the definition ()1.51() . 

Thus, we see that I-surfaces are in a one-to-one correspondence with 
M-nets in L^^^, i.e., with Lorentz-harmonic h^^^ -valued functions. 

Let us address the problem of minimal data which determine an isother- 
mic surface (i.e., an M-net in L^"^^) uniquely. Guided by an analogy with 
the case of K-surfaces, one is tempted to think that two arbitrary curves sf-B. 
in L^^^ would be such data. However, as a consequence of the fact that 
now we are dealing with the light cone L^"*"^ = {(s, s) = 0} rather than with 
the sphere = {(n, n) = 1} as a quadric where M-nets live, we cannot find 
an expression for qi2 in terms of the first derivatives of s anymore; instead, 
one has: 

_ {dfs,d2s) _ {d^s,dis) 
{dis,dis) {d2S,d2s)' 

This shows that the coordinate curves s [3. are not arbitrary but rather 
subject to certain further conditions. We leave the question on correct initial 
data for an isothermic surface open. 

Darboux pairs of isothermic surfaces are characterized in terms of their 
lifts as follows. 

Theorem 1.26 The lifts s, s+ : ^ L^+^ of a Darboux pair of isothermic 
surfaces f,f^ : — > are related by a Moutard transformation, i.e., 
there exist two functions pi,p2 : R'^ ^ R such that 

dis^ -\- dis = pi{s~^ — s), d2S~^ — 828 = p2{s~^ + s). (1.58) 

Conversely, for an M-net s in L^"^-^, any Moutard transform s+ with values 
in L^"*"^ is a lift of a Darboux transform /+ of the isothermic surface f. 
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Note that the quantity (s, s~^) is constant (does not depend on u G M^), and 
is related to the parameter c of the Darboux transformation: (s, s^) = — c/2. 
The formulas 

_ {s,dis+)-{dis,s+) _ {d,s,s+) 

make it apparent that a Moutard transform is completely determined by 
prescribing its value s"^(0) at one point. Indeed, eqs. (|1.58|) with coefficients 
()1.59() form a compatible system of first order differential equations for s"*" : 
course, data (Di^2) are encoded in s~''(0) in a straightforward 

manner. 

We summarize the considerations of these chapter in the following table: 

0-net / in ^ conjugate net / in ~ P(L^+^), 

A-net / in ^ M-net n in M^, 

K-net / in ^ M-net n in S^, 

I-net / in ^ M-net s in L^+^ 

The next chapter will be devoted to discretizing all these relations. 



Appendix: Mobius-geometric formalism 

The classical model of the A^-dimensional Mobius geometry, which allows 
for a linear representation of Mobius transformations, lives in the Min- 
kowski space M^+^'^, i.e., in an (A^ -|- 2)-dimensional space with the basis 
{ei, . . . ,e]\f+2}, equipped with the Lorentz scalar product in which are 
pairwise orthogonal and 




1, l<i<N + l, 
-1, i = N + 2. 



(Although we use the same symbol for the Lorentz scalar product in R^"*"^'^ 
and for the Euclidean scalar products in R-'^ and in R'^"'"^, its concrete 
meaning should be always clear from the context.) 

Points. The space of points in the Mobius geometry is P(L^"^^) - the 
space of straight line generators of the light cone 



L^+i = {eGR^+''^ : (^,0=0} 



(1.60) 
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The sphere § C (where we regard as spanned by ei, . . . , e^v+i) 

is identified with a section of L^^"*^ by the affine hyperplane eAr+2) = — 1: 

^ = {e e L^+i : = 1}, 

which is a copy of shifted by eAr+2: 

7n:S^3y ^ y = ysph = y + eiv+2 eQ^. (1.61) 

Similarly, the Euclidean space may be identified with the section of 
l^N+i |-|y ^j^g affine hyperplane (^,e7v+2 — ^N+i) = —1, 

:^ = G L^+i : ^iv+i + = 1}, 
(Euclidean metric d^f + . . . + d£^j^ being induced from the ambient M^^"*^'^): 

vTo : 9 a; 1-^ x = xeuc = 2; + ^(1 - gat+i + ^(1 + [xp) eAr+2 

= x + eo + |x|2eoo gQo • (1-62) 

Here the following notations are introduced: 

eo = + gn+i) and eoo = - ^n+i)- 

Thus, the space is modelled as a paraboloid in an {N + l)-dimensional 
affine subspace through eo spanned by ei, . . . , e^r, eoo- An important prop- 
erty of the Euclidean identification H1.62() is: 

X2) = — — X2p, Vxi,X2GlR^. (1.63) 

Note that the correspondence between Q-'^ and Qg^ along the straight line 
generators of L''^"'"^ induces the stereographic projection a : — > , 

1 + Ixp 1 + 

In particular, the generators of L^"*"^ through the points eo and eoo corre- 
spond to the north pole yo = eN+i and the south pole yoo = —^N+i on 
and to the zero and the point at infinity in M^, respectively. 

Spheres. A hypersphere S in the conformal A^-sphere is the (non- 
empty) intersection of P(L'^+^) with a projectivized hyperplane. Thus, S 
can be put into a correspondence with the point S G P(M^+-'^'^) polar to the 
above mentioned hyperplane with respect to the light cone. This point is 
space-like, i.e., any its representative s £ M^+^'-^ has {s,s) > 0. There are 



1.6. ISOTHERMIC SURFACES 



35 



various choices of this representative which have nice geometric interpreta- 
tions. 

Fixing the ejv+2-component of s leads to the choice 

s = ssph = s + eAr+2, with s € M^+\ {s,s)>l. (1.64) 

This is related to a description of the hypersphere S C §^ as the intersection 
of with the hyperplane (s,y) = 1 in M^+^. Indeed, the latter equation 
is equivalent to {s,y) = for y from eq. The point s lies outside 

of and S C S'^ is the contact set of S''^ with the tangent cone to 
with the apex s. Also, 5 C is the intersection of with the orthogonal 
A^-sphere in R^+^ with the center s and the radius p = {{s,s) — l)-*^/^. See 
Fig.E2l 

Fixing the eo-component of s leads to the choice 

s = seuc = c -|- eo + (|cp — r^)eoo, where c € M^, (1.65) 

related to the Euclidean description of the hypersphere S C M^. Indeed, 
{s,x) = for X from eq. (|1.62() is equivalent to \x — cp = r^. This is the 
equation for points x of the sphere S with the Euclidean center c and the 
Euclidean radius r. 

Still another choice is to fix the Lorentz norm of s: 

s = SMob = ±-(s + eiv+2) = ±-(c + eo + ([cp - r^)eoo) € L^+\ (1.66) 
p r 

where for any k > the following quadric is introduced: 

L^+i = U G M^+i'i : (e,0 = '^'}- (1-67) 

Actually, this represents oriented hyperspheres, each choice of the sign it cor- 
responding to one of the two possible orientations of a given hypersphere. 
For any two (oriented) hyperspheres ^i, 5*2 the scalar product of their rep- 
resentatives SMob is a Mobius invariant: if k = 1, then 

1 1 

(si,S2) = ((si,S2) - l) = 7, {rf+rl - |ci - C2p), 

PiP2 2rir2 

is the cosine of the intersection angle of 5i, 5*2, if they intersect, and the 
inversive distance between Si, S2, otherwise. 

Transformations. Elements of the group 0~''(A^ -|- 1,1), i.e., Lorentz 
transformations preserving the time-like direction, induce Mobius transfor- 
mations on ~ S^. Therefore we identify 0^(A^ -|- 1, 1) with the group 
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M(A'^) of Mobius transformations of S . Similarly, Lorentz transformations 
which fix Coo, induce Euclidean motions on Q-^ ~ M^, and therefore we 
identify the corresponding isotropy subgroup 0^(iV + 1, 1) with E(A^), the 
group of Euclidean motions of . 

It is convenient to work with spinor representations of these groups. Re- 
call that the Clifford algebra Qi{N+l, 1) is an algebra over M with generators 
ei, . . . , eAr+2 € M^^'-'^'^ subject to the relation 

^r? + r/C = -2(C,r?)l = -2(^,7?), V^,?? e M^+^'^ 

This implies that = -(^,^), therefore any vector ^ G \ L^+^ has 

an inverse = — ^/(^, ^). The multiplicative group generated by invertible 
vectors is called the Clifford group. We need its subgroup generated by unit 
space-like vectors: 

g = Pin+(iV + 1, 1) = {V; = 6 • • -^n : 4' = -1}, 

and its subgroup generated by vectors orthogonal to Boo'. 

Soo = Pmi(iV + 1, 1) = = ^1 . . . : = -1, iCi, eoo) = 0}. 

These groups act on M^^+^'-^ by twisted conjugations: A^{r]) = {—ly-tl^^^rjip. 
In particular, for a vector ^ with = — 1 one has: 
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which is the reflection in the hyperplane orthogonal to ^. Thus, S is gener- 
ated by reflections, while 9oo is generated by reflections which fix eoo, and 
therefore leave Qq invariant. Actually, S is a double cover of 0+(A''-|-l, 1) ~ 
M(A^), while Soo is a double cover of 0+(A^ + 1, l)oo ^ E(A^). Orientation 
preserving transformations from S, Soo form the subgroups 

:K = Spin+(iV + l,l), J{oo = Spin+(iV + l,l), 

which are singled out by the condition that the number n of vectors in 
the multiplicative representation of their elements ip = ■ ■ ■ ^s. even. The 
Lie algebras of the Lie groups and Ji^o consist of bivectors: 

[) = spin(A/" + 1, 1) = spanjeie^- : i,j e {0,1, ... ,N,oo}, i ^ 
[)oo = spinoo(A/"-|- 1,1) = spanjejej : i,je{l,...,N,oo}, 
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Chapter 2 



Discrete differential 
geometry 

For functions on Z*^, we define translation and difference operators in a 
standard manner: 

(T,/)(n) = f{u + a), {6if){u) = f{u + a) - f{u), 

where Cj is the i-th coordinate vector of Z*^. We use the same notation for 
(discrete) s-dimensional coordinate planes, 

= {li G Z*^ : = for i / ii, . . . , is}, 

as in the continuous case. 

2.1 Discrete conjugate nets 

The following definition is due to Sauer |Saj for M = 2, and to Doliwa and 
Santini |Do?^ for general M. See [Doll lDo2l lDo3l IdHKI IDoSMl IMPS) for 

further relevant developments. 

Definition 2.1 A map f : Z*^ is called an AI -dimensional Q-net 

(quadrilateral net, or discrete conjugate net) inW^ , if any of its elementary 
quadrilaterals is planar, i.e., if at any u € Z^^ and for all pairs 1 < i 7^ j < 
M the four points f, Tif , Tjf, andT-iTjf are co-planar. 

Note that this definition actually belongs to the projective geometry, as it 
should. To understand what restrictions does this condition impose, we 
consider various values of M. 
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M=2: discrete surface parametrized by conjugate lines. Suppose 
two coordinate lines, f\-B^ and /f^j > on a Q-surface / : — > are given. 
To extend the surface into the quadrant Z^, say, one proceeds by induction 
whose step consists of choosing /12 in the plane spanned by /, /i and /2, 
provided the latter three points are known (here we write /j, fij for Tj/, 
TiTjf, etc.). The planarity condition is equivalent to the relation 

^l^2f = C2lSlf + €1262/. 

So, one has two free real parameters C21, C12 on each such step. It is con- 
venient to think of these parameters as attached to the elementary square 
(n, n + ei,u + ei + 62, u + 62) of the lattice Z^. Thus, one can define a 
Q-surface / by prescribing its two coordinate lines / , / ; two 
real- valued functions ci2,C2i defined on all elementary squares of Z^. 

Actually, the combinatorics of Q-surfaces may well be more complicated 
than that of Z^. Indeed, Definition 12.11 can be literally extended to maps 
/ : V{1)) M^, where V{1)) is the set of vertices of an arbitrary quad- 
graph T). A quad-graph is a strongly regular cell decomposition of a surface 
with all quadrilateral faces. We will later need also the notation E{D) and 
F(D) for the sets of edges, resp. faces, of a quad-graph D. As we will show, 
the integrable nature (multi-dimensional consistency) of the Q-nets gives an 
opportunity to realize D as a surface in some Z*^ and to work only on this 
larger (but simpler) definition domain. 

M=3: basic 3D system. Suppose that three coordinate surfaces of 
a three-dimensional Q-net / are given, that is, / f^^j > / \'B23 ^iid / f S13 • 
Of course, each one of them is a Q-surface. To extend the net into the 
octant Zi^, one proceeds by induction whose step consists of determining 
/123, provided /, /j and fij are known for all 1 < i 7^ j < 3. The point 
/123 has to lie in three planes TiUjk (1 < i < 3), where Tjlljfc is the plane 
passing through three points {fi, fij, fik)- This condition determines /123 
uniquely. Indeed, all three planes rin23, r2ni3 and T3ni2 belong to the 
three-dimensional affine space through the point / spanned by the vectors 
^if (1 ^ i ^ 3), and therefore generically these planes intersect at exactly 
one point. An elementary construction step of a three-dimensional Q-net 
out of its three coordinate surfaces, i.e., finding the eighth vertex of an 
elementary hexahedron out of the known seven vertices, is symbolically rep- 
resented on Fig. 12.11 This is the picture we have in mind when thinking 
(and speaking) about three-dimensional systems. Of course, one can also 
give an analytic formulation of this picture. This is done as follows. The 
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characteristic property of a Q-net is: 



5i6jf = Cji6if + Cij6jf, i / j. 



(2.1) 



Here, as before, functions Cij, Cji, as well as equation H2.1|) itself, are thought 
of as defined on elementary squares of parallel to the coordinate plane 
"Bij. Six such equations, attached to six facets of an elementary cube of 
Z^, form the three-dimensional system encoded on Fig. 12.11 Here the num- 
bers {cjk} on the facets adjacent to / are considered as known, while the 
numbers {riCjk} on the facets adjacent to /123 are uniquely defined by the 
compatibility of eqs. (|2.1|) on all six facets. In other words, it is required 
that 6i{5j5kf) does not depend on the permutation {i,j,k) of the indices 
(1,2,3). This compatibility condition gives: 

SiCjk = {TkCij)cjk + {TkCji)cik - iTiCjk)cik, i ^ j ^ k i. (2.2) 

Note that equations for Cjk split off from equations for /; they constitute a 
system of 6 (linear) equations for 6 unknown variables TiCjk in terms of the 
known ones Cj^. The resulting map {cjk} 1— > {TiCjk} is birational. Sometimes 
it is this map that is considered as the fundamental 3D system encoded on 



M > 4- consistency. Turning to the case M > 4, we see that one 
can prescribe all two-dimensional coordinate surfaces of a Q-net, i.e., /t-Sj^ 
for all 1 < i < j < M. Indeed, these data are clearly independent, and 
one can construct the whole net from them. In doing so, one proceeds by 
induction, again. The inductive step is essentially three-dimensional and 
consists of determining provided /, fi and fij are known. However, 



Fig.im 





Figure 2.1: 3D system on an elementary cube 
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this inductive process works, only if one does not encounter contradictions. 
To see the possible source of contradictions, consider in detail the case of 
M = 4; higher dimensions do not add anything new. From /, /j and fij 
one determines all fijk uniquely. After that, one has, in principle, four 
different ways to determine /1234, from four 3-dimensional cubes adjacent 
to this point; see Fig. 12.21 A remarkable property of Q-nets is that these 
four values for /1234 automatically coincide. We call this property the 4D 
consistency. 

Definition 2.2 A 3D system is called 4D consistent, if it can be imposed 
on all three-dimensional facets of an elementary hypercube of I?'. 




Figure 2.2: 4D consistency of 3D systems 



The following fundamental theorem is due to Doliwa and Santini |DoSlj . 

Theorem 2.3 The 3D system governing Q-nets is 4^D- consistent. 

Proof. Actually, the statement of the theorem is about the properties of the 
map {cjk} ^ {TjCjfc}. For such maps with the fields on 2D plaquettes (here 
each plaquette carries two fields) the 4D consistency means that the two 
values Ti{TjCki) and Tj{TiCM) coincide for any permutation {i,j,k,i) of the 
indices (1, 2, 3, 4). However, an algebraic proof of this claim could be hardly 
performed without help of a computer system for symbolic computations. 
A geometric approach, dealing with this system augmented by the fields 
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/ € R on the vertices, allows us to perform a conceptual proof, free of 
computations. Indeed, the map {cjk} i— > {najk} does not depend on the 
dimension N of the space where / lies. Therefore, we are free to assume that 
iV > 4. This will enable us to use geometric "general position" arguments. 
In the construction above, the four values in question are 

/l234 = TiT'2^34, H riT3n24 H riT4n23 , 

and three other ones obtained by cyclic shifts of indices. In general position, 
the affine space 111234 through / spanned by 6if (1 < i < 4), is four- 
dimensional. We prove that all six planes TiTjUki in this four-dimensional 
space intersect generically at exactly one point (which will be then /i234). It 
is easy to understand that the plane TiTjUi^i is the intersection of two three- 
dimensional subspaces TiUjki and TjHiki- Here, of course, the subspace 
TiUjki is the one through the four points (fi, fij, fik, he), or, equivalently, 
the one through fi spanned by 5jfi, 5k fi and difi. Now the intersection in 
question can be alternatively described as the intersection of the four three- 
dimensional subspaces rin234, T2ni34, T3ni24 and r4ni23 of one and the 
same four-dimensional space 111234. This intersection consists in the generic 
case of exactly one point. □ 

The M-dimensional consistency for M > 4 is defined and proved analo- 
gously. Actually, it follows from the 4-dimensional consistency. 

On the level of formulas we have for M > 4 the system (|2.1|) . (|2.2|) . 
where now all indices i,j,k vary between 1 and M. This system consists 
of interrelated three-dimensional building blocks: for any triple of pairwise 
different indices {i,j,k) the equations involving these indices only form a 
closed subset. The M-dimensional consistency of this system means that all 
three-dimensional building blocks can be imposed without contradictions. 
A set of initial data which determines a solution of the system ()2.1() . (|2.2() . 
consists of 

(Q^) values of / on the coordinate axes for 1 < i < M; 

{Q2) values of Cij, cji on all elementary squares of the coordinate plane "Bij, 
for 1 < i < j < M. 

Quadratic reduction of Q-nets. An important observation made 
by Doliwa |.Do2^ is that quadrilateral nets can be consistently restricted to 
an arbitrary quadric in M^. The importance of this resides on the fact 
that many of the geometrically relevant nets turn out to be reductions of 
quadrilateral nets to some quadrics or to intersections of quadrics. The 
quadratic reduction is based on the following fundamental claim: 
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Theorem 2.4 // seven points f , fi and fij (1 < i < j < 3) of an elemen- 
tary hexahedron of a quadrilateral net belong to a quadric Q C M^, then so 
does the eighth point /i23- 

Proof. A deep reason for this is the fact, well known in the projective geom- 
etry of the 19-th century: for seven points of a three-dimensional projective 
space in a general position, one can find the eighth associated point which 
belongs to any quadric through the original seven points (these quadrics 
form a two-parameter linear family, which contains Q). Under conditions 
of Theorem 12. 4| we have three (degenerate) quadrics through the original 
seven points: the pairs of planes Hjk U Tjlljfc for i = 1,2,3. Clearly, their 
eighth intersection point is /123 = Tin23 n T2n3i n r3ni2, and this has to be 
the associated point. □ 

Corollary. If the coordinate surfaces fl'Sij of a Q-net f : M.^ lie 

in a quadric Q, then so does the whole of f . 

Alternative analytic description of Q-nets. In a complete analogy 
with the smooth case, one can give a (non-local) description of discrete con- 
jugate nets, with somewhat simpler equations. Given the plaquette functions 
Cij , define quantities hj , attached to the edges parallel to the j'-th coordinate 
axes, as solutions of the system of difference equations 



whose compatibility is assured by eq. (|2.2|) . Introduce vectors vj = h~^6jf, 
attached to the same edges. There follows from (|2.1() and 1)2. 3|) that these 
vectors satisfy the following difference equations: 



Define the discrete rotation coefficients (attached to plaquettes) as 



6ihj = 



^ijhj , 



(2.3) 




(2.4) 




(2.5) 



Then we end up with the following system: 




Sif 

5iVj 



hiVi , 



(2.6) 
(2.7) 
(2.8) 
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Rotation coefficients satisfy a closed system of difference equations {discrete 
Darhoux system), whicli follows from eqs. (|2.2|) upon substitution (|2.5|) . or 
otherwise can be derived as compatibility conditions of the linear difference 
equations (|2.7j) : 

^Aj = {TjPki)Pij , iy^ j ^i. (2.9) 

This system, which is considerably simpler than H2.2|) . was first derived in 
BogK| without any relation to geometry. A geometric interpretation in 
terms of Q-nets was found in |DoSlj . It should be mentioned that eqs. ()2.8|1 
and (|2.9|1 are implicit, but can be easily solved for the shifted variables, 
resulting in 

Tihj = 

The last formula defines an explicit rational 3D map {Pkj} '-^ {Tif^kj}- Like 
the map {ckj} '-^ {TiCkj} for the local plaquette coefficients, the map ()2.11|) 
is 4D-consistent, but now this can be checked by an easy computation "by 
hands" . 

Transformations of Q-nets. A natural generalization of Definition 
11.21 would be the following one. 

Definition 2.5 A pair of m- dimensional Q-nets f,f^ : TT^ is called 

a Jonas pair, if four points f, Tif , f'^ and Tif^ are co-planar at any point 
u € and for any 1 < i < m. The net f^ is called a Jonas transform of 
the net f . 

But actually this relation can be re-phrased as follows: set F{u,0) = f{u) 
and 1) = then F : x {0, 1} ^ is an M-dimensional 

Q-net, where M = m + \. Thus, in the discrete case there is no difference 
between conjugate nets and their Jonas transformations. The situation of 
Definition 12.51 is governed by the equation 

5^f^ = ai5if + h{f+-f), (2.12) 

where coefficients Oj, hi are nothing but ai = 1 + cmi , h = oim- These 
coefficients are naturally attached to elementary squares of parallel to 
the coordinate plane "Bj^f . It is also convenient to think of them as attached 
to edges of Z"^ parallel to (to which the corresponding "vertical" squares 



1 Pij(3ji 
Pkj + PkiPij 
1 Pijf3ji 



(2.10) 
(2.11) 
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are adjacent). Equations of the system (|2.2|) with one of the indices equal 
to M give: 

5iaj = {Tjbi){aj - 1) + {TjUi - Tiaj)cij, (2-13) 
^^bJ = c±bj + cp, - {Tib,)bi, (2.14) 
= {Tjai)cij + {Tjbi){aj - I). (2.15) 

Following data are needed to specify a Jonas transform of a given 
m-dimensional Q-net /: 

(Jf ) value of /+(0); 

(J2^) values of Oj, bi on all edges of the respective coordinate axis Sj, for 
1 < i < m. 

Alternative description of discrete Jonas transformations. We 

give here a discrete version of the Eisehart's formulation of the Jonas trans- 
formation. One derives from eqs. (|2.13|) " H2.15|1 the following formulas: 

l + nC-^))(l + -)=l + {l + c^A + (1 + (2.16) 
V a, / / V a, / ai o,- 



1 + nbjj {l + bij=l + {l + c±)bj + (1 + cpbi. (2.17) 

The symmetry of their right-hand sides implies that there exist functions 
<j), (j)^ : Z"^ M (associated to points of Z"^) such that 

— - = 1 H , — — = l + 6i, l<i<m. (2.18) 

These functions are defined uniquely up to constant factors, which can be 
fixed by requiring 0(0) = (0) = 1. Moreover, eqs. (|2.16|) . (|2.17|) imply 
that the functions (j), (j)'^ satisfy the following equations: 

SiSjcf) = Cij6j4> + Cji6i4>, (2-19) 
SiSjcj)-^ = ct.6j<l)+ + c±5icl)+, (2.20) 

for all 1 < i 7^ j < m. Thus, like in the smooth discrete Jonas 

transformation yields additional scalar solutions (p and of the equations 
describing the nets / and respectively. The solution (j) is directly spec- 
ified by the initial data (J^). Introduce the functions g : Z"^ and 
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ip : Z'" — > M by the formulas ()1.19|) . so that also the classical represen- 
tation (|1.14|) remains valid. A direct computation based on eqs. (|2.12j) . 
(|2.13|1 - H2.15|) . and ()2.18|) shows that the following equations hold: 

Sig = aidif, (2.21) 
6itp = ai6i(t>, (2.22) 



where 



"i = ■ (2-23) 



Next, one checks that the quantities Oj satisfy the equations 

5iaj = Cij{Tjai - Tiaj). (2.24) 

The same argument as in the smooth case shows that the data (J^) yield 
the values of (/>"'", and thus the values of Oi, on the coordinate axes Sj. This 
uniquely specifies the solutions of the compatible linear system (|2.24|) . 
which, in turn, allows for a unique determination of the solutions g, of 
eqs. (fOTT) . (l2:22ll with the initial data ^(0) = /+(0) - /(O) and ip{0) = 1. 

Continuous limit. Observe that eqs. (|2.1j) . (|2.2|) are quite similar to 
eqs. (|1.2() characterizing smooth conjugate nets. We will demonstrate 

in Sect. 01 that the status of this similarity can be raised to that of a mathe- 
matical theorem about approximation of smooth conjugate nets by discrete 
ones. More precisely, we will show how to choose initial data for a discrete 
system (with M = m and a small mesh size e) so that it approximates a 
given m-dimensional smooth conjugate net as e ^ 0. 

Analogously, eqs. (|TT^ - (|TT^ are similar to eqs. (frTU)) - (P~n|) . Ac- 
cordingly, initial data of a discrete system with M = m + 1 can be choosen 
so that, keeping one direction discrete, one arrives in the limit at a given 
smooth conjugate net and its Jonas transform. 

For M = m + 2 and M = m + 3, keeping the last two, resp. three 
directions discrete, one proves the permutability properties of Jonas trans- 
formations formulated in Theorem II. 31 Thus, permutability of Jonas trans- 
formations, which is a non-trivial theorem of differential geometry, becomes 
an obvious consequence of the multidimensional consistency of discrete con- 
jugate nets, combined with the convergence result mentioned above. 



2.2 Discrete orthogonal nets 

Two-dimensional circular nets (M = 2) were introduced in |MPSl INuj 
as discrete analogs of the curvature lines parametrized surfaces. A dis- 
cretization of triply orthogonal coordinate systems {M = 3, N = 3) was 
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first proposed in |Bobj . The next crucial step was done in |CDSj . where 
discrete orthogonal nets were generalized to arbitrary dimensions by con- 
sidering them as a reduction of Q-nets. For further developments, see 
jP oMS , KoSdH [KVhKV. PoS2, RobHej . 

Definition 2.6 A map f : —>■ M.'^ is called an M -dimensional discrete 
0-net (discrete orthogonal net, or circular net) in R^, if any of its ele- 
mentary quadrilaterals is circular, i.e., if at any u € and for all pairs 
^ 1^ i ^ j 1^ M the four points f, Tif , Tjf and TiTjf are concircular. 

It is important to observe the Mobius invariance of Definition 12.61 To un- 
derstand restrictions imposed on Q-nets by the circularity condition, we 
consider again various values of M. 

M=2: discrete surfaces parametrized along curvature lines. 

Suppose two coordinate lines / and / [■32 on a discrete 0-surface / 
are given. An elementary inductive step for extending the 0-surface to the 
quadrant Z^, consists of choosing /12 on the circle through /, fi and /2. In 
doing so, one has the freedom of choosing one real parameter at each such 
step, for instance, the cross-ratio of four points qi2 = q{f, /i, fi2, f2), which 
is naturally attached to the elementary square (n, n-|-ei,u-|-ei + e2,u + 62). 
Thus, to define a discrete 0-surface /, one needs to prescribe the coordinate 
lines /ts^ and /fsj and one real- valued function qi2 defined on elementary 
squares. 

However, it turns out to be technically more convenient to use other 
functions on elementary squares of characterizing the form of circular 
quadrilaterals, namely, discrete rotation coefficients. Introduce discrete met- 
ric coefficients hi = \6if\ and unit vectors Vi = h~^6if, so that 

Tif = f + hvi. (2.25) 

Then the rotation coefficients (3i2 and /32i are defined by the formula 

'Ti'^j = ^Ji'i'^j + l^jiVi), (2.26) 

which holds due to planarity of the elementary quadrilateral, with Vji be- 
ing appropriate normalization coefficients. From (|2.25|) . (|2.26|1 there easily 
follows a formula for metric coefficients hj-. 

Tihj = iy-^{hj + hiPij). (2.27) 

It is natural to assume that the variables Vi, hi are attached to the edges of 
Z^ parallel to the coordinate axis Sj, while the rotation coefficients /3i2,/?2i 
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are attached to elementary squares of Z^. Elementary considerations of 
the geometry of a circular quadrilateral on Fig. 12.1^1 show that the discrete 
rotation coefficients are related by 

Pij+Pji + 2{vi,vj) = 0, (2.28) 

while the normalization coefficients satisfy 

4 = 4 = 1- A,/?,,. (2.29) 

For an embedded quadrilateral there holds Uij = Vji > 0, while for a non- 
embedded one there holds Uij = —Vji- Thus, formula (|2.28j) expresses the 




Figure 2.3: An elementary quadrilateral of a discrete orthogonal net 

circularity constraint, and reflects the fact that the number of independent 
plaquette functions necessary to define a discrete 0-surface is equal to 1. 
For instance, it is enough to prescribe /3i2, or /32i, or /3i2 — /32ii or any other 
function of /3i2, (321, independent on /3i2 + /?2i- 

Equations 1)2. 26() . ()2.27|) give the evolution of the edge variables Vj, hj, 
known the plaquette variables P12, ^21- Together with eq. (|2.25j) . this allows 
us to reconstruct a discrete 0-surface /, provided the following initial data 
are given: two coordinate lines /ts^, /Tsj (o^' equivalently, the point /(O) 
and the functions Vi, hi on the edges of the coordinate axes Sj, i = 1,2), 
and the function P21 (say) defined on all elementary squares of Z-^. 

As in the case of Q-surfaces, the combinatorics of discrete 0-surfaces may 
be more complicated than that of Z^, because Definition 12 . 61 can be literally 
extended to an arbitrary quad-graph. This possibility is very important 
from the geometric point of view, since vertices of valence different from 4 
serve as a model for umbilic points of a smooth surface parametrized along 
curvature lines. 
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Figure 2.4: An elementary hexahedron of a discrete orthogonal net 



M=3: basic 3D system. Suppose that three coordinate surfaces of a 
three-dimensional discrete 0-net / are given, that is, /["2i2; / 1^223 s-^id /tsis- 
An inductive extension step consists of determining /123, provided /, and 
fij are known for all 1 < i < j < 3 and satisfy the circular condition on 
the three corresponding squares. Such a step is possible due to the Miguel 
theorem from the elementary geometry: 

Theorem 2.7 Given three circles Cij and seven points f,fi,fij in (1 < 
i < j < 3), such that each quadruple {f,fi,fj,fij) lies on Cij, define three 
new circles TiCjk as those passing through the triples {fi, fij, fik), respec- 
tively. Then these new circles intersect at one point: 

/123 = ^1^23 n T2C31 n T3C12 . 

The claim of the Miquel theorem is equivalent to the following one: if seven 
points /, fi, fij of an elementary hexahedron of a Q-net / : — > lie on 
a two-sphere 5^, then so does the eighth point /i23. To see the equivalence, 
note that the four points /, fi determine the sphere S"^ uniquely. Then 
fij G is equivalent to fij € Cij = Hij fl S"^. Now, the intersection point 
of the three circles TiCjk = Tjlljfc n S'^ belongs to and coincides with the 
unique intersection point of the three planes 

/123 = ^11123 n T2n3i n r3ni2. 

This construction is illustrated on Fig. 12.41 Clearly, this is a particular 
issue of Theorem 12.41 Thus, finding the eighth point of a three-dimensional 
discrete 0-net out of the seven known ones is a 3D system in the sense of 
Fig.O 
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The Miquel theorem guarantees that circularity constraint propagates 
in the construction of a Q-net from its three coordinate surfaces. This was 
first observed by Dohwa and Santini |CDSj : 

Theorem 2.8 If three coordinate surfaces /t2i2; f\'B2a' f^'Bia '^f Q-net 
f : Zi^ ^ M.^ are discrete 0-surfaces, then the whole of f is a discrete 0-net. 

This is a discrete analog of the statement which holds also in the smooth 
context, and is a sort of inversion of the classical Dupin theorem. 

On the level of formulas, for discrete 0-nets the system encoded on 
Fig. EU consists of eqs. (jS^SJ, ^TM . itOTll for all 1 < i 7^ j < 3, with 
the normalization coefficients given in (|2.29|) . which have to be augmented 
by evolution equations for the rotation coefficients. These latter ones repre- 
sent the solvability of the system of linear equations (|2.26|) and are derived 
from the requirement Ti{TjVk) = Tj{TiVk). Like in the case of Q-nets, these 
equations are decoupled from the rest ones, and generate a well defined map 
{Pkj} ^ {nl3kj}- 

nPkj = {i^kinjY^ {Pkj + (3ki(3ij). (2.30) 

However, this evolution of the discrete rotation coefficients (3kj can be cou- 
pled to the evolution ()2.26|) of the unit vectors Vi by means of the circularity 
constraint 1)2. 28() . It follows from the Miquel theorem, and can be easily 
checked analytically, that this constraint propagates in the coupled evolu- 
tion. 

M > 4- consistency. The 4D consistency of discrete 0-nets is a 
consequence of the analogous property of Q-nets, since the 0-constraint is 
compatible with the Q-property. 

On the level of formulas we have for M > 4 the system (|2.25|) . (|2.26|) . 
(|2.27l) . (|2.3U|) . augmented by the circularity constraint 1)2. 28() . where now 
all indices i,j,k vary between 1 and M. For any triple of pairwise differ- 
ent indices {i,j,k), equations involving these indices solely, form a closed 
subset. The M-dimensional consistency of this system means that all three- 
dimensional building blocks can be imposed without contradictions. Initial 
data which allow for a unique solution of this system consist of: 

(of) values of / on the coordinate axes (or, what is equivalent, value 
/(O) and values of Vi and hi on all edges of the coordinate axes Sj) 
for 1 < i < M; 

(of) values of M{M — l)/2 functions (say) on all elementary squares 
of the coordinate planes "Bij for 1 < i < j < M. 
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At this point, it will be useful to observe the similarity of eqs. ()2.25|) . 
^riE^ . ^rm . ^nWHi and the constraint fTM with the system (rni|l - (fnTll 

and the constraint ()1.24|1 governing smooth orthogonal nets. Like in the 
case of conjugate nets, we will demonstrate that this analogy can be given a 
qualitative content, so that smooth 0-nets can be approximated by discrete 
ones. However, there is a substantial obstruction in performing this. We 
think of smooth rotation coefficients as being approximated by discrete ones. 
But since the discrete rotation coefficients (3kj only have i ^ k,j as evolu- 
tion directions (that is, they are plaquette variables attached to elementary 
squares parallel to 23^^), there is seemingly no chance to get an approxima- 
tion of such smooth quantities as diPij involved in the smooth orthogonality 
constraint ()1.24() . In order to be able to achieve such an approximation, we 
need some discrete analogs of the smooth rotation coefficients which would 
live on edges. Such analogs will be introduced with the help of the notion of 
a frame of a discrete 0-net. Technical means for defining and constructing 
frames are given by the apparatus of Mobius differential geometry. 

Transformations of discrete O-nets. A natural generalization of 
Definition 11.51 would be the following one. 

Definition 2.9 A pair of m- dimensional discrete O-nets f,f^ : 
is called a Ribaucour pair, if four points f, Tif, f^ and Tif'^ are concircular 
at any u G and for any 1 < i < m. The net is called a Ribaucour 
transform of the net f . 

But this simply means that, if we set F{u,0) = f{u), F{u, 1) = f^{u), then 
F : X {0, 1} — > is an M-dimensional discrete 0-net, where M = m+l. 
So, in the discrete case there is no difference between orthogonal nets and 
their Ribaucour transformations. To specify a Ribaucour transform /+ of a 
given m-dimensional discrete 0-net /, one clearly needs the following data: 

(Rf) value of /+(0); 

(R^) values of ^Mi (say) on "vertical" elementary squares attached to all 
edges of the coordinate axes for 1 < i < m. 

Mobius-geometric description of discrete O-nets. Putting dis- 
crete O-nets in the Mobius-geometric model Q-^ of the Euclidean space M^, 
we observe first of all that Theorem 1 1 . 71 admits an almost literal generaliza- 
tion to the discrete case. 

Theorem 2.10 A Q-net f ^ is a discrete 0-net, if and only 

if l/p satisfies the same equation \2.1\) as f does, in other words, if the 
corresponding f : Z*^ — > Qq' is a Q-net in R^"*"^'-^. 
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(The first claim here is due to |KoSchj .^ Thus, in the Mobius-geometric 
picture discrete 0-nets are a quadratic reduction (to Qq) of Q-nets. 

Let Vi be the (Lorentz) unit vectors parallel to 5if. For the discrete 
metric coefficients hi = \Sif\ there holds also hi = \dif\- One readily verifies 
that from (/, /) = there follows (rj/ + f,Vi) = and hi = —2{f,Vi). 
Therefore, Vi are interpreted in the Mobius-geometric picture as reflections 
taking / to nf: 

Tif = f + hiVi = Vifvi. (2.31) 

Due to Theorem I2.1(J1 vectors Vi satisfy the same linear relations (|2.26jl as 
Vi, with the same discrete rotations coefficients l3ji. 

Theorem 2.11 (Spinor frame of a discrete O-net) For a discrete O- 
net f : Z^^ ( and the corresponding Q-net f : Qq ), there exists 

a function ip : Z*^ — !Koo ( called a frame of f) such that 

f = V'-'eoV, (2.32) 

satisfying the system of difference equations: 

nij; = -Sitl^Vi , l<i<M. (2.33) 

Proof. This theorem is due to |BobHej . Circularity of an elementary quadri- 
lateral of / implies that the angle between 5if and Ti5jf and the angle 
between 6jf and Tj6if sum up to tt: 

{TjVi,Vj) + {vi,TiVj) = 0. (2.34) 

The same relation holds for the vectors Vi. Moreover, since for these vectors 
there holds also eq. ()2.26|) with z/jj = Vji, there follows that eq. 1)2. 34() for 
Vi is equivalent to 

VjiTjVi) + Viinvj) = 0. (2.35) 

But this is exactly the compatibility condition of the system ()2.33|) of linear 
difference equations for tp. Indeed, the solvability condition TjTiip = TiTjip is 
written as 

eiejipVj{TjVi) = ejeiipVi{TiVj), 

which is equivalent to (|2.35|) . Thus, a solution to (|2.33j) exists and is uniquely 
defined by the choice of 'ip{0) G 3<oo- Choose 1^(0) so that eq. (|2.32|) holds 
at M = 0. Then eqs. (|2.33|) . (|2.31|) imply that eq. (|2.32|) holds everywhere 
on Z^. □ 



54 



CHAPTER 2. DISCRETE DIFFERENTIAL GEOMETRY 



Now introduce vectors Vi = TpViip~^, so that the frame equations ()2.33|) 
take the form {Tiip)ip~^ = —eiVi. Expanding these vectors with respect to 
the basis vectors e^, we have a formula analogous to ()1.34|1 : 

Vi = IpVii)''^ = (^i^i Pki^k + hiGoo- (2.36) 

The fact that the eoo-component here is equal to hi, is easily demonstrated. 
Indeed, from eq. (|2.33|) there follows that Tj/ — / = hiVi = hi{Tiip)~^eiip. 
Now eq. ()2.32() allows us to rewrite this equivalently as [eo, {Ti^)Tp~^] = hi^i, 
which proves the claim above. Observe also the normalization condition 

Coefficients p^i are edge variables analogous to smooth rotation coeffi- 
cients. Indeed, vectors Vi are defined on edges of Z*^ parallel to the coor- 
dinate axis Sj, but they do not immediately reflect the local geometry near 
these edges. Rather, they are obtained by integration of the frame equations 
1)2. 33() . and thus are of a non-local nature. Thus, in the discrete case we have 
two different analogs of the rotation coefficients: local plaquette variables 

for 1 < i 7^ j < M, defined on elementary squares of Z*^ parallel to Sj-,-, 
and non-local edge variables pki for 1 < i < M, l<k<N,k^i, defined 
on edges of Z,'^'^ parallel to "Bj. 

Evolution equations for Vi are obtained from (|2.26p and the frame equa- 
tions (lOHl) : 

In the derivation one uses the identity Vi{vj + (3jiVi)vi = Vj + PijVi , which 
follows easily from ()2.28|) . The resulting evolution equations for the edge 
variables pkj read: 

nPkj = J^Jl^iPkj + PkiAj), (2-38) 
npij = vj^{-pij + 2(7 i(3ij). (2.39) 

Here 1 < « 7^ j < M, 1 < A; < A^, and k ^ The circularity constraint 
(|2.28|) can be now written as 

Pij + Pji = CFiPij + ajpji - ^ X] PkiPkj^ (2-40) 
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and gives a relation between local plaquette variables fiij and non-local edge 
variables pkj. The system consisting of (|2.|-{8|) . (|2.l'{9|) and (|2.4()j) can be 
regarded as the discrete Lame system. 

Continuous limit. In Sect. 13.41 we will show how to choose these data 
on the lattice with mesh size e in m directions and mesh size 1 in the M-th 
direction, in order to achieve a simultaneous approximation of a smooth or- 
thogonal net and its Ribaucour transformation. Performing a smooth limit 
so that two or three of the directions remain discrete, one can derive the 
permutability properties of Ribaucour transformations formulated in Theo- 
rem 11.61 This way, permutability of Ribaucour transformations becomes a 
simple consequence of properties of discrete orthogonal nets, combined with 
the convergence result mentioned above. 



2.3 Discrete Moutard nets 

Discrete Moutard nets were introduced in |NiSchj . 

Definition 2.12 A map f -.l? ^ is called a two-dimensional discrete 
M-net, if it satisfies the discrete Moutard equation 

riT2f + f = au{nf + T2f) (2.41) 

with some ai2 : — > M. 

Initial data that can be used to determine a discrete M-net are: 

(M^) values of / on the coordinate axes Si, 'B2; 

(M^) function ai2 defined on elementary squares of Z^. 

Like in the smooth case, discrete M-nets are not Q-nets, and there are no M- 
dimensional discrete M-nets with M > 3. There is a construction, analogous 
to the smooth case, relating discrete M-nets to a special class of Q-nets. Let 
: — > M be any solution of the same discrete Moutard equation ()2.41|) 
(for instance, any component of the vector /), then y = u^^f : —>■ is 
a Q-net: 

y\2-y = — -, — ■ — r m -y) + — — ■ — r (2/2 - y)- 

However, there exists a relation to discrete conjugate nets of a quite different 
flavor, which is of a purely discrete nature and has no smooth analogs. To 
describe it, perform in (|2.43|) the change of variables 

f{u) ^ (-1)«2 f{u), u = (m, ns) G I?. (2.42) 
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Then the new function / satisfies the following modified form of the discrete 
Moutard equation: 

Tir2f - f = ai2{T2f - Tif ). 

This equation admits a multidimensional generalization: 

Definition 2.13 A map f : Z^'^ — > is called an AI -dimensional T-net 
(trapezoidal net), if for any u G T,^^ and for any pair of indices i ^ j there 
holds the discrete Moutard equation 

T^TJf - f = aij{Tjf -nf), (2.43) 

with some aij : Z*^ — M, in other words, if all the elementary quadrilaterals 
{f,Tif,TiTjf,Tjf) are planar and have parallel diagonals. 

Of course, coefficients Oij have to be skew-symmetric, aij = —aji. As usual, 
we will consider these functions as attached to elementary squares of Z^^. 

T-nets, unlike discrete M-nets, form a subclass of Q-nets. The condition 
of parallel diagonals is expressed as Cij + Cjj + 2 = for the coefficients Cij 
of a Q-net, the skew-symmetric coefficients of the T-net being Oj,- = Cjj + 1. 



M=2: T-surfaces. To define a two-dimensional T-net f : 1? — 
one can prescribe two coordinate curves, /f^i and f 1^,-2 ■, and a real- valued 
function ai2 on elementary squares of W?. 

M=3: basic 3D system. We show that three-dimensional T-nets are 
described by a well-defined three-dimensional system. An inductive con- 
struction step of the net / is as follows. Suppose that /, fi and fij are 
given for all 1 < i 7^ j < 3, satisfying eq. (|2.43j) . Three equations (|2.43j) for 
the facets of an elementary cube on Fig. 12.11 adiacent to /123, lead to con- 
sistent results for /123 for arbitrary initial data, if and only if the following 
conditions are satisfied: 

inajk)aij + {Tjaki){ajk + Oij) = —1, 

where (i, j, k) is an arbitrary permutation of (1, 2, 3). These conditions con- 
stitute a system of 6 (linear) equations for 3 unknown variables Tittjk in terms 
of the known ones Ojk- It turns out that this system is not overdetermined 
but admits a unique solution: 

^ = , (2.44) 

With TiOjk so defined, eqs. (|2.43|) are fulfilled on all three quadrilaterals 
adjacent to /123. 
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Eqs. ()2.44l) represent a well-defined birational map {ajk} i— > {Tiajk}, 
which can be considered as the fundamental 3D system related to T-nets. 
It is sometimes called the "star-triangle map". Moreover, this system is a 
reduction of the system describing Q-nets: 

Theorem 2.14 If three coordinate surfaces f\'Bi2> /ts23' f\'Bi3 of a Q-net 
f -.l? ^ are T-nets, then f is a T-net. 

Proof. Let three quadrilaterals {f,fi,fij,fj) be planar and have parallel 
diagonals. The planarity of the three quadrilaterals (/j, fij, fijk, fik) defines 
the point /123 as the intersection point of three planes planes Tjlljfc. Then 
these three quadrilaterals automatically have parallel diagonals. Indeed, by 
the above argument, there exists a point /123 with this property, and it has 
to coincide with the one defined by the planarity condition. □ 

M > 4- consistency. The 4D consistency of T-nets is a consequence 
of the analogous property of Q-nets, since T-constraint is compatible with 
the Q-property. On the level of formulas we have for T-nets with M > 4 
the system (|2.43|) . (|'2.44() . All indices i,j, k vary now between 1 and M, and 
for any triple of pairwise different indices {i,j, k), equations involving these 
indices solely, form a closed subset. Initial data which allow for a unique 
solution of this system consist of values of / on the coordinate axes "Bj for 
1 < i < M, and values of M{M — l)/2 functions Oij on all elementary 
squares of the coordinate planes Sjj for 1 < i < j < M. 

Transformations of discrete M- and T-nets. Because of the multi- 
dimensional consistency, transformations of T-nets do not differ from the 
nets themselves, and are described by the discrete Moutard equations. We 
would like, however, to consider transformations of discrete M-nets, which 
have geometric applications. For this, we consider a three-dimensional T-net 
F : 7? X {0, 1}, perform the change of variables 

F{u) ^ {-ir^F{u), U = (U1,U2,U3), 

and then set / = F(-,0) , /+ = F{-, 1). Eqs. (lO!! for = (1,2) turn 

into 1)2. 41() . so that /, /+ are discrete M-nets, while for = (1,3) and 
(2, 3) eqs. turn into 

ri/+ - / = 6i(/+ - n/), r2/+ + / = 62(/+ + r^f), (2.45) 

where hi = 013 , 62 = 032. The quantities ftj, defined on the "vertical" 
plaquettes of T? x {0, 1}, parallel to ^43, can be also associated to the edges 
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of Z'^ parallel to the coordinate axes Sj. Eqs. ()2.44|) with one of the indices 
equal to 3 express compatibility of eqs. (|2.45|) with eq. (|2.43|) : 



T2&1 _ _ "l2 _ 1 

bi 62 (112 {h + 62)012 - bib2 



(2.46) 



Eqs. (|2.45|) . 1)2. 46(1 define a discrete Moutard transformation of /. To specify 
a Moutard transform of a given discrete M-net /, one can prescribe the 
following data: 

(MTf) value of /+(0); 

(MT^) values of bi on "vertical" elementary squares attached to all edges 
of the coordinate axes for i = 1,2. 

The permutability properties of Moutard transformations are governed by 
the discrete Moutard equations (|2.43j) with i,j > 3. 

Due to the first equation in 1)2. 46() . there exists a function : ^ M 
(associated to the points of Z^) such that 



Tit 



1,2. (2.47) 



The last equation in ()2.46|) implies that this function is a scalar solution of 
the discrete Moutard equation (|2.41|) . This solution is specified by its values 
on the coordinate axes 23j, which are immediately obtained, via eq. (|2.47l) . 
from the data (MT^). Recall that /+ is a solution of the discrete Moutard 
equation 1)2. 41() with the transformed potential a^2- The second equation in 



(|2.46|) gives a representation of in terms of 6: 

(rig)(r2g) _ +_1 

This is a discrete analog of the classical formulation of the Moutard trans- 
formation. 



2.4 Discrete asymptotic nets 

The following definition is due to Sauer |Saj for M = 2, when it describes a 
discrete analog of surfaces parametrized along asymptotic lines ( A-surfaces) . 
For M > 3, see IDoJNSTl lNie| . 
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Definition 2.15 A map f : M.'^ is called an M -dimensional discrete 

A- net (discrete asymptotic net) in M^, if for any u G all the points 
f^uzizCi) lie in some plane through f{u). 

Note that this definition belongs to the projective geometry. For M = 3, the 
geometry of an elementary hexahedron of a discrete A-net is exactly that 
of a Mohius pair of tetrahedra, i.e., a pair of tetrahedra which are inscribed 
in each other. In other words, each vertex of each tetrahedron lies in the 
plane of the corresponding facet of the other one (eight conditions). In our 
case the two tetrahedra are those with the vertices (/, /12, /23, /31) and with 
the vertices (/i, /2, /s, /i23)- Such pairs were introduced by Mobius |Moj . 
who demonstrated that eight conditions mentioned above are not indepen- 
dent: any one of them follows from the remaining seven. Mobius pairs are 
remarkable and well-studied objects of the projective geometry, cf. |Bobj . 

For a non-degenerate discrete A-net, all quadrilaterals (/, r,/, TiTjf, Tjf) 
are non-planar. It would be in principle possible to consider discrete A- nets 
in M.^ with > 3, however it would not lead to an essential generaliza- 
tion. Indeed, for any fixed u E and for any pair of indices i ^ j from 
{!,... ,M}, consider the three-dimensional affine subspace of through 
/ = f{u) which contains Tj/, Tjf and TiTjf. A simple induction shows that 
the whole net / lies in this subspace. 

For a discrete A-net / : Tj^'^ — > M.^ we have a field of tangent planes 
y : — > Gr2(3), and therefore a well-defined normal direction at every 
point of Z^^. A remarkable way to fix a certain normal field is given by the 
discrete Lelieuvre representation |KoPj which states: 

Theorem 2.16 For a non- degenerate discrete A-net f , there exists a nor- 
mal field n : Z^^ such that 

5if = Tinxn, i = l,...,M, (2.49) 

called a Lelieuvre normal field. It is uniquely defined by a value at one point 
uq G Z*^. All other Lelieuvre normal fields are obtained by n{u) ^ an{u) 
for \u\ = til -|- . . . -|- um even, and n{u) 1-^ a~^n{u) for \u\ odd, with some 
a G M (black-white rescaling). 

It follows from eq. ()2.49|) immediately that {TiTjU — n) x (Tin — Tju) = 0, 
that is, the Lelieuvre normal field satisfies the discrete Moutard equations 

TiTjU — n = aij{Tjn — Tin) (2.50) 

with some Oij : Z''^ — > M. Conversely, given a T-net n : Z*^ M^, formula 
(imH) produces a discrete A-net / : Z^ ^ M^. 
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Theorem 2.17 Discrete A-nets in are in a one-to-one correspondence, 
via the discrete Lelieuvre representation \2.4y\) , with T-nets in M^. 

In particular, the initial data which determine a discrete A-net are analogous 
to the data (Mf^g) for the Lelieuvre normal field: 

(A^) values of n on the coordinate axes 23 j for 1 < i < M; 

{A^) values of M(M - l)/2 functions Qij on all elementary squares of the 
coordinate planes for 1 < i < j < M . 

Transformations of discrete A-nets. A natural generalization of 
Definition 11.51 would be the following one. 

Definition 2.18 A pair of discrete A-nets f,f^ ■ — > is called a 
Weingarten pair, if, for any u E Z™, the line segment [f{u),f^{u)] lies in 
both tangent planes to f and at the points f{u) and f~^{u), respectively. 
The net f~^ is called a Weingarten transform of the net f . 

But, as usual, this definition means simply that the net F : x {0, 1} ^ 
with F(n, 0) = f{u) and F{u, 1) = f^{u) is an M-dimensional discrete A- 
net, where M = m + 1. So, once again, transformations of discrete A-nets do 
not differ from the nets themselves. The Lelieuvre representation of the net 
F can be written now as a relation between the Lelieuvre representations of 
the nets f,f^: 

f+-f = n+xn. (2.51) 

Definition 12.181 makes sense for an arbitrary M; however, we will be 
mainly interested in the case M = 2, which is an immediate discretization 
of the smooth A-surfaces. The change of variables 1)2.42(1 for the Lelieuvre 
normal field, 

n(n) (-1)''2 n(n), u = {ui,U2) e , (2.52) 

leads to a replacement of the general Lelieuvre formulas ()2.49|) by 

5if = Tin X n = 5in x n, ^2/ = n x T2n = n x 52n. (2.53) 

For a Moutard transformation of the Lelieuvre normal field there hold for- 
mulas ^ZB- 

Tin~^ — n = bi{n~^ — Tin), T2n~^ + n = b2{n~^ + T2n). (2.54) 

Thus, a Weingarten transform f^ of a given discrete A-net / is determined 
by a Moutard transform n"*" of the Lelieuvre normal field n, and in order to 
specify the latter, one can prescribe the following data: 
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(Wf) value of n+(0); 

(W^) values of bi on "vertical" elementary squares attached to all edges of 
the coordinate axes Hi for i = 1,2. 

Permutability of discrete Weingarten transformations is a direct consequence 
of the 4D consistency of T-nets. Permutability of smooth Weingarten trans- 
formations (Theorem I1.16() will follow, if one combines the discrete per- 
mutability with the convergence results of Sect. 13.61 

2.5 Discrete K-nets 

In discretizing the K-surfaces and their transformations, we take as a starting 
point the characterization of Theorem II. 181 

Definition 2.19 A discrete A-net f : Z*^ — s- is called an M-dimensional 
discrete K-net, if for any elementary quadrilateral {f,Tif,TiTjf,Tjf) there 
holds: 

\nTjf - Tjf\ = \Tif - f\ and \TiTjf - Tif \ = \Tjf - /I, 

in other words, if for any i = 1, . . . , M , the function at = \6if \ (defined on 
the edges parallel to the coordinate axes Hi) depends on Ui only. 

This notion is due to Sauer |Saj in the case M = 2 and to Wunderhch |WI 
in the case M = 3. A study of discrete K-surfaces within the framework of 
the theory of integrable systems was performed in Bo bPlj . 

A characterization of the Lelieuvre normal field of a discrete K-net is 
analogous to the smooth case. 

Theorem 2.20 The Lelieuvre normal field n : Z*^ ^ M'^ of a discrete K- 
net f : Z^^ —>■ M.^ takes values, possibly upon a black-white rescaling, in some 
sphere C . In case when the radius of this sphere is equal to 1 ( which 
we associate to the discrete Gaussian curvature K = — Ij, the Lelieuvre 
normal field coincides with the Gauss map. Conversely, any T-net n in the 
unit sphere is the Gauss map and the Lelieuvre normal field of a discrete 
K-net f with K = —1. The functions f3i = \5in\ depend on Ui only, and are 
related to the functions Oi = \Sif \ by Oi = (3i{l — [3f /'i)^/'^ . 

Thus, discrete K-nets are in a one-to-one correspondence with T-nets in 
S^. We proceed to the study of the latter object. 

M=2: basic 2D system. If two coordinate curves of the Gauss map 
n of a discrete K-surface / are given, that is, n and n f^j) then one 
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can extend n to the whole of Z^. An inductive step in extending n to 
consists of computing TiT2n = n + ai2(r2n — Tin), where the coefficient au 
(attached to every elementary square of the discrete surface) is determined 
by the condition that TiT2n G S^, so that 



ai2 



{n,Tin - T2n) 
(nn, T2n) 



1 



This elementary construction step, i.e., finding the fourth vertex of an ele- 
mentary square out of the known three vertices, is symbolically represented 
on Fig. E31 



n2 



"-12 



n 



ni 



Figure 2.5: 2D system on an elementary quadrilateral 



M > 3: consistency. Turning to the case M > 3, we see that one can 
prescribe all coordinate lines of a T-net n in S^, i.e., n|"s. for all 1 < i < M. 
Indeed, these data are independent, and one can, by induction, construct 
the whole net from them. The inductive step is essentially two-dimensional 
and consists of determining TiTjU, provided n, r^n and tju are known. In 
order for this inductive process to work without contradictions, equations 

TiTjU -n = aij[Tjn - Tin), aij = ^ (2.55) 

1 - {nn, Tju) 

must have a very special property. To see this, consider in detail the case 
of M = 3; higher dimensions do not add anything new. From n and rii 
one determines all riij uniquely. After that, one has, in principle, three 
different ways to determine ni23, from three squares adjacent to this point; 
see Fig. 12.61 These three values for ni23 have to coincide, independently of 
initial conditions. 

Definition 2.21 A 2D system is called 3D consistent, if it can be imposed 
on all two-dimensional faces of an elementary cube ofl?. 
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Figure 2.6: 3D consistency of 2D systems 



Theorem 2.22 The 2D system \2. 55\) governing T-nets in §^ is 3D con- 
sistent. 

Proof. This can be checked by a tiresome computation, which can be 
however avoided by the following conceptual argument. Note that T-nets in 
are a result of imposing of two admissible reductions on Q-nets, namely 
the T-reduction and the restriction to a quadric S^. This reduces the effective 
dimension of the system by 1 (allows to determine the fourth vertex of 
an elementary quadrilateral from the three known ones), and transfers the 
original 3D equation into the 3D consistency of the reduced 2D equation. 
Indeed, after finding ni2, n23 and ni3, one can construct ni23 according to 
the Q-condition (as intersection of three planes). Then both the T-condition 
and the S^-condition are fulfilled for all three quadrilaterals adjacent to 
rai23. Therefore, these quadrilaterals satisfy our 2D system. Note that this 
argument holds also for T-nets in an arbitrary quadric (not necessarily in 



Thus, we have for M > 3 a system consisting of interrelated two- 
dimensional building blocks (|2.55|) . for all pairs of indices i,j between 1 
and M. The M-dimensional consistency of this system means that all two- 
dimensional building blocks can be imposed without contradictions. A set 
of initial data which determines a solution of the system (|2.55|) consists of 

(K^) values of n on the coordinate axes for 1 <i < M. 



§2). □ 



Transformations of discrete K-nets. A natural generalization of 
Definition II. 191 would be the following one. 
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Definition 2.23 A Weingarten pair of discrete K-nets f,f^ : Z"* is 
called a Backlund pair, if the distance |/^ — /| is constant, i.e., does not 
depend on u ^ Z™. 

But, comparing tiiis with Definition 12.191 we see that the net F : 77^ x 
{0, 1} with F(n,0) = f{u) and F{u, 1) = f+{u) is an M-dimensional 

discrete K-net, where M = m+1. Once more, transformations of discrete K- 
nets do not differ from the nets themselves. To specify a Backlund transform 
/"*" of a given m-dimensional discrete K-net /, or a Moutard transform n+ 
of the Gauss map n, one can prescribe the following data: 

(B^) value of n+(0). 

We will be mainly interested in transformations of discrete K-surfaces. After 
the usual change of variables H2.52|) we will have for the Gauss map the 
equation 

/ , ^ {n,Tin + T2n) 

TiT2n + n = ai2[Tin + T2n), au = ——, r, (2.56) 

1 + {Tin,T2n} 

while for its transformation we will have the equations 
Tin~^ — n = 6i(n^ — Tin), 

T2n~^ + n = b2{n~^ + T2n), 

Of course, this has to be supplied with the Lelieuvre formulas (|2.53|) and 
1)2. 54() . Permutability of Backlund transformations for discrete K-surfaces is 
a direct consequence of the 3D consistency of T-nets in S^, and the smooth 
result (Theorem II. 2U() will be derived in Sect. 13.71 by the continuous limit. 



h 
b2 



{n. Tin — ri^) 
1 — (rin, n+) ' 
{n,T2n + n+) 
1 + (r2n,n+) ' 



(2.57) 
(2.58) 



2.6 Discrete isothermic nets 

Discrete isothermic surfaces were introduced in |RohP2j . See also |HeHP[ 
IHell lHe2l ISchlj for further developments. 

Definition 2.24 A discrete 0-surface / : — > is called a discrete 
I-surface (discrete isothermic surface), if the cross-ratios of its elementary 
quadrilaterals satisfy 

q{f,fi,fi2,f2) = -—, (2.59) 
a2 

where the functions ai depend on Ui only (i = 1,2), and the usual notations 
fi = Tif, fi2 = TiT2f are used. 
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It is natural to consider the functions a-i as defined on edges of 7? parallel 
to the coordinate axes Sj, with the property that any two opposite edges of 
any elementary square carry the same value of the corresponding Oj. Such 
functions will be called a labelling of the edges of Z^. Actually, it will be 
more convenient to associate —02 to the edges parallel to 232- To determine 
a labelling, one has to prescribe it on the coordinate axes 23 j. Thus, discrete 
I-surfaces are governed by a 2D equation (|2.59|) . and the following initial 
data determine such a surface completely: 

(I^) values of / in the coordinate axes Sj {i = 1, 2), i.e., two discrete curves 
fl'Bi with a common intersection point /(O); 

(I^) two functions Oj : ^ M on the edges of the coordinate axes for 
i = 1,2, which yield an edge labelling, ai = ai{ui). 

A discrete analog of the function s from eq. (|1.5H1 is defined on the vertices 

of Z2. 

Theorem 2.25 For a discrete isothermic surface f , there exists a function 
s : ^ M (positive if both ai > 0), such that 

\f^ - f\^ = aiSSi (i = l,2). (2.60) 

It is defined uniquely, up to a black-white rescaling (which is fixed by pre- 
scribing s at one point, say at u = 0). Conversely, existence of such a 
function s for a discrete 0-surface f implies that f is isothermic. 

The following property actually characterizes discrete isothermic surfaces. 

Theorem 2.26 (Dual discrete I-surface) Let / : — > be a discrete 
isothermic surface. Then the M.'^ -valued discrete one-form 6f* defined by 

mfr SSl |52/|2 SS2 

is closed. Its integration defines (up to a translation) a surface f* -.1? ^ 
, called dual to the surface f , or Christoffel transform of the surface f . 
The surface f* is discrete isothermic, with 

q{rJiJi2J2) = --- (2.62) 
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Proof. We rewrite eq. H'2.59() in several equivalent forms. Recall that in 
this equation we regard / as belonging to the Clifford algebra S^(M^), and 
that the cross-ratio is defined by eq. ()1.56|) . One checks straightforwardly 
that eq. (|2.59|) is equivalent to: 

(«i + a2)(/i2 - /)-' = ai(/i - + a2(/2 - /)~\ (2.63) 
or, interchanging the roles of / and /12, to: 

(ai + a2)(/i2 - = ai(/i2 - f2)'^ + a2(/i2 - /i)"^ (2.64) 
Taking into account that (/j - f)'^ = -{fi - f)/\fi - /P, we find: 

fi — f , f2 — f /12 - /2 , /12 — /l 
"117 7T2+«2r7 7T2=«1T7 712 + "277 TI2 , (2.65) 

1/1 -/r 1/2 -/r 1/12 -/2r l/i2-/ir 

or, due to 

/i — / , /2 — / /12 — /2 , /12 — /l cc\ 
1 = 1 , (2.66) 

ssi SS2 S2S12 S1S12 

This is equivalent to the first statement of the theorem. The second one 
follows readily from the Clifford algebra expressions fi — f* = ai(/i — f)~^ 
and /I-/* = -a2(/2-/)-^ Note that \f*-r\^ = a]\fi-f\-^ = ai{ss,)-\ 
so that s* = s^^ can be taken as the analog of the function s for the dual 
surface /*. □ 

Definition 2.27 (Discrete Darboux transformation) A Ribaucour 
transform : ^ of a given discrete isothermic surface f :1? ^ 
is called a Darboux transform, if the cross-ratios of its elementary quadri- 
laterals can he likewise factorized: 

q{f^JtJt2Jt) = --- (2.67) 

Oi2 

It can be demonstrated that, apart from the trivial case when /"*" is a Mobius 
transform of /, the Darboux transform is given by the following formulas: 

g(/, /i, /+) = ^, q{f, /2, f}, /+) = (2.68) 

where c G M is its parameter. Clearly, following data determine a Darboux 
transform /"*" of a given discrete isothermic surface / uniquely: 

(Df ) a point /+(0); 
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(D2 ) a real number c, the parameter of the transformation. 

Comparing formulas (|2.59|) . (|2.67j) . and (|2.68|) . we see that iterating a Dar- 
boux transformation is nothing but adding a third dimension to a two- 
dimensional discrete isothermic net, = /s. The parameter c plays the 
role of the function a^, attached to all edges parallel to the third lattice 
direction. Formula (|2.6fl|l for the third direction reads: 

I/+ _ /|2 = (2.69) 

(which literally coincides with the corresponding formula (jl.55|) for the 
smooth case) . All this is possible due to the following fundamental statement 
about multidimensional discrete isothermic nets: 

Theorem 2.28 The 2D cross-ratio equation 

9(/,/^,^-,/i) = - (2.70) 
aj 

is 3D consistent for any labelling Oi of the edges. 

This statement is due to |HeHPj : its generalization, where S^(M^) is replaced 
by an arbitrary associative algebra, was proven in |BobSu2j . and will be given 
in Theorem 14.41 

Remark. In view of the fundamental importance of the consistency 
property which holds if the right-hand side of the cross-ratio system is fac- 
torized as in (|2.7fl|) , one might wonder why we introduced the minus sign in 
eq. ()2.59|) for discrete I-surfaces, which propagated also into eq. ()2.68|) for 
the Darboux transformations. This choice is motivated by the convenience 
of passing to the continuous limit only. The fundamental factorization form 
(|2.7fl|) is restored by the change 02 — > —0:2- 

Mobius-geometric characterization of discrete isothermic sur- 
faces. Cross-ratio of four concircular points is a Mobius invariant quantity. 
Therefore, the notions of discrete isothermic surfaces and of their Darboux 
transformations are also Mobius invariant. To give their characterization 
within the Mobius-geometric formalism, we note first of all that, according 
to Theorem 12.101 eq. ()2.66|) holds also for the image / of the net / in the 
quadric Q'^: 

fl - f _|_ f2 - f _ fi2 - f2 _|_ /12 - fl 

SSl SS2 S2S12 S1S12 

This yields a discrete analog of Theorem ll.251 
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Theorem 2.29 The lift s = s"^/ : ^ L^+i of f to the light cone of 
gdiigji^g ifig discrete Moutard equation 

T1T2S + S = ai2{TlS + T2S), (2-71) 

with ai2 = {tiT2S~^ + s~^)/(ris~^ + T2S~^). 

Conversely, given a discrete M-net s : — > L^"*'"'^ in the light cone, let 
the functions s ^ M+ and f ^ he defined by 

s = s"^(/ + eo + l/peoo) 

(so that s~^ is the BQ-component, and s~^f is the -part of s in the basis 
ei, . . . ,6^,60,600). Then f is a discrete isothermic surface. 

Eq. H2.71() is a 2D equation in L^"^^ which determines T1T2S from s, tis, 
T2S. This yields: 

au = . (2.72) 

{TiS,T2S) 

iFvom this one readily sees that for a discrete M-net in L^+^ the quantities 
(s, Tjs) depend on Ui only. They are related to the labelling aj of Definition 
12.241 bv {s,Tis) = —ai/2. Thus, the labelling is already encoded in the lift s 
of a discrete I-surface /. 

For the consistent cross-ratio system ()2.7U|) one can introduce the scalar 
function s : Z*^ R (not necessarily positive) according to eq. ()2.6U|) with 
1 < i < M, which yields a multidimensional version of the lift s = s~^f : 

— > L^^"*^. It satisfies a discrete Moutard equation with minus signs: 

TiTjS — s = aij{TjS — Tis). (2.73) 

This is an instance of 3D consistent T-nets in a quadric (see the proof of The- 
orem r2.2'^ . and describes Darboux transformations of discrete I-surfaces and 
their permutability properties. To perform a smooth limit a usual change 
of signs (see Sect 12. 3() is required. This results in eq. (|2.71() for the surface 
and in 

tJ+ -s = bi{s+ -Tis), r2S+ + s = 62(5^ + T2s) (2.74) 



for its Darboux transformations. 
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2.7 Discrete surfaces made of spheres 

It is common to think of discrete surfaces as of maps / : ^ M^. However, 
it is not the only possibihty. One of the further ideas is to "blow up" the 
points /, making them to (small) hyperspheres in M^. We discuss here 
two realizations of this idea, which generalize the notions of discrete 0-nets 
and of discrete I-surfaces, maintaining important features of these notions, 
including transformations with permutability properties. 

S-orthogonal nets. We first generalize M-dimensional discrete 0-nets 
/ in M.^ . Blowing up the points to spheres is done in the following way: 

Definition 2.30 An S-circular net is a map 

S : |non-oriented hyperspheres in R^} (2.75) 

such that, for any elementary square of l?"^ , the corresponding four hyper- 
spheres S, TiS, TjS and TiTjS have a common orthogonal circle. 

Clearly, this reduces to Definition [THl if the radii of all hyperspheres become 
infinitely small. Observe also that the plane of the common orthogonal circle 
contains the Euclidean centers of all four hyperspheres. 

Recall (see Appendix to Chapter ^ that hyperspheres in can be 
represented as elements of P(M^,^^'^), where 

M^+i'i = {J G M^+i'i : {s,s)>0} (2.76) 

is the space-like part of M^+^'^. Thus, blowing up points to spheres cor- 
responds to the step away from L^"^'^ to M^^J^"'^'^. To establish a charac- 
terization of S-circular nets similar to that of Theorem I2.1U1 the following 
statement is needed. 

Theorem 2.31 Four points of M.^^^'^ are linearly dependent if and only if 
the corresponding four hyperspheres in M.^ 

(i) have a common orthogonal circle, or 

(a) intersect along an {N — 3) -sphere, or else 

(Hi) intersect at exactly one point. 

Case (Hi) can be regarded as a degenerate case of both (i) and (ii). 
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The distinction between these three cases depends on the signature of the 
quadratic form obtained as the restriction of the Minkowski scalar product 
to the {N — l)-dimensional (Lorentz-)orthogonal complement of the three- 
dimensional vector subspace of M^+^'^ spanned by our four points. This 
form is positive-definite in the case (i), it has the signature (A^ — 2,1) in 
the case (ii), and it is degenerate in the case (iii). It is understood that for 
N = 3 the 0-sphere in the case (ii) is a point pair. 

Since it is difficult to control the signature of this quadratic form, it 
makes sense to introduce a notion which is somewhat wider than that of a 
S-circular net. 

Definition 2.32 An S-orthogonal net is a map \2. 75\) such that, for any 
elementary square of Z*^, the points s, TiS, tjs and TiTjs, corresponding to 
four spheres S, TiS, tjS and TiTjS, are linearly dependent. 

For = 3, S-orthogonal nets with all intersections of type (ii) are nat- 
ural discrete analogs of sphere congruences parametrized along principal 
lines, because four infinitesimally neighboring spheres of such a congruence 
intersect this way, the pairs of intersection points comprising two enveloping 
surfaces of the congruence (see, e.g., |E2j ) . 

Theorem 2.33 Six circles (corresponding to the faces) of an elementary 
cube of a S-circular net lie on a 2-sphere. For N = 3 these spheres comprise 
an S-orthogonal net. 

A convenient choice of representatives s of hyperspheres S" in a fixed 
affine hyperplane of M^^J^^'^ is given by eq. Therefore: 

Theorem 2.34 A map \2. 75| ) is an S-orthogonal net, if and only if the 
corresponding map 

s : Z'" ^ M^+^'^ n {eo = 1}, 1 = c + eo + {\c\^ - r^)e^, 

is a discrete Q-net inMJ^~^^'^ , i.e., if the centers c : Z^'^ M.^ of the spheres 
S form a discrete Q-net in M^, and their radii r : Z*^ M_(. are such that 
the function |cp — satisfies the same equation \2.1\) as the centers c. 

A smooth version of this theorem for arbitrary sphere congruences in 
parametrized along principal lines can be found, e.g., in IE2 I . §99. 

The multidimensional consistency of discrete quadrilateral nets can be 
immediately transferred into the corresponding property of S-orthogonal 
nets, with the following reservation: given seven points s, TjS, TiTjs in 
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^out^'^ {?o = 1}) the Q-property (planarity condition) uniquely defines 
the eighth point T1T2T3S in R-^"^^'^ n {^0 = 1}, which, however, might get 
outside of M'^yjj^'^, and therefore might not represent a real hypersphere. 
Thus, the corresponding discrete 3D system is well-defined only on a proper 
open subset of the definition domain. As long as it is defined, it can be used 
to produce "S-Ribaucour" transformations of S-orthogonal nets, with usual 
permutability properties. 

S-isothermic surfaces. Recall that, according to Theorem 12.291 dis- 
crete I-surfaces are characterized as discrete M-surfaces (or T-surfaces) in 
l^N+i^ In this context, we blow up points to spheres by replacing L^+^ by 
L^+i (seeeqs. (fTHHl) . (jOTTl) !: 

Definition 2.35 A map 3:1?^ {oriented hyperspheres in M^} is called 
an S-isothermic surface, if the corresponding map s : — s- L^"^^ is a 
discrete M- surface (T- surf ace). 

S-isothermic surfaces, along with their dual surfaces were originally intro- 
duced in [BobP3 for the special case of touching spheres. The general class 
of Definition l2.3.^ together with Darboux transformations and dual surfaces, 
is due to |Hoj . 

S-isothermic surfaces are governed by the equation 1)2. 71(1 with 

{S,TIS + T2S) 
+ {tiS,T2S} 

and have the property that the quantities = {s, tis) depend on Ui only. 
If radii of all hyperspheres become uniformly small, r(n) ~ ks{u), k ^ 0, 
then in the limit we recover discrete isothermic surfaces, as characterized in 
Theorem 

Consistency of discrete T-nets in L'^^^ implies, in particular, Darboux 
transformations for S-isothermic surfaces, governed by eq. 1)2. 73(1 . A Dar- 
boux transform : Z"^ L^+^ of a given S-isothermic surface s is uniquely 
specified by a choice of one of its spheres s"'"(0). 

Geometric properties of S-isothermic surfaces are the following. First 
of all, S-isothermic surfaces form a subclass of S-orthogonal nets. Further, 
the quantities {s, Tis) depend on Ui only, and have the meaning of cosines of 
the intersection angles of the neighboring spheres (resp., of their so called 
inversive distances if they do not intersect). An important characterization 
is the following generalization of Theorem 12.261 



Theorem 2.36 (Dual S-isothermic surface) 
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Let S ^ {oriented hyperspheres in M^} he an S-isothermic surface. 
Denote the Euclidean centers and radii of S hy c :!? ^ and r : — > M, 
respectively. Then the M.^ -valued discrete one- form 5c* defined by 

dic = , O2C = (2.78) 

rri rr2 

is closed, so that its integration defines (up to a translation) a function 
c* : — M^. Define also r* : 7? ^ M. by r* = r^^. Then the hyperspheres 
S* with the centers c* and radii r* form an S-isothermic surface, called dual 
to S. 

Proof. Consider eq. ()2.71|) . in terms of 

s = r~^ (c + eo + (|cp - r^)eoo) • 

Its eo-part yields: ai2 = {r^2 + i'^^) / {"i^i^ + ^2^"*^)- This allows us to rewrite 
eq. (|T7T|) as 

(rfi + r^'){su + s) = {r^^ + r-^){h + h)- (2.79) 

A direct computation shows that the M^-part of this equation can be rewrit- 
ten as 

Cl — C Co — C Cto — Co Cio — Ci 

J: + ^ = _ti £ + (2.80) 

rri rr2 r2r\2 r\r\2 

which is equivalent to closeness of the form 5c* defined by H2.78() . In the 

same way, the e^o-part of eq. (|2.79|1 is equivalent to closeness of the discrete 

form bw defined by 

5i(|c|2-r2) 52(|c|2-r2) 

0\W = , O2W = . 

m rr2 
For similar reasons, the second claim of the theorem is equivalent to the 
closeness of the form 



<5i(|c*|2_(r*)2) _ 52(|c*r-(r 



Siz = — , 52Z 



I I ^ r*r2 
But one easily checks that 

ri r r2 r 

6lW = (Ci - C*, Cl + C) 1 , S2U) = (Co - C*, C2 + c) H 

r ri r r2 

6lZ = (Cl - C, cT + C*) - — + — , 62Z = (C2 - C, Co + c*) + — - — . 

ri r r2 r 

The sum of these two one-forms is closed: 

5i{w + z) = 2(ct, Cl) - 2(c%c), 62{w + z) = 2(c^, C2) - 2(c*, c), 

therefore they are closed simultaneously. □ 



Chapter 3 

Approximation 



3.1 Discrete hyperbolic systems 

To formulate the most general scheme covering all the situations encoun- 
tered so far, we have to put our hyperbolic systems into the first order 
form. It should be stressed that this is necessary only for general theoretical 
considerations, and will never be done for concrete examples. 

Definition 3.1 A hyperbolic system of first order partial difference equa- 
tions is a system of the form 

SiXk = gk,iix), i G £fc- (3.1) 

for functions '■ with values in Banach spaces X^. For each x^, 

eqs. VJ.1\) are given for z G £fc C {1, . . . , M}, the evolution directions of x^- 
The complement 8^ = {1, . . . , Af} \ 8.^. consists o/ static directions of x^- 

We think of the variable Xk{u) as attached to the elementary cell of 
dimension #Sfc adjacent to the point u G and parallel to Sg^, : 

Cfc = |n + ^ IXiCi : Hi G [0, 

Here, recall, 

■Bs = {-u G Z^'^ : Ui = if 
for an index set S C {1, . . . , M}. 
Definition 3.2 

1) A local Goursat problem for the hyperbolic system iS. 1\) consists of 
finding a solution x^ for all k and for all cells within the elementary 
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cube ofL at the origin from the prescribed values Xfc(O). A global Goursat 
problem consists of finding a solution of ll'^.ll) on subject to the following 
initial data: 

2;fcf'Bs^=^fc, (3.2) 

where : Bg^. ^ are given functions. 

2) The system is called consistent, if the local Goursat problem for 

it is uniquely solvable for arbitrary initial data x[{)). 

The following rather obvious but extremely important statement holds: 

Theorem 3.3 A Goursat problem for a consistent hyperbolic system 
has a unique solution x on all ofL^^. 

Consistency conditions read: 5j5iXk = 6i5jXk for all i / j. Substituting 
eqs. one gets the following equations: 

Sj9k,i{^) = Sigk,j{x), i ^ j, (3.3) 

or gk,i{x + gj{x)) - gk,i{x) = gkj{x + gi{x)) -gkjix), where gi{x) is a vector 
function whose i-th component is equal to ge^i{x), if i G Ei, and is undefined 
otherwise. 

Lemma 3.4 For a consistent system of hyperbolic equations itl.l\) . the func- 
tion gk^i depends on those components xi only for which 8^ C 8^ U {i}. 

Proof. Eqs. (|3.3j) have to hold identically in x. This implies that the 
function gk^i can depend on those components Xi only, for which 6jX£ is 
defined, i.e., for which j £ 6^. As (|3.3|) has to be satisfied for all j € E,k, 
j ^ i, one obtains that for these i there holds 8.k \ {i} C £.£. □ 

It follows from Lemma [3.4l that for any subset 8 C {1, . . . , M}, equations 
of for k with 8fc C 8 and for i £ § form a closed subsystem, in the 
sense that gk^i depend on xi with 8^ C 8 only. 

Definition 3.5 The essential dimension d of the system US. 1\) is given by 

(i= l + max(#8fc). (3.4) 

k 

If (i = M, system has no lower-dimensional hyperbolic subsystems. If 
d < M, then d-dimensional subsystems corresponding to 8 with #8 = d are 
hyperbolic. In this case, consistency of system is a manifestation of a 
very special property of its d-dimensional subsystems, which we suggest to 
treat as the discrete integrability (at least under some further conditions, 
excluding certain non-interesting situations, like trivial evolution in some of 
the directions). Sect. I^will be devoted to giving a solid background for this 
suggestion in the case d = 2. 
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3.2 Discrete approximation in hyperbolic systems 

To handle with approximation results for discrete geometric models, we need 
to introduce small parameters into hyperbolic systems of partial difference 
equations. The definition domain of our functions becomes 



If ej = for some index i, the respective component in "B*^ is replaced by 
R. For instance, if e = (0, ...,0), then = R^^ . So the domains 
posses continuous and discrete directions, with mesh sizes depending on 
the parameters e^. Definitions of translations and difference quotients are 
modified for functions on in an obvious way: 



If €i = 0, then Si is naturally replaced by the partial derivative di. For a 
multi-index a = (ai, . . . , oa/ ), we set 5" = 5"^ . . . . 

The definition of elementary cells Cfc, carrying the variables Xk, is mod- 
ified as follows: 



(so that the cell size shrinks to zero in the directions with = 0). We see how 
the discreteness helps to organize the ideas: in the continuous case, when all 
= 0, all the functions Xk live at points, independently on the dimensions 
#Sfc of their static spaces. In the discrete case, when all ej > 0, one can 
clearly distinguish between functions living on vertices (those without static 
directions), on edges (those with exactly one static direction), on elementary 
squares (those with exactly two static directions), etc. 

Having in mind the limit e — > 0, we will treat only the case when the 
first m < M parameters go to zero, ei = • • • = = e, while the other 
M — m ones remain constant, e^+i = • • • = em = 1- Thus, in this case 
■B^ = (eZ)™ X Z*^-™, and we set S = 2° = x Z^^"™. Assuming that 
the functions gk,i = gl^on the right-hand sides of t^'A.l^ depend on e smoothly 
and have limits as e — > 0, we will study the convergence of solutions x"^ of 
the difference hyperbolic system (|3.1|) towards solutions x^ of the limiting 
differential(-difference) hyperbolic system 



S'^ = eiZ X ••• X eAfZ. 



{Tif){u) = fiu + eiCi), {Sif){u) = —{f{u + eiei) - /(n)). 




diXk 
diXk 



9k,iix), l<i<m, 
gl/x), m + l<i<M. 



(3.5) 
(3.6) 
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Naturally, p.Sf) . H3.6|) describe the respective m-dimensional smooth geom- 
etry with M — m permutable transformations. 

Throughout this section, a smooth function g : T) ^ X is one that is 
infinitely often differentiable on its domain, g G C°°{D). For a compact set 
3C C T), we say that a sequence of smooth functions converges toward a 
smooth function g with the order 0(e) in C°°(3C), if 

11/ - S'llc"^(x) < c^e 

with suitable constants C£ for any £ G N. Convergence in C°°{D) means 
convergence in C°°{%) for all compact sets IK C D. 

Convergence of discrete functions (defined on lattices 'B^ with different 
e) is understood as follows. We say that a family of discrete functions 
x*^ : S"^ — >■ X converges to a function x : — > X with the order 0(e) in 
C°°('B), if for any multi- index a = (ai, . . . , ajv/) there holds 

sup |(5"(x"^ — x){u)\ < Ca€ 

with some constants Cq. The symbol x on the left-hand side is understood 
as a restriction of x to C 23. 

Finally, we mention that we are mainly concerned with local problems, 
so that actually we work with bounded domains of the lattices 'B^, 

■B'(r) = {n G 2' : Ui G [0,r], uj G {0, 1}, 1 < i < m< j < M}. 

Each 23 '^(r) contains only finitely many points (though their number grows 
infinitely with e ^ 0, if r remains fixed). For functions ti"^ defined on a 
bounded lattice domain 23 "^(r) only, the notion of convergence is modified 
in an obvious way: the supremum is taken only over those lattice sites u, 
where the respective difference quotient 6"x^{u) exists. 

Theorem 3.6 Consider a Coursat problem hH. ij)) for a hyperbolic 

system of difference equations. Suppose that: 

i) the discrete system h'J.l\) is consistent for all e > 0; 

a) functions g\^ on the right-hand sides of eqs. converge with the 

order 0(e) in C^{Xk) to smooth functions g^^; 

Hi) the Coursat data converge with the order 0(e) in C°°('B§(-fc)) to 
smooth functions X^. 
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Then, for e > small enough, solution of the Goursat problem exists and 
is unique on S'^(r) for a suitable (e-independent) r > 0; moreover, solutions 
x'' converge to a smooth function x^ with the order 0(e) in C°°(S(r)); this 
function is a unique solution on S(r) of the Goursat problem for \S. ,5)) . 
\3. 61) with the Goursat data X^. 

Proof of this theorem is technical, based on the discrete Gronwah-type 
estimates. Details can be found in .BobMaSl^. □ 

In condition ii), convergence g'j.^^ g^^'m. C°°(Xfc), i.e., on every com- 
pact subset of Xfc, is assumed for simplicity of presentation only. In applica- 
tions, functions g^. ■ are often defined on certain subdomains D^, C X^, with 
the property that is open and dense in X^. In such a case, one requires 
in ii) the convergence 9ki ^ 9ki '^^ C°°(D*^). Then conclusions of Theorem 
1^.61 hold for generic initial data. 

As for condition Hi), smooth data : Sg^ — > Xfc are usually given 
a priori, and discrete data X^ are obtained by restriction to the lattice: 
X^ = . In such a situation, condition Hi) is fulfilled automatically. 

3.3 Conjugate nets 

Discretization of a conjugate net. Recall that a smooth conjugate net 
/ : R™ ^ is determined by the initial data (Qi,2) (see Sect. [TTll . 
while a discrete Q-net : (eZ)"^ is determined by the initial data 

(Qfa) (see Sect. |23I- We now demonstrate how to produce from the data 
(Qi,2) certain discrete data (Q^2)i which will assure the convergence of the 
corresponding discrete Q-nets to the smooth conjugate net. 

Define the discrete curves /"^ f^f by restricting the curves / [3. to the 
lattice points: 

f^u) = f{u), uG-Bl l<i<m. 

Similarly, define the plaquette functions cjjfs|. by restricting Cjj|"s-^ to the 
lattice points: 

Cij{u) = Cij{u), u € "Blj, l<i^ i <m. 

An even better option is to read off the values of cij f at the middlepoints 
of the corresponding plaquettes of 23 ■^•: 

= + f + f Cj), u e 'Blj, 1 < i / J < m. 
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Either choice gives the data (Qi 2) which define an e-dependent family of 
discrete Q-nets : (eZ)"^ — s- M^, called canonical discrete Q-nets corre- 
sponding to the initial data (Qi,2)- 

Theorem 3.7 For some r > 0, the canonical discrete Q-nets : 'B^(r) — > 
converge, as e — > 0, to the unique conjugate net f : 'B{r) with the 

initial data (Qi,2)- Convergence is with the order 0(e) in C°°(23(r)). 

This follows directly from Theorem 13.61 since eqs. H2.1() . (|2.2() and 
(|1.2j) are manifestly hyperbolic (and can be easily re-written in the first 
order form). 

Discretization of a Jonas pair. Recall that a Jonas transform of 
a given conjugate net is determined by the initial data (Ji,2) (see Sect. 
II. We now produce out of these the initial data {jf'2) (see Sect. 12. 1|) 
for an e-dependent family of Jonas transforms of canonical discrete Q-nets 
corresponding to the initial data (Qi,2)- 

Take the point f~^{0) from (Ji). Define the edge functions a^fsfj ^ifsf 
by restricting the functions ajf^-, to the lattice points, or, better, to 

the middlepoints of the corresponding edges of This gives the data set 
(J^2); along with the data (Q^2) produced above this yields in a canonical 
way an e-dependent family of discrete Q-nets F"" : (eZ)*" x {0, 1} —>■ M^, 
which will be called the canonical ones for the initial data (Qi,2)) (Ji,2)- 

Theorem 3.8 The canonical Q-nets (/')+ = F^{-,1) : 'B^(r) con- 
verge to the unique Jonas transform : 'B(r) — s- of f with the initial 
data (Ji,2)- Convergence is with the order 0(e) in C°°(S(r)). 

Again, this follows directly from Theorem 13.61 applied to the hyperbolic 
systems consisting of eqs. (|2.13p - ()2.15() in the discrete case and of eqs. 
(|l.ll|) - (|1.13j) in the smooth case. Note that the discrete equations are 
implicit, and their solvability for e small enough is guaranteed on the set 
{aj 7^ : 1 < j < m}, which is open dense in the phase space. 

3.4 Orthogonal nets 

We start with approximation of a single orthogonal net / : M™' —^ M^. For 
the approximating discrete 0-nets, we have M = m and all e^ = e. In all 
formulas of Sect. 12.21 one has to replace the lattice functions hi, Pij, pkj by 
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ehi, ef3ij, epkj, respectively. Observe that formulas H2.29() . H2.37() become 
i^ji = V,, = (1 - e^f3ijl3jifl^ = 1 + 0(e2), 

Under this re-scaling, eqs. (|T25|) . (|T26|) . (|T27|) and can be put into 

the standard form ()3.1() with the functions on the right-hand sides approx- 
imating, as e — > 0, the corresponding functions in eqs. (|1.6j) - (|1.9j) with the 
order 0(e). 

Nevertheless, formally speaking. Theorem 13.61 cannot be applied to or- 
thogonal nets. The reason for this is that the full system of differential 
equations describing orthogonal nets, consisting of eqs. ()1.6|) ~ H1.9() and the 
constraint (|1.35j) . is non- hyperbolic. Its non-hyperbolicity rests on the fact 
that the constraint H1.35() is not resolved with respect to the derivatives 
difiij. Note, however, that constraint (|1.35|) does not take part in the evolu- 
tion of solutions starting with the data given in the coordinate planes S jj : it 
is satisfied automatically, provided it is fulfilled for the coordinate surfaces 
f\'hij- Therefore, we will obtain a convergence result for orthogonal nets as 
soon as it will be established for coordinate surfaces. 

Discretization of an O-surface. Initial data for a smooth 0-surface 
/ : Si2 ^ are: 

(i) two smooth curves f\'Zi {i = 1;2), intersecting orthogonally at /(O); 

(ii) a smooth function 712 : 23i2 M, whose designated meaning is 712 = 
\{dif3i2-d2(32i). 

Let / I" 3. be the images of the curves /l"^- in the Mobius-geometric model Qq^. 
Let hi = \dif \ and Vi = h^^dif be the metric coefficients and unit tangent 
vectors of the coordinate curves. Choose an initial frame ip{0) € 3ioo such 
that 

/(O) = V-'(0)eoV(0), ^)^(0) = V'-i(0)e,^(0) {i = 1, 2). 

Define the frames ^f) : Hi ^ !Koo of the curves /t^i as the solutions of eqs. 
p. 32)1 for i = 1,2 (considered as ordinary differential equations) with the 
initial value V'(O)- Rotation coefficients of the curves f\'Bi are the functions 
(3ki : Si ^ M defined by the formula (fT^ for i = 1, 2. 

Define the discrete coordinate curves by restricting the functions 

/ f-B, to the lattice points. Let hi = \5if'^\ and vl = {h'i)~^5if'^ be the 
discrete metric coefficients and unit vectors along the discrete curves. Define 
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the frame : "Bj ^ "Koo by iterating the difference equation H2.33() for 
i = 1,2 with the initial condition '^'^(0) = ^(0). Then canonical rotation 
coefficients of the discrete curves are the coefficients pf.- : — > M in 

the expansions 

Finally, let the plaquette function : 'BJg ^ M be obtained by restrict- 
ing 712 to the lattice points (or to the middlepoints of the corresponding 
plaquettes of "Big). 

Thus, we get valid Goursat data for a hyperbolic system of first-order 
difference equations for the variables f^, vf, /if, consisting of eqs. (|2.25j) . 
TTM . (|T77j) . f^T^ and with distinct i,j G {1,2} andl < k < N, 

where the following expressions should be inserted: 

P12 = <^\P\2 - Pk\Pk2-l\2) , P2I = ^2P21 -^{\Y1 PklPk2+ll2) ■ 

k>2 k>2 

The discrete nets f' : 2 ^ defined as solutions of the Goursat problem 
just described are discrete 0-sur faces, since they fulfill the circularity con- 
straint H2.4U() . They will be called canonical discrete 0-surfaces constructed 
from the above initial data. 

Theorem 3.9 For some fixed r > 0, the canonical discrete 0-surfaces 
t ■■ 'Bl2(?') Qo" converge, with the order 0(e) in C~('Bi2(r)), to the 
unique 0-surface f : 'Bi2(r) — > Qq^ with the initial data / [2,; ft = 
and ^{diPi2 — O2P21) = 712 • Edge rotation coefficients and plaquette 
rotation coefficients (3^2, P21 of the discrete 0-surfaces f^ converge to the 
corresponding rotation coefficients Pki of the 0-surface f. 

Proof. First, we show the convergence of the frames, ip^ — > ip, and of the 
rotation coefficients, p^.^ Pki , along the discrete curves This follows 

from two observations. First, vf{0) = Vi{0)-\- ^{diVi){0)-\-O{e'^), so that there 
holds: 

in - 1W{0) = -iei^P{0){dm){0) + 0(e2). 

Second, combining frame equations on two neighboring edges of "B^, one 
finds that everywhere on there holds: 

in - rr^W = -e^V'^(l - rr^Wi = -ee,r{d,Vi) + O(e'). 
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The claim follows by standard methods of the ODE theory. 

Now an application of Theorem IIS . HI shows that functions /"^ : Qq 
converge to the functions f : 'B12 ^ Qo 

which solve the Goursat problem 
for the hyperbolic system of first order differential equations, consisting of 
eqs. (frB|) - ((0|) with distinct i,j G {1,2} and 1 < A: < iV, and 

diPi2 = -^^/3fci/?fc2 +712, 92/321 = -^^/3fci/?fc2 -712- 

k>2 k>2 

The solutions Pki satisfy the orthogonality constraint (|1.24|) and the relation 

UdiPu - 82^21) = 712. □ 

Discretization of an m-dimensional orthogonal net. Given the 
initial data (Oi^2) for an m-dimensional orthogonal net (see Sect. 11.2(1 . we 
can apply the procedure described in the previous paragraph, with an initial 
frame ip{0) G "Koo such that 

/(O) = V~'(0)eoV'(0), ViiO) = ij-\0)e,i^iO) (1 < i < m), 

to produce, in a canonical way, the discrete 0-surfaces f^E, and their 
plaquette rotation coefficients /3|j. Thus, we get data (Of^2) (see Sect. 12.2(1 
for an e-dependent family of discrete 0-nets f'' : (eZ)™ Qq. These nets 
will be called the canonical discrete 0-nets corresponding to the initial data 

(Ol,2). 

Theorem 3.10 The canonical discrete 0-nets f*" : S'^(r) — > converge, 
as e — > 0, to the unique orthogonal net f : 'B{r) —>■ with the initial data 
(01^2) • Convergence is with the order 0(e) in C°°{'B{r)). 

Proof. The data (0^2) yield a well-posed Goursat problem for the hy- 
perbolic system of first-order difference equations for the variables f^, vf, 
hi, 131 consisting of eqs. (tO^l) . (tOHl) . (tOTl) . (tOnl) . The convergence 
of these Goursat data is assured by Theorem IH.9I Now the claim of the 
theorem follows directly from Theorem 13.61 □ 

Discretization of a Ribaucour pair. Given the initial data (Ri,2) 
for a Ribaucour transform of an orthogonal net (see Sect. 11.2(1 . define the 
plaquette rotation coefficients on the "vertical" plaquettes along the 
edges of the coordinate axes 23^ by restricting the corresponding functions 
eOi to lattice points or, alternatively, to middlepoints of the corresponding 
edges of S^: 

I3lji{u) = eOiiu) or eOiiu + f ), G l<i<m. 
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Thus, we get the data (R^2) (^^^ Sect. 12. 2p . which, together with (Of^2)) 
allow us to construct in a canonical way discrete 0-nets : (eZ)™ x {0, 1} 
M^. They will be called the canonical ones corresponding to the initial data 

(Ol,2), (Rl,2). 

Theorem 3.11 The canonical discrete 0-nets (f^)^ = -F^(-,l) : B'^(r) 

converge to the unique Ribaucour transform : 'B(r) ^ of f with 
the initial data (Ri,2)- Convergence is with the order 0(e) in C°°{'B{r)). 

Proof. Define f^/(0) as the unit vector parallel to 5/(0) = /^(O) — /(O), 
and set /i^(0) = |(^/(0)|. These data along with /J^^^ on the coordinate axes, 
added to the previously found ones /^(O), vf, hj, (3fj for I < i, j < m, form 
valid Goursat data for the system fTIKii . fTMi . fTFll . fTMi . 

The circularity constraint (|2.28|) implies that /3|^ = —2{vl,v\^) — e9i on 
all edges of S^. Perform the substitution 

vli = y+0{e), h%j = e+0{e), (31,, = ee,+0{e^), = -2(z;,, y)+0(e) 

in eqs. (|2.26|) . (|2.27j) . I|2.H()|) with one of the indices equal to M. Taking into 
account that in this limit one has 

one sees that the limiting equations coincide with eqs. (|1.26j) . (|1.27|) . (|1.28|) . 
A reference to Theorem 13.61 finishes the proof. □ 

3.5 Moutard nets 

Discretization of an M-net. Given the initial data (Mi 2) for an M-net 
(see Sect. II. 3() . we produce initial data (M^2) fo'^ e-dependent family of 
discrete M-nets with ei = €2 = e. Discrete curves f^ fs^ are obtained from 
the smooth curves /fs. by restricting to the lattice points: 

f%u)=f{u), ue'Bl i = l,2. 

The plaquette function a|2 : (e^)^ — > R is obtained from the function qi2 
restricted to the lattice points: 

al2(n) = l + ie2gi2(n), u e (eZf 

(one could also restrict qi2 to middlepoints of the corresponding plaque- 
ttes). Now discrete M-nets : (eZ)^ —>■ are defined as solutions of the 
difference equation 1)2. 41() with the above data (Mf^g)- 
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Theorem 3.12 Canonical discrete M-nets : 53'^ (r) — s- M converge, as 
e ^ 0, to the unique M-net f : 'B(r) with the initial data (Mi^2)- 

Convergence is with the order 0(e) in C°°('B(r)). 

Proof. Eq. (|2.4H) is manifestly hyperbolic (and can be easily put in the 
first order form). It approximates eq. H1.36() . because it can be re- written 
as 

Sid2f = hlninf + T2f) = qu{f + f <5i/ + f <52/). 
Now Theorem 13.61 can be applied. □ 

Discretization of a Moutard pair. Let the initial data (MTi 2) for a 
Moutard transformation be given. Define the edge variables 6f from the 

i 

functions Pilsi restricted to the lattice points: 

6i(ui,0) = l + epi(ui,0), bl{0,U2) = l + ep2{0,U2), ui e eZ 

(one could restrict Pil's- to the middlepoints of the corresponding edges, as 
well). This gives us the data (MT^g); which canonically generate discrete 
M-nets (/^)+ : {eZf . 

Theorem 3.13 Canonical discrete M-nets {f'^)'^ : 'B'^(r) M.^ converge to 
the unique Moutard transform /+ : S(r) of f with the initial data 

(MTi^2)- Convergence is with the order 0(e) in C°°('B(r)). 

Proof. The system consisting of eqs. (j2.45|) . (|2.46|) is hyperbolic. Upon 
substituting bi = 1 + epi and ai2 = 1 -|- ^e^qu , these equations can be 
re-written as 

<5i/+ + Sif = piif+ - nf), 62f+ - 52f = P2(/+ + r2f), 

and 

1 + eT2Pi ^ 1 + eripa ^ 1 + (e72) g+ ^ 1 

1 + epi l + ep2 l + (e2/2)(?i2 I + e^{qi2 - P1P2) + 0{e^)' 

Clearly, they aproximate, as e ^ 0, eqs. ((07|) ~ (fOH|) and ((0^ - (fOn|) . 
respectively. It remains to apply Theorem 13.61 □ 

3.6 A-surfaces 

Discretization of an A-surface. Initial data (Ai^2) for an A-surface (see 
Sect. II. 4() . are nothing but initial data (Mi^2) for the Lelieuvre normal field 
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n :M? —fW^. In Theorem 13.121 we described the canonical construction of 
the initial data (A^2)> which give a converging family of the discrete Lelieu- 
vre normal fields : (eZ)^ — > M^. Eqs. (|2.53l) define the discrete A-surfaces 
: (eZ)^ — > R^, called the canonical discrete A-surfaces corresponding to 
the initial data (Ai^2)- 

Theorem 3.14 Canonical discrete A-surfaces f": ^"^(r) M'^ converge, as 
e ^ 0, to the unique A-surface f : H^r) with the initial data (Ai^2)- 

Convergence is with the order 0(e) in C°°('B(r)). 

Proof. Eqs. (|2.53)) are hyperbolic and approximate eqs. (|1.43|) . Theorem 
13.61 can be applied to prove the convergence of after the convergence of 
n'' has been already proved. □ 

Discretization of a Weingarten pair. Initial data (Wi,2) for a Wein- 
garten transformation are nothing but initial data (MTi^2) foi" a Moutard 
transformation of the Lelieuvre normal field. We already described, in The- 
orem E^SI the canonical construction of the initial data (Wf^2) ^ con- 
verging family of the discrete Lelieuvre normal fields (n'^)'*' : (eZ)^ — >■ M^. 
Now the transformed A-surfaces (Z*^)"^ : (eZ)^ are obtained by eq. 

(l2TlTl . 

Theorem 3.15 Canonical discrete A-nets (Z^)"'' : 'B^(r) — > converge to 
the unique Weingarten transform : S(r) — > of f with the initial data 
(Wi^2)- Convergence is with the order 0(e) in C°°('B(r)). 

Proof follows by comparing the (identical) formulas H1.45() and 1)2. 51() . after 
the convergence of ra^ and (n*^)"^ has been proved. □ 

3.7 K-sur faces 

Discretization of a K-surface. Given the initial data (K) for a K-surface 
(see Sect. II. 5() . we define initial data (K^) (see Sect. 12. 5|) for an e-dependent 
family of discrete K-surfaces with ei = e2 = e by restricting n to the 
lattice points, as for general A-surfaces. Define discrete M-nets n*^ : (eZ)^ — > 
as solutions of the difference equations (|2.56|) with the initial data (K^). 
Finally, define the discrete K-surfaces : (eZ)^ with the help of the 

discrete Lelievre representation ()2.53() . These will be called the canonical 
discrete K-surfaces corresponding to the initial data (K). 

Theorem 3.16 Canonical discrete K-surfaces f": 'B'^(r) — > converge, as 
e ^ 0, to the unique K-surface f : 'B(r) — > R^ with the initial data (K). 
Convergence is with the order 0(e) in C°°{'B{r)). 
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Proof. We have for n = rf: 



{n,Tin + T2n) 2 + e{n,5in + d2n) 



1 + (rin,T2n) 2 + e(n, 6in + 52n) + e2(5in, (52n) ' 



Since {n,5in) = 0(e), we find that 



al2 = l-^e^{6in,62n) + 0{e^). 



Comparing this with eq. (|1.48|) . we see that Theorem 13.61 can be apphed in 
order to prove approximation of n : ^ by n"^ : (eZ)^ — > S^. Finally, 
approximation of / by follows exactly as for general A-surfaces. □ 

Discretization of a Backlund pair. Let the initial data (B) for a 
Backlund transformation of a given K-surface /, i.e. the point n~^{0), be 
given. Take it as the initial data (B^) for the discrete Backlund transforma- 
tions (Z^)"'' : (eZ)^ of the family of discrete K-surfaces constructed 
in Theorem 13. 161 

Theorem 3.17 Canonical discrete K-surfaces (/'^)~'" : 23'^(r) M'^ converge 
to the unique Backlund transform : 'B(r) — > M'^ of f with the initial data 
(B). Convergence is with the order 0(e) in C°°('B(r)). 

Proof. In eqs. (^37]) . we have: 



{n, Tin — n"*") 



1 + e 



e ; — - 

1 — (n, n+) 



{6in, n"*") 



+ 0(e2), 



1 — (rin, n+) 
{n,T2n + n+) 



= 1 - 



+ 0(e2). 



1 + (r2n, n+) 



1 + (n,n+) 



Comparing this with eqs. ()1.46|) - H1.47|) and applying Theorem 13. 61 we prove 
approximation of the Gauss maps. □ 
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Chapter 4 

Consistency as integrability 



4.1 From 3D consistency of discrete 2D systems to 
Backlund transformations and zero curvature 
representations 

Equations of a discrete 2D hyperbolic system are associated to elementary 
squares of 1? . Such a system may contain fields on vertices and/or on edges 
of 1?. In our differential-geometric considerations we encountered two types 
of such systems: the cross-ratio system, with fields on the vertices and the 
parameters on the edges of Z^, and discrete Lorentz-harmonic maps into 
quadrics, with fields on the vertices of 1? . One can imagine further types of 
discrete hyperbolic 2D systems, for instance, those with fields on edges only. 
We will mainly consider systems of the cross-ratio type, and will discuss 
other types briefly. 

2D systems with fields on vertices and with labelled edges. A 

typical representative of this class of equations is the cross-ratio system: 

?(/,/l,/l2,/2) = -. (4.1) 
Oil 

The variables j : 1? ^ e£(M^) are associated to the vertices of Z^, while 
the variables : — > M satisfy the labelling property 

b^ax = 0, ha2 = 0. (4.2) 

Therefore a, play the role of parameters naturally associated to edges of 
T? parallel to and constant along the strips in the complementary direc- 
tion. (As a disclaimer, we mention that further edge variables have to be 
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introduced in order to put this equation in the first-order form.) One more 
example of such a system with vertex variables and parameters sitting on 
edges and having the labelling property is given by the Hirota equation, 



fl2 
f 



"2/1 - al/2' 



(4.3) 



see |BobPll IBobMaSl] for the geometric interpretation of this system in 
terms of discrete K-surfaces. A general system of this class consists of equa- 
tions 

Q(/,/i,/i2,/2;ai,a2) = 0, (4.4) 

see Fig. 14.11 In the simplest situation one has complex fields / : — > C, 
and complex parameters ai on the edges of parallel to Sj, satisfying the 
labelling condition (|4.2I) . We require that eq. 1)4. 4|) be uniquely solvable for 
any one of its arguments /, /i, /2, /12 G C. Therefore, the solutions have to 
be fractional- linear in each of their arguments. This naturally leads to the 
following assumption. 

Linearity. The function Q{x, u, y, v; a, (3) is a polynomial of degree 1 in 
each argument (affine linear): 

Q{x, u, y, V] a, (3) = ai{a, (3)xuyv + • • • + ai6(a, /?). 

Note that for the cross-ratio equation (|4.1j) with complex-valued argu- 
ments, Q = P{x — u){y — v) — a{u — y){v — x) , while for the Hirota equation 
1)4. 3|1 . Q = a{xu + yv) — (3{xv + yu). 




O 

/ "1 /i 
Figure 4.1: Elementary quadrilateral; fields on vertices 



Actually, this setup admits an important generalization: elementary 
quadrilaterals carrying eqs. (|4.4)) can be attached one to another with the 
combinatorics more complicated then that of T?. 
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Definition 4.1 A quad-graph T) is a strongly regular polytopal cell decom- 
position of a surface with all quadrilateral faces. 

We denote by V{D), E{T>), F{D) the sets of vertices, of edges and of faces 
of T>, respectively. We consider eq. ()4.4|) for fields / : V{T>) — > C, with 
a : E{T>) C being a labelling of edges of T>, i.e., a function taking equal 
values on any pair of opposite edges of any quadrilateral from F{T)). In the 
context of equations on quad-graphs, there are no distinguished coordinate 
directions, nevertheless it will be convenient to continue to use notations 
of eq. (|4.4j) . with the understanding that indices are used locally (within 
one quadrilateral), and do not stand for shifts into the globally defined 
coordinate directions. So, /, /i, /12, /2 can be any cyclic enumeration of the 
vertices of an elementary quadrilateral. Eq. (|4.4j) should not depend on the 
enumeration of vertices, therefore the following assumption is natural when 
considering equations on general quad-graphs. 

Symmetry. The function Q has the symmetry properties 

Q{x, u, y, v; a, (3) = eQ{x, v, y, u; (3, a) = aQ{u, y, v, x; (3, a), e, a = ±1. 

(For the complex cross-ratio equation and for the Hirota equation the func- 
tions Q given above possess this property with e = a = —1). 

Assume now that eq. (|4.4j) possesses the property of 3D consistency. 
Recall that this means that this equation can be consistently imposed on 
all 2D faces of a combinatorial cube on Fig. 14.21 with the parallel edges of 
the i-th direction all carrying the parameter Oj. More precisely, given the 
initial data f,fi, one determines fij in virtue of eq. (|4.4|) . and then three 
equations for /123 (for three faces adjacent to this vertex) lead to identical 
results. This is the case for the complex cross-ratio equation, with 



/123 



(ai - a2)/i/2 + («2 - az)f2fz + (^3 - «i)/3/i 
(02 - ai)/3 + ("3 - "2)/! + («! - a3)/2 



(4.5) 



as well as for the Hirota equation, with 



/123 



azjoil - ai)/i/2 + aijal - ai)/2/3 + «2(ai - aDfsfi 
asiaj - af )/3 + ai(a| - a|)/i + 02(0^ - a§)/2 



(4.6) 



We will demonstrate that this condition automatically leads to two basic 
structures associated in the soliton theory with integrability: Backlund 
transformations and zero curvature representation. 
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Figure 4.2: 3D consistency of 2D systems 



Theorem 4.2 For any solution f : V{D) C of a 3D consistent equation 
\4---i\) on the quad- graph D, there is a two-parameter family of solutions 
: V{Ti) C of the same equation, satisfying 



QifJi,ft,f^;ai,^) = o 



(4.7) 



for all edges {f,fi) G E{T>). Such a solution is called a Bdcklund 
transform of f , and is determined by its value at one vertex of T) and by 
the parameter A. 



Proof. Extend the planar quad-graph D into the third dimension. Formally 
speaking, we consider the second copy !D+ of D and add edges connecting 
each vertex / G V{T>) with its copy /"*" G V{1)~^). (We slightly abuse the 
notations here, by using the same letter / for vertices of the quad-graph 
and for the fields assigned to these vertices.) On this way we obtain a "3D 
quad-graph" D, with the set of vertices 

V{U) = y(D) UF(D+), 

with the set of edges 

E{B) = E{V) U E{T)^) U {(/, /+):/€ V{V)}, 
and with the set of faces 



F(D) = F(a3) U F{D+) U {(/, A, /+, /+) : /, /i G V{D)]. 

Elementary building blocks of D are combinatorial cubes (/, /i, /i2> /2j 
f~^ , fi , fi2, f2)^ shown on Fig. 14.31 The labelling on E{Y)) is defined 
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in the natural way: each edge {f~^,fi~) € E^D^) carries the same label ai 
as its counterpart {f,fi) G -E'(I')) while all "vertical" edges {f,f^) carry 
one and the same label A. Clearly, the content of Fig. 14.31 is the same as 
of Fig. 14.21 up to notations. Now, a solution : V{D^) ^ C on the first 
flour of D is well defined due to the 3D consistency, and is determined by 
its value at one vertex of and by A. We can assume that is defined 
on V{D) rather than on V{D~^), since these two sets are in a one-to-one 
correspondence. □ 




Figure 4.3: Elementary cube of D 



Theorem 4.3 A 3D consistent system i4.4\ ) on the quad-graph T) admits a 
zero curvature representation with spectral parameter dependent 2x2 
matrices: there exist matrices attached to directed edges ofD, 

L{t, a(c); A) : ^(D) ^ GL2(C)[A], (4.8) 

such that for any quadrilateral face (/, /i, /12, /2) G -^(^) there holds 

-^^(/i2,/i,a2;A)L(/i,/,ai; A) = L(/i2, /2, ai; A)L(/2, /, 02; A), (4.9) 

identically in A. 

Proof. Due to the linearity assumption, equations 1)4. 7(1 can be solved 
for /.^ in terms of a Mobius (fractional-linear) transformation of f~^ with 
coefficients depending on /, ff. 



ft = L{f\,f,a,;X) [/+], 



(4.10) 
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Here we use the standard matrix notation for the action of Mobius trans- 
formations: 

L[z] = {az + b){cz + d)-^, where ^ = (^^ (^-H) 

Now 3D consistency for yields that for any there holds: 

L{fi2, /i, 02; A)L(/i, /, ai; A) [/+] = L{fu, f2, ai; A)L(/2, /, as; A) [/+]. 

(4.12) 

Therefore, eq. (|4.9() holds at least projectively, i.e., up to a scalar factor. 
A normalization of determinants of L (or any other suitable normalization) 
allows one to achieve that eq. (|4.9jl holds in the usual sense. □ 

As an example, we derive a zero curvature representation for the complex 
cross-ratio equation 1)4.1(1 . It will be convenient to re-define the spectral 
parameter in this case by A 1-^ A^^, so that equations on the vertical faces 
of Fig. 14.31 read: 

i — — Act? . 

This gives the Mobius transformation (|4.1U() with 

L(/„/,a„A)=/ + -^( -/J^) (4.13) 

The determinant of this matrix is constant (equal to 1 — Aa^), therefore 
no further normalization is required. A more usual form of the transition 
matrices of the zero-curvature representation for the complex cross-ratio 
equation is obtained by the gauge transformation 

L(/„/,a,;A) ^ A-\n)Linj,af, X)A{f), A{f) = 

which leads to the matrices 

^^••^■"•^^' = (a„.(/-/.)-' ^^O' 

These matrices ()4.14|) are interpreted as matrices of the Mobius transforma- 
tions of the shifted quantities: 

ft-h = mj,af,X) [/+-/]. 
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Thus, we have seen that 3D consistency of a 2D equation on a quad-graph 
with complex fields at vertices and with labelled edges implies existence of 
Backlund transformations and of the zero curvature representation. This is 
not a pure existence statement but rather a construction: both attributes 
can be derived in a systematic way starting with no more information than 
the equation itself, they are in a sense encoded in the equation provided it 
is 3D consistent. 

Non-commutative equations. The truth contained in the last para- 
graph is by no means restricted to the situation for which it was demon- 
strated. For instance, in |BobSu2] it was extended to equations with fields 
on vertices taking values in some associative but non-commutative algebra 
A with unit over a field %, and with edge labels with values in %. The 
transition matrices of the zero curvature representation are in this case 2x2 
matrices with entries from A. They act on A according to eq. (|4.11|) . 
where now the order of the factors is essential. The examples worked out in 
|BobSu2] include the non-commutative analogs of the Hirota equation (|4.3|1 
and of the cross-ratio equation ()4.1() (recall that the latter equation with 
fields in yi = C£(M^) and with parameters ai from IK = M governs discrete 
isothermic surfaces in M^). In particular: 

Theorem 4.4 The cross-ratio equation in an associative algebra A is 3D 
consistent. It possesses a zero curvature representation with the transition 
matrices \4-M^ ! where the inversion is treated in A. 

Proof. The transition matrices are obtained from the equations them- 
selves, in a full analogy with the scalar case. Verification of the zero cur- 
vature representation thus obtained is a matter of a direct computation. 
The non-commutative 3D consistency is a consequence of the zero curvature 
representation, see |BobSu2j . □ 

We consider here one more equation of this kind: 



with the vertex variables / taking values in A and with the edge labels a 
from %. Clearly, in the case when A = C^(M^), solutions of this equation 
are special T-nets in M^: 



In the latter vector form this equation was introduced in |Schlj under the 
name of discrete Calapso equation. 



(/l2 - /)(/2 - /l) =02-01 



(4.15) 




(4.16) 
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Theorem 4.5 Equation j-^^. j5| ) in an associative algebra A is 3D consistent. 
It possesses a zero curvature representation with the transition matrices 

^/i,/,«i;A)= ( ( ^""j/^^)- (4-17) 
Proof. Equations H4.15() on the vertical faces of Fig. 14.31 read: 

= / + (A - - hr' = Lin, f, A) [/+]. 

This gives the transition matrices, which can then be used to prove the 3D 
consistency, cf. |BobSu2j . □ 

2D systems with fields on vertices. Regarding T-nets in quadrics, 
we encounter an equation with vertex variables / : ^ Q = {/ € : 
(/> /) = ^^}, aiid with no edge variables: 

f f (f f \ (/, /1-/2) 2(/, /1-/2) 1Q^ 

For the quantities Oj = 2(/, /j) the labelling property ()4.2|) is fulfilled, but 
now they are functions of the vertex variables / rather than parameters of 
equation. Comparing eq. 1)4. 18() with eq. 1)4. 16() . we see that the former can 
be regarded as a particular instance of the latter. 

Theorem 4.6 Equation \4.1t^ , describing T-nets in quadrics, is 3D consis- 
tent. It possesses a zero curvature representation with the transition matrices 
with entries from S^(M^).- 

mj;X)=(^{ ^tjf)- (4-19) 

Proof. 3D consistency has been proven geometrically in Theorem 12.221 As 
for the transition matrices, we can take those from eq. (|4.17)) with 

ai = 2{f,fi) = -ffi-fif. 

Note the geometrical meaning of the spectral parameter: A = 2{f,f^) for 
the Backlund transformation from which the transition matrices are 
constructed. □ 

2D systems with fields on edges. Another large class of 2D systems 
on quad-graphs build those with fields assigned to the edges, see Fig. 14.41 
In this situation it is natural to assume that each elementary quadrilateral 
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Figure 4.4: Elementary quadrilateral; fields on edges 

carries a map i?:XxXi— >XxX, with X being the set where the fields x, y 
take values, so that (x2,yi) = R{x,y). The 3D consistency of such maps 
can be encoded in the formula 

i?23 o i?i3 o i?i2 = Rl2 o Rl3 o i?23 , (4.20) 

where each Rij : X^ i— > X^ acts as the map R on the factors i,j of the 
cartesian product X^, and acts identically on the third factor. This equa- 
tion should be understood as follows. The fields x, y are supposed to be 
attached to the edges parallel to the 1st and the 2nd coordinate axes, re- 
spectively. Additionally, consider the fields z attached to the edges parallel 
to the 3rd coordinate axis. Fig. 14.51 illustrates eq. 1)4. 2U() . its left-hand side 
corresponding to the chain of maps along the three rear faces of the cube: 

(x2,yi), {x2,z) ^ {x23,zi), (yi,zi) ^ {yi3,zi2), 

and the right-hand side corresponding to the chain of maps along the three 
front faces of the cube: 

(y, z) {ys, Z2), {x, Z2) ^ (x3, 212), (xa, y^) ^ (X23, 2/13)- 

So, eq. (|4.2flj) assures that two ways of obtaining (X23, yis, -^12) from the 
initial data {x,y,z) lead to the same results. Maps with the property (|4.2U|) 
were introduced by Drinfeld |Drj under the name of set-theoretical solutions 
of the Yang-Baxter equation, an alternative name Yang-Baxter maps was 
proposed by Veselov in the recent study p^, see also references therein. The 
notion of the zero curvature representation makes perfect sense for Yang- 
Baxter maps, and is expressed as 

L{yi;X)L{x2; A) = L(x; A)L(y; A). 

There is a construction of zero curvature representations for Yang-Baxter 
maps with no more input information than the maps themselves |SuVj , close 
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Figure 4.5: Yang-Baxter relation 



in spirit to Theorem 14.31 Consider parameter- dependent Yang-Baxter maps 
R{a,f5), with the parameters a,(3 £ C assigned to the same edges of the 
quadrilateral on Fig. 14.41 as the fields x, y, opposite edges carrying the same 
parameters. Although this can be considered as a particular case of the 
general notion, by introducing X = X x C and R{x, a; y, (3) = R{a, P){x, y), 
it is convenient for us to keep the parameter separately. Thus, on Fig. 14.51 
all edges parallel to the x (resp. y, z) axis, carry the parameter a (resp. 
f3,^), and the corresponding version of the Yang-Baxter relation reads: 

i?23(A7)^i3(«,7)^i2(a,/3) = i?i2(«,/3)i2i3(«,7)^23(/3,7)- (4.21) 

Theorem 4.7 Suppose that there is an effective action of the linear group 
G = GLjv(C) on the set X (i.e., A £ G acts identically on X only if A = I), 
and that the Yang-Baxter map R{a, j3) has the following special form: 

X2 = B{y,[3,a)[xl yi = ^(x, a, /3)[y]. (4.22) 

Here A,B :XxCxC^G are some matrix-valued functions on X depending 
on parameters a and [3 and A[x\ denotes the action of A €z G on x X. 
Then, whenever {x2,yi) = R{a, P){x,y), there holds: 

A{yi,P,X)A{x2,a,X) = A{x,a, X)A{y, X), (4.23) 
B{x2,a,X)B{yi,P,X) = B{y, P, X)B{x,a, X). (4.24) 

In other words, both A{x,a,X) and B~^{x, a, X) form zero curvature repre- 
sentations for R. 
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4.2 Classification of 3D consistent systems 

The notion of 3D consistency proves useful also in various classification prob- 
lems of the integrable systems theory. We give here a presentation of results 
of [AdBSlj on the classification of integrable quad-graph equations (|4.4I) 
with complex fields on vertices and a complex- valued labelling of edges, and 
those of |AdBS2j on the classification of integrable equations with complex 
fields on edges (Yang-Baxter maps). 

4.2.1 2D systems with fields on vertices and labelled edges 

The classification of integrable equations of the type (|4.4j) is performed in 
|AdBSlj under the linearity and symmetry assumptions of the previous sub- 
section, and one additional assumption. This latter one is less natural but it 
is fulfilled for the vast majority of interesting examples, including the com- 
plex cross-ratio and Hirota equations. An attentive look at eqs. (|4.5() . (|4.6p 
tells that for the cross-ratio and for the Hirota equations there holds: 

Tetrahedron property. The value /123, existing due to 3D consistency, 
depends on /i, /2 and /s, but not on /. 

Thus, the fields /i, /2, /s and /123 (vertices of a white tetrahedron on 
Fig. 14. 3|) are related by a well-defined equation. Of course, due to symme- 
try the same holds for the fields /, /12, /23 and /13 (vertices of the black 
tetrahedron) . 

Theorem 4.8 The 3D consistent quad-graph equations i4-4\ ) with the lin- 
earity, symmetry, and tetrahedron properties are exhausted, up to common 
Mohius transformations of the variables f and point transformations of the 
parameters a, by the following three lists Q, H, A (we use the abbreviations 
X = f, u = fi, V = f2, y = fi2, a = ai, /3 = 02^. 
List Q: 

(Ql) a{x - v){u -y) - f3{x -u){v -y) + 5'^aP{a - /5) = 0, 

(Q2) a{x — v){u — y) — (3{x — u){v — y) + aj3{a — I3){x + y + u + v) 

-af3{a - /3)(a2 - a(3 + (3^) = 0, 



(Q3) sin(a)(3;M + vy) — sm{l3){xv + uy) — sin(a — (3){xy + uv) 

sin(a — /?) sin(a) sin(/3) = 0, 
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(Q4) sn(a)(xM + vy) — sn{P){xv + uy) — sn{a — (3){xy + uv) 

+sn(a — /3)sn(a)sn(/3)(l + k'^xyuv) = 0, 

where sn(a) = sii(a; k) is the Jacobi elliptic function with the modulus k. 
List H: 

(HI) {x-y){u-v)+l3-a = Q, 

(H2) {x - y){u - v) + {15 - a){x + y + u + v) + 15'^ - = 
(H3) a{xu + vy) - P{xv + uy) + d{a^ - 0^) = 0. 
List A: 

(Al) a{x + v){u + y) - f5{x + u){v + y)- 6^aP{a - /?) = 0, 

(A2) sm(a)(xf + uy) — sm{P){xu + vy) — sm(a — + xyuv) = 0. 

Remarks. 

1) Parameter 6 in eqs. (Ql), (Q3), (H3), (Al) can be scaled away, so 
that one can assume without loss of generality that 6 = or 6 = 1. 

2) If one extends the transformation group of equations by allowing 
Mobius transformations to act on the variables x, y differently than on n, v 
(white and black subgraphs of a bipartite quad-graph), then eq. (Al) is 
related to (Ql) by the change u — > —u, v — > —v, and eq. (A2) is related to 
(Q3)5=o by the change u ^ 1/u, v ^ 1/v. So, really independent equations 
are given by the lists Q and H. 

3) Note that the above lists contain the complex versions of the cross- 
ratio equation (Ql)5=0) of the Hirota equation (H3)5=o, and of the discrete 
Calapso equation (HI). The latter two equations are, probably, the oldest 
ones in our lists, they can be found in the work of Hirota |Hij . Eqs. (Ql) and 
(Q3)5=o S° back to |QNCL| . Eq. (Q4) was found in |Adj (in the Weierstrass 
normalization of an elliptic curve; the observation that the formulas become 
much nicer in the Jacobi normalization is due to J. Hietarinta). Eqs. (Q2), 
(Q3)5=i, (H2) and (H3)5=i seem to have appeared explicitly for the first 
time in |AdBSlj . 

The classification in Theorem 14.81 is performed modulo simultaneous 
Mobius transformations of all the variables. Applying this theorem, one 
would like to identify a 2D equation at hand with one of the equations of 
the lists Q, H, A. This task is greatly simplified by using some further infor- 
mation. For all 3D consistent equations with the linearity, symmetry and 
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the tetrahedron properties, the following holds. The symmetric biquadratic 
polynomial 

g{x,u;a,(3) = QQyy - QyQy (4-25) 
admits a representation 

g{x,u;a,(3) = k{a, (3)h{x,u;a), (4-26) 

where the factor k is antisymmetric, k{f3, a) = —k{a, (3), and the coefficients 
of the polynomial h{x, u; a) depend on a single parameter a in such a way 
that its discriminant 

r{x) = hl- 2hKu (4.27) 

does not depend on a at all. Thus, we can characterize 3D consistent equa- 
tions Q = of Theorem 14.81 by much less complicated objects, namely the 
biquadratic symmetric polynomials h{x, u; a) of two variables and the quar- 
tic polynomials r{x) of one variable. The action of the simultaneous Mobius 
transformations on the variables x {ax + b) /{cx + d) transforms the poly- 
nomials h, r as follows: 

,/ X / a9/ + ^ au + b 
n{x, u; a) {cx + d) {cu + d) n\ -, -; a 

Y CX I CL CUj I (Ji 

/X / j-,A fax + b 

r(x] I— > (cx + a) r[ 

^ ^ ^ ^ \cx + d 

Using such transformations one can put the polynomial r{x) into one of the 
following canonical forms, depending on the distribution of its zeroes: either 
r{x) = 0, or r{x) = 1 (one quadruple zero), or r{x) = x (one simple zero 
and one triple zero), or r{x) = (two double zeroes), or r{x) = l — x"^ (two 
simple zeroes and one double zero), or, finally, r{x) = (1 — x^){\ — k'^x^) 
with /c^ 7^ 1 (four simple zeroes). So, the first step in identifying a 3D 
consistent equation is computing the polynomial r{x) and putting it by a 
Mobius transformation into one of the canonical forms above. After that, one 
has to identify the polynomial h{x, u; a) with one of those corresponding to 
equations of our lists. This might require a further Mobius transformation in 
the cases r{x) = 0, r{x) = 1, and r{x) = x^ . In other three cases, r(x) = x, 
r(x) = 1 — x^, and r(x) = (1 — x^)(l — A;^x^) (corresponding to equations 
(Q2), (Q3) and (Q4)), the polynomials /i(x, u; a) are uniquely determined by 
r(x), as a one-parameter family of symmetric biquadratic polynomials with 
the discriminant r(x). In the non-degenerate case r(x) = (1 — x'^)(l — k'^x'^) 
this family is given by 

h{x,u;a) = — ^—-(x'^+u^-2cn{a)dn{a)xu-sn'^{a){l+k'^x'^u^)). (4.28) 
2sn(a) V / 
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One can recognize this polynomial as the addition theorem for the elliptic 
function sn(x; k); more precisely, h{x, u;a)=0 if and only if x = sn(^; k) and 
u = sn(7y; k) with — rj = ita. This is the origin of the elliptic parametriza- 
tion of the equation (Q4). Eqs. (Q1)-(Q3) are obtained from (Q4) by 
degenerations of an elliptic curve into rational ones. Similarly, eqs. (Hl)- 
(H2) are limiting cases of (H3). One could be tempted to spoil down the lists 
Q, H to one item each. However, the limit procedures necessary for that are 
outside of our group of admissible (Mobius) transformations, and, on the 
other hand, in many situations the "degenerate" equations (Q1)-(Q3) and 
(H1)-(H2) are of interest for themselves. This resembles the situation with 
the list of six Painleve equations and the coalescences connecting them. 

4.2.2 2D systems with fields on edges 

Consider Yang-Baxter maps R : X x X — > X x X, {x,y) i— > {u,v) in the 
following special framework. Suppose that X is an irreducible algebraic 
variety, and that R is a birational automorphism of X x X. Thus, the 
birational map R~^ : X x X — > X x X, {u,v) i—^ {x, y) is defined. This is 
depicted on the left square on Fig. 14.61 Further, a non-degeneracy condition 
is imposed on R: rational maps u{-,y) : X ^ X and v{x, •) : X — > X should 
be well defined for generic x, resp. y. In other words, birational maps 
^ : X X X — > X X X, (x, t;) I— > (u, y) and R^^ : X x X ^ X x X, {u,y) i—>- 
{x,v), called companion maps to R, should be defined. This requirement is 
visualized on the right square on Fig. 14.61 Birational maps R satisfying this 
condition are called quadrirational |AdBS2j . 



u u 




X X 



Figure 4.6: A map on X x X, its inverse and its companions 



It turns out to be possible to classify all quadrirational maps in the case 
X = CP^; we give a short presentation of the corresponding results |AdBS2j . 
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Birational isomorphisms of CP^ x CP^ are with necessity of the form 

^ _ aiy)x + bjy) ^ ^ _ A{x)y + B{x) ^ ^^^^^^ 
c{y)x + d{y) ^ C{x)y + D{x)^ 

where a{y), . . . , d{y) are polynomials in y, and A{x), . . . , D{x) are polyno- 
mials in X. For quadrirational maps, the degrees of all these polynomials are 
< 2. Dependent on the highest degree of the coefficients of each fraction in 
(jUSni), we say that the map is [1:1], [1:2], [2:1], or [2:2]. The most rich and 
interesting subclass is [2:2]. A necessary condition for a map of this subclass 
to be quadrirational is that quartic polynomials 5{y) = a{y)d{y) — b{y)c{y) 
and A(x) = A{x)D{x) — B{x)C{x) are simultaneously of one of the following 
five types: they have either (I) four simple roots, or (II) two simple and two 
double roots, or (III) two double roots, or (IV) one simple and one triple 
root, or, finally, (V) one quadruple root. 

Theorem 4.9 Any quadrirational map [2:2] on CP"*^ x CP^ is equivalent, 
under some change of variables acting by Mobius transformations indepen- 
dently on each field x, y, u, v, to exactly one of the following five maps: 

{l-/3)x + P-a+{a-l)y 



(Ri): 


u = 


-- ayP, 


V = PxP, P 


(Rll): 


u - 


a 


X 


(Rlll): 


u 


= '-P, 
a 


X 

»=^p. p 


(Riv): 


u 


= yP, 


V = xP, P = 


(Rv): 


u -- 


= y + P, 


V = X + P, 



/3(1 — a)x + (q — (3)yx + a{(3 — l)y ' 
ax — f3y + P — a 

x-y 
ax — (3y 
x-y 

1 + ^, 
x-y 

x-y' 

with some suitable constants a, /3. 

Each one of these maps is an involution and coincides with its companion 
maps, so that all four arrows on Fig. 14. 61 are described by the same formulas. 
Note also that these maps come with the intrisically built-in parameters a, (5. 
Neither their existence nor a concrete dependence on parameters is pre- 
supposed in Theorem 14.91 A geometric interpretation of these parameters 
can be given in terms of singularities of the map; it turns out that the 
parameter a is naturally assigned to the edges x, u, while [3 is naturally 
assigned to the edges y, v. 
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The most remarkable fact about the maps (Ri)-(Rv) is their 3D con- 
sistency. For T = I, II, III, IV or V, denote the corresponding map Rj of 
Theorem 14.91 bv R'j{a,(5), indicating the parameters exphcitly. Moreover, 
for any ai,a2,a3 € C, denote by Rij = the corresponding maps 

acting nontrivially on the z-th and the j-th factors of (CP-*^)^. 

Theorem 4.10 For any T = I, II, III, IV or V, the maps Rij satisfy the 
Yang-Baxter equation \4-2U^ . 

Actuahy, the 3D consistency of quadrirational maps on CP^ x CP^ takes 
place not only for the normal forms R<j{ai,aj) but under much more general 
circumstances. The only condition for quadrirational maps [2:2] consists in 
matching singularities along all edges of the cube (see details in |AdBS2]). 
Similar statements hold also for quadrirational maps [1:1] and [1:2], so that 
in the case X = CP^ the properties of being quadrirational and of satisfying 
the Yang-Baxter equation (|4.20|) are related very closely. 

The maps R'j of Theorem 14 . 91 admit a very nice geometric interpretation. 
Consider a pair of nondegenerate conies Qi, Q2 on the plane CP'^, so that 
both Qi are irreducible algebraic varieties isomorphic to CP^. Take X E 
Qi, Y € Q2, and let £ = XY be the hue through X,Y (well-defined if 
X ^ Y). Generically, the line i intersects Qi at one further point U X, 
and intersects Q2 at one further point V ^ Y. This defines a map 3' : 
{X,Y) ^ {U, V), see Fig. |0 for the picture. The map J : Qi x Q2 ^ 
Qi X Q2 is quadrirational, is an involution and moreover coincides with its 
both companions. This follows immediately from the fact that, knowing 
one root of a quadratic equation, the second one is a rational function of 
the input data. Intersection points X G Qi fl (^2 correspond to the singular 
points {X, X) of the map J'. 

Generically, two conies intersect in four points, however, degeneracies 
can happen. Five possible types I — V of intersection of two conies are 
described in detail in |Berj : 

I: four simple intersection points; 

II: two simple intersection points and one point of tangency; 
III: two points of tangency; 

IV: one simple intersection point and one point of the second order tan- 
gency; 

V: one point of the third order tangency. 
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Figure 4.7: A quadrirational map on a pair of conies 



Using rational parametrizations of the conies: 

CP^ 9 X ^ X{x) eQiC CP2, resp. CF^ B y ^ Y{y) £ Q2 C CP^ 
it is easy to see that 3" pulls back to the map F : {x,y) 1— > {x2,yi) which 



is quadrirational on CP 



One shows by a direct computation that 



the maps F for the above five situations are exactly the five maps listed in 
Theorem 14.91 Now, we obtain the following geometric interpretation of the 
statement of Theorem 14.1 













" / z 1 










\^23 





















Figure 4.8: 3D consistency at a linear pencil of conies 
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Theorem 4.11 Let Qi, i = 1,2,3 be three non- degenerate members of a 

linear pencil of conies. Let X ^ Qi, Y & Q2 and Z £ he arbitrary points 
on these conies. Define the rn,aps 3'ij as above, corresponding to the pair 
of conies {Qi,Qj). Set {X2,Yi) = 3^u{X,Y), {X^,Z^) = 'Sxz{X,Z), and 
{Y3,Z2) = 3^23{Y,Z). Then 

X23 = X3Ysnx^eQi, Y^3 = x^3nY[zleQ2, 
Zi2 = YiZ^r]X^eQ3. 

4.3 Integrability of discrete Laplace type equations 

There exist discrete equations on graphs which are not covered by the theory 
we developed up to now. 

Definition 4.12 Let 9 be a graph, with the set of vertices F(S) and the set 
of edges -E(S)- Discrete Laplace type equations on the graph S for a 
function f : V{S) C read: 

Yl Hfi^o), fix); u{xo, x)) = 0. (4.30) 
x^xo 

There is one equation for every vertex xq G ^(S); the summation is extended 
over the set of vertices of S connected to xq by an edge; the function (j) 
depends on some parameters v : -E(S) — > C. 

The classical discrete Laplace equations on S are a particular case of this 
definition: 

Yu{xo,x){f{x)-f{xo))=0, (4.31) 

X^Xo 

with some weights u : E{S) — > M+ attached to (undirected) edges of 9- We 
use the notation star(a:o) = star(xo; S) for the set of edges of S incident to 
Xq. In view of this, discrete Laplace type equations live on stars of S. 

The notion of integrability of discrete Laplace type equations is not well 
established yet. We propose here a construction which works under addi- 
tional assumptions about the graph S and is based on the relation to 3D 
consistent quad-graph equations (with fields on vertices). The additional 
structural assumption on S is that it comes from a strongly regular polytopal 
cell decomposition of an oriented surface. To any such S there corresponds 
canonically a dual cell decomposition S* (it is only defined up to isotopy, 
but can be fixed uniquely with the help of the Voronoi/Delaunay construc- 
tion). If one assigns a direction to an edge e € E{3), then it will be assumed 
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that the dual edge e* G E{S*) is also directed, in a way consistent with the 
orientation of the underlying surface, namely so that the pair (e, e*) is posi- 
tively oriented at its crossing point. This orientation convention implies that 
e** = — e. The double D is a quad-graph, constructed from S and its dual 
g* as follows. The set of vertices of the double V is V{V) = V{9) U ^(3*). 
Each pair of dual edges, say e = {xo,xi) € E{S) and e* = {yo,yi) G -^(S*), 
defines a quadrilateral {xo,yo,xi,yi). These quadrilaterals constitute the 
faces of the cell decomposition (quad-graph) D. Thus, a star of a vertex in 
xo S ^(S) generates a flower of adjacent quadrilaterals from F{T>) around 
xo, see Figs. 1121 Eni 

This construction can be reversed. Start with a bipartite quad-graph 
D, whose vertices V{D) are decomposed into two complementary halves 
("black" and "white" vertices), such that the ends of each edge from E(D) 
are of different colours. For instance, any quad-graph embedded in C, or in 
an open disc, is automatically bipartite. Any bipartite quad-graph produces 
two dual polytopal (not necessarily quadrilateral) cell decompositions S and 
g*, with V{9) being the "black" vertices of D and V{S*) being the "white" 
ones, and edges of g (resp. of g*) connecting "black" (resp. "white") vertices 
along the diagonals of each face of D. 



In order to extract from this geometric construction some consequences 
for discrete Laplace type equations, we need the following deep and some- 
what mysterious property, which quite often accompanies the 3D consis- 
tency of quad-graph equations with fields on vertices ()4.4|) . We write here 




Figure 4.9: The star of the vertex 
xq in the graph g. 



Figure 4.10: Faces of D around 
the vertex xq. 
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this equation in slightly modified notations as 

Q{xo,yo,xi,yi;ao,ai) = 0. (4.32) 

For notational simplicity, vertices x stand here for the corresponding fields 
f{x); the edges {xo,yo), (xo,yi) carry the labels oq, ai, respectively. 

Definition 4.13 An equation (^3^ possesses a three-leg form centered 
at the vertex xq, if it is equivalent to the equation 

ip{xo,yo;ao) - tp{xo,yi; ai) = (j){xo,xi;ao,ai) (4.33) 

with some functions tp,(p. The terms on the left-hand side correspond to 
"short" legs {xQ,yQ), {xQ,yi) € E(T)), while the right-hand side corresponds 
to the "long" leg (xo,xi) G -E'(S). 

A summation of quad-graph equations for the flower of quadrilaterals 
adjacent to the "black" vertex xq € V{9) immediately leads, due to the 
telescoping effect, to the following statement. 

Theorem 4.14 a) Suppose that eq. i4.i^i^ on a bipartite quad-graph T) 
possesses a three-leg form. Then a restriction of any solution f : V{T)) — > C 
to the "black" vertices V{S) satisfies the discrete equations of the Laplace 
type, 

n 

^(j){xo,Xk;ak-i,ak) = 0, (4.34) 

k=l 

where n is the valence of the vertex xq in Q. 

b) Conversely, given a solution f : V{S) ^ C of these Laplace type equa- 
tions on a graph S coming from a cell decomposition of a simply- connected 
surface, there exists a one-parameter family of its extensions f : VifD) C 
satisfying eq. J^.^j ) on the double T). Such an extension is uniquely defined 
by a value at one arbitrary vertex from V{S*)- 

Sometimes it is more convenient to write the three-leg equation (|4.33jl in the 
multiplicative form: 

^'(xo, yo; ao)/^(a;o, yi; ai) = ^{xq, xi;ao, ai) (4.35) 

with some functions ^, <I>, so that the Laplace type equations (|4.34|) also 
become multiplicative: 

n 

Y\_^{xo,Xk;ak-i,ak) = I. (4.36) 
fc=i 
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This relation between integrable quad-graph equations and Laplace type 
equations was discovered in |BobSul] . The three- leg forms for all quad-graph 
equations of Theorem l4.8l were found in jAdBS2j . The next theorem provides 
three-leg forms for all equations of the lists Q and H (the results for the list 
A follow from these ones). The functions (j){x,y; a, f3), resp. ^{x,y; a, (3), 
corresponding to the "long" legs, yield integrable additive (resp. multiplica- 
tive) Laplace type equations on arbitrary planar graphs. 



Theorem 4.15 Three-leg forms exist for all equations from Theorem \^.^ 
(Ql)5=o- Additive three-leg form with 0(x, y; a, /?) = '0(x, y;a — 13), 

OL 

tl^{x,u;a) = . (4.37) 

X — u 

(Ql)<5=i.' Multiplicative three-leg form with ^{x, y; a, (3) = ^'(x, y; a — (3), 

^{x,u;a) = ^ + "~^ , (4.38) 

X — a — u 

(Q2); Multiplicative three-leg form with <I>(x, y; a, 13) = ^'(x, y;a — (3), 

^{x,u;a) = \- ^ , 4.39 

(A — aj^ — u 

where x = . 

(Q3)5=o-' Multiplicative three-leg form with ^{x,y; a, f3) = ^(x, y; a//3), 

^{x,u;a) = . 4.40 

X — au 

(Q3)5=i; Multiplicative three-leg form with ^{x,y;a, /3) = ^{x,y;a — P), 

^ , . sm(X -\- a) -u / . \ 

m{x,u;a) = , (4.41) 

sm(A — a) — u 

where x = sin(X). 

(Q4); Multiplicative three-leg form with <I>(x, y; a, (3) = ^'(x, y;a — (3), 

QiX + a) sn(X + a) — u , . . 

^{x,u-a) = -^ ' 4.42 

fc)(A — a) sn(A — a) — u 

where x = sn(X), and Q{X) is the Jacobi theta- function. 
(HI); Additive three-leg form with 

a- (3 

(l)(x,y;a, f3) = , ip(x,u;a) = x u. (4.43) 

x-y 
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(H2).- Multiplicative three- leg form with 



X — y + a — (3 
X — y — a + (3 



^(x, u;a) = X + u + a. 



(4.44) 



(H3)/ Multiplicative three- leg form with 



^{x,y;a,P) 



f3x — ay 



^'(x, u; a) = XU + 5a. 



(4.45) 



ax — f3y ^ 



One sees that there are only six "long" legs functions (j){x,y; a, (3), resp. 
^{x,y] a, (3), leading to integrable Laplace type equations on arbitrary pla- 
nar graphs. Three of them are rational in x,y. Each of the corresponding 
Laplace type equations admits two extensions to a 3D consistent quad-graph 
equation: one from the list Q, where the "short" legs essentially coincide with 
the "long" ones - (Ql)5=o, (Ql)5=i, and (Q3)5=o, and another one from the 
list H, with different "short" legs - (HI), (H2), and (H3). Other three func- 
tions $ are rational in y only, and require for a uniformizing change of the 
variable x. Corresponding Laplace type equations admit only one extension 
to 3D consistent quad-graph equations - (Q2), (Q3)5=i, and (Q4). Thus, 
restriction to a subgraph allows us to derive from 3D consistent quad-graph 
equations much more general systems, which clearly inherit the integr ability. 
See a detailed realization of this idea for the (Q4) equation in |AdSj . 

Remark. It should be mentioned that existence of a three-leg form 
allows us to derive (and, in some sense, to explain) the tetrahedron property 
of Sect. 14.2.11 Indeed, consider three faces adjacent to the vertex /123 on 
Fig. 1121 namely the quadrilaterals (/i, /12, /123, /13), (/2, /23, /123, /12), and 
(/s) /i3) /i23 5 /23)- A summation of the three- leg forms (centered at /123) of 
equations corresponding to these three faces leads to the equation 

0(/i23,/i;a2,a3) +'^(/i23,/2; 03,01) + 0(/i23 , /s ; «! , "2) = 0. (4.46) 

This equation relates the fields at the vertices of the "white" tetrahedron on 
Fig. 14.21 It can be interpreted as a discrete Laplace type equation coming 
from a spatial flower with three petals. 

4.4 Geometry of boundary value problems for 3D 
consistent equations 

We discuss here several aspects of the problem of embedding of a quad-graph 
into a regular multi-dimensional square lattice. As a combinatorial problem, 
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it was studied in a more general setting of arbitrary cubic complexes in 
IDoShU IDoSli2| IShShj . We are interested here in its relation to integrable 
equations. 

We start with the question about correct initial value problems for dis- 
crete 2D equations on quad-graphs. Let P be a path in the quad-graph D, 
i.e., a sequence of edges tj = {vj,Vj+i) £ E(T>). We denote by E{P) = {cj} 
and V{P) = {vj} the set of edges and the set of vertices of the path P, 
respectively. One says that the Cauchy problem for the path P is well- 
posed, if for any set of data fp : V{P) —i- C there exists a unique solution 
/ : V{V) ^ C such that f\v(P)= fp- 

The answer to this question was given in |AdVj with the help of the 
notion of a strip in D. 

Definition 4.16 A strip in T) is a sequence of quadrilateral faces qj € F[Tl) 
such that any pair qj-i, qj is adjacent along the edge tj = qj-i H qj, and tj, 
ej+i are opposite edges of qj . The edges tj are called traverse edges of the 
strip. 

So, any strip in D can be associated to a label a sitting on all its traverse 
edges Cj. 

Theorem 4.17 Consider one of the 3D consistent systems i4-4^ of Theo- 
rem \4.^ Let D be a finite simply connected quad-graph without self-crossing 
strips, and let P be a path without self- crossings in D. Then: 

i) If each strip in T> intersects P exactly once, then the Cauchy problem 
for P is well-posed. 

a) If some strip in T) intersects P more than once, then the Cauchy prob- 
lem for P is overdetermined (has in general no solutions). 

Hi) If some strip in D does not intersect P, then the Cauchy problem for 
P is underdetermined (has in general more than one solution). 

It should be mentioned that this theorem is not valid for equations without 
the 3D consistency property. One of the proofs of the statement i) in |AdVj 
is based on the embedding of T> into the unit cube of Z^, where N is the 
number of edges in P (the number of distinct strips in D). In this embedding 
the path P turns into the path (0,0,0, ... ,0), (1,0,0, . . . ,0), (1,1,0, . . . ,0), 
. . . , (1, 1, 1, . . . , 1). It is clear that for a 3D consistent equation the Cauchy 
problem for this path is well-posed. 

A different aspect of embeddability of the quad-graph D into a regu- 
lar multi-dimensional cubic lattice was studied in |BobMeS] . based on the 
following result |KeSj . 
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Theorem 4.18 A quad-graph T) admits an embedding in C with all rhombic 
faces if and only if the following two conditions are satisfied: 

i) No strip crosses itself or is periodic. 

ii) Two distinct strips cross each other at most once. 

Given a rhombic embedding p : V{Ti) — > C, one defines the fohowing func- 
tion on the directed edges of D with values in = {0 € C : \9\ = 1}: 



(9(x, y) = p{y) - p{x), V{x, y) G E{Ti). 



(4.47) 



This function can be called a labelling of directed edges, since it satisfies 
9{—t) = —0{t) for any e € E{D), and the values of 6 on two opposite 
and equally directed edges of any quadrilateral from F{T)) are equal. See. 
Fig. 14.111 For any labelling 9 : E{T)) §^ of directed edges, the function 
a = 9"^ : E(T)) — > S-*^ is a labelling of (undirected) edges of T> in our usual 
sense. 




Figure 4.11: Labelling of directed edges 



Definition 4.19 A rhombic embedding p : V{D) — s- C o/ a quad-graph D is 
called quasicrystallic, if the set of values of the function 9 : E{T)) 
defined by \4-4'l\j , is finite, say = {±0i, . . . , ^9d\. 

It is of a central importance that any quasicrystallic rhombic embedding p 
can be seen as a sort of a projection of a certain two-dimensional subcomplex 
(combinatorial surface) of a multi-dimensional regular square lattice Z*^. 
The vertices of JId are given by a map P : V{T>) TL^ constructed as 
follows. Fix some xq G Vi^^, and set P{xq) = 0. The images in Z'^ of all 
other vertices of D are defined recurrently by the property: 

for any two neighbors x,y € Yi^), if p(y) — p[x) = ±6*4 € G, then 
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P{y) — P{x) = zbcj , where Cj is the i-ih. coordinate vector of Z*^. 

Edges and faces of correspond to edges and faces of D, so that the 
combinatorics of is that of D. 

To exploit possibihties provided by the 3D consistency, extend the la- 
beHng 6 : E{T)) — s- §^ to all edges of l/, assuming that all edges parallel 
to (and directed as) Cj^ carry the label Oj.. This gives, of course, also the 
labelling a = 9^ oi undirected edges of l/. Now, any 3D consistent equation 
can be imposed not only on ri^, but on the whole of U^: 

Qif, fj, fjk, fk-, aj, = 0, l<j^k<d. (4.48) 

Here indices stand for the shifts into the coordinate directions. Obviously, 
for any solution / : Z"^ ^ C of eq. (ITIHl) . its restriction to ¥{^1,) ~ V{D) 
gives a solution of the corresponding equation on the quad-graph 2). As 
for the reverse procedure, i.e., for the extension of an arbitrary solution of 
eq. (|4.4|) from D to Z*^, more thorough considerations are necessary. An 
elementary step of such an extension consists of finding / at the eighth ver- 
tex of an elementary 3D cube from the known values at seven vertices, see 
Fig. 14.121 This can be alternatively viewed as a flip (elementary transfor- 
mation) on the set of rhombically embedded quad-graphs D, or on the set 
of the corresponding surfaces Qt> in Z*^. The 3D consistency assures that 
any quad-graph T) (or any corresponding surface Od) obtainable from the 
original one by such flips, carries a unique solution of eq. (|4.48j) which is an 
extension of the original one. 




Figure 4.12: Elementary flip 



Definition 4.20 For a given set V C Z'^, its hull "K^V) is the minimal set 
!K C Z*^ containing V and satisfying the condition: if three vertices of an 
elementary square belong to "K, then so does the fourth vertex. 
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One shows by induction that for an arbitrary connected subcomplex of Z'^ 
with the set of vertices V, its hull is a brick 

Ila,b = {n = {ni,...,nd) eZ'^ : ak<nk<bk, k = l,...,d}, (4.49) 

where 

ttk = ak{V) = minnk, = bk{V) = maxuk, k = l,...,d, (4.50) 

and in the case that iik are unbounded from below or from above on V, we 
set ak{V) = — oo, resp. bk{V) = oo. 

Combinatorially, all points of the hull 'K{V (Q-x))) can be reached by 
flips of Fig. I4.12( starting from fij). However, there might be obstructions, 
having nothing to do with 3D consistency, for extending solutions of eq. 
(|4.4|) from a combinatorial surface (two-dimensional subcomplex of Z'^) to 
its hull. For instance, the surface O shown on Fig. 14. 1^^ supports solutions of 
eq. (|4.4jl which cannot be extended to solutions of eq. (|4.48j) on the whole 
of !K(y(r2)): the recursive extension will lead to contradictions. The reason 
for this is a non-monotonicity of 17: it contains pairs of points which cannot 
be connected by a monotone path, i.e., by a path in Q with all directed 
edges lying in one octant of Z"^. However, such surfaces Q, do not come from 
rhombic embeddings, and in the case of Q^) there will be no contradictions. 

Theorem 4.21 Let the combinatorial surface Q-j) in come from a rhom- 
bic embedding of a quad-graph D, and let its hull be 3i(y(r2x))) = ^a,b- 
arbitrary solution of eq. j^.^D on Q.'X) can be uniquely extended to a solution 
of eq. U4^ on Iia,b- 




Figure 4.13: A non-monotone surface in 



Note that intersections of Vlxi with bricks correspond to combinatorially 
convex subsets of D, as defined in |Me2j . 



Chapter 5 



Discrete complex analysis. 
Linear theory 



5.1 Basic notions of the discrete linear complex 
analysis 

Many constructions in the discrete complex analysis are parallel to the dis- 
crete differential geometry in the space of the real dimension 2. 

Recall that a harmonic function n : ~ C ^ M is characterized by the 
relation: 

A conjugate harmonic function u : ~ C ^ M is defined by the Cauchy- 
Riemann equations: 

dv 
dy 

Equivalently, f = u + iv : 
Cauchy-Riemann equation 



The real and the imaginary parts of a holomorphic function are harmonic, 
and any real-valued harmonic function can be considered as a real part of a 
holomorphic function. 

A standard classical way to discretize these notions, going back to Fer- 
rand ^ and Duffin |Dulj . is the following. A function n : ^ M is called 
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du dv du 
dx ' dx dy 

~ C ^ C is holomorphic, i.e., satisfies the 
dy dx 
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discrete harmonic, if 

A natural definition domain of a conjugate discrete harmonic function v 
^ M is the dual lattice: 



1 




1 


> 



The defining discrete Cauchy-Riemann equations read: 



Vl 



Ul 



6 



VQ 



Ul 



vo o- 



Uo 



Vl — Vo = Ul — Uo 



vi-vo = -{ui - Uo) 



IS 



The corresponding discrete holomorphic function / : U (Z^)* —> 
defined on the superposition of the original square lattice and the dual 

one (Z^)*, by the formula / = ^ * ' which comes to replace the smooth 

one f = u + iv. Remarkably, the discrete Cauchy-Riemann equation for / 
is one and the same for the both pictures above: 
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fl 



/3 



fl 



/3 



fi- f2= i{f3 - fl)- 



This discretization of the Cauchy-Riemann equations apparently preserves 
the most number of important structural features, and the corresponding 
theory has been developed in [F | IDolj . A pioneering step in the direction 
of a further generalization of the notions of discrete harmonic and discrete 
holomorphic functions was undertaken by Duffin |Do2j . where the combina- 
torics of I? was given up in favor of arbitrary planar graphs with rhombic 
faces. A far reaching generalization of these ideas was given by Mercat 
|Melj . who extended the theory to discrete Riemann surfaces. 

Discrete harmonic functions can be defined for an arbitrary graph S with 
the set of vertices V{S) and the set of edges -E'(S)- 

Definition 5.1 For a given weight function v : i?(S) — on edges of S; 
the Laplacian is the operator acting on functions f : V{S) ^ C by 



where the summation is extended over the set of vertices x connected to xq 
by an edge. A function f : V{3) — > C is called discrete harmonic (with 
respect to the weights u), if Af = 0. 

The positivity of weights v in this definition is important from the analytic 
point of view, as it guarantees, e.g., the maximum principle for the discrete 
Laplacian under suitable boundary conditions (so that discrete harmonic 
functions come as minimizers of a convex functional). However, from the 
pure algebraic point of view, one might consider at times also arbitrary real 
(or even complex) weights. 

If S comes from a cellular decomposition of an oriented surface, let 9* 
be its dual graph, and let the quad-graph T> be its double, see Sect. 14., SL 
Extend the weight function to the edges of S* according to the rule 




(5.1) 




(5.2) 
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Definition 5.2 A function f : VCD) C is called discrete holomorphic 

(with respect to the weights v), if for any positively oriented quadrilateral 
{xo,yo,xi,yi) G F{D) (see Fia. \5.1\) there holds: 

f{yi) - f{yo) ■ / N 1 Q^ 
— — — = iu{xo, xi) = — -. (5.3) 

f{xi)-f{xo) i^{yo,yi) 

These equations are called discrete Cauchy-Riemann equations. 




Figure 5.1: Positively oriented quadrilateral, with a labelling of directed 
edges 

The relation between discrete harmonic and discrete holomorphic func- 
tions is the same as in the smooth case. It is given by the following statement 
which is a particular case of Theorem 14. 141 

Theorem 5.3 a) If a function f : V{T)) C is discrete holomorphic, then 
its restrictions to V{S) and to V{S*) are discrete harmonic. 

b) Conversely, any discrete harmonic function f : V{S) C admits a 
family of discrete holomorphic extensions to V{T>), differing by an additive 
constant on V{S*)- Such an extension is uniquely defined by a value at one 
arbitrary vertex y £ V{S*). 

5.2 Moutard transformation for Cauchy-Riemann 
equations on quad-graphs 

Observe that discrete Cauchy-Riemann equations (|5.3|) is formally not dif- 
ferent from the Moutard equations 1)2. 43() for T-nets. One only has to fix 
the orientation of all quadrilateral faces {xo,yo,xi,yi) G F{D). We assume 
that it is inherited from the orientation of the underlying surface. 

One can apply now the Moutard transformation of Sect. 12.31 to discrete 
holomorphic functions. To this aim, one has to choose the orientation of 
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all elementary quadrilaterals on Fig. 15.21 This can be done, for example, as 
follows: for the quadrilaterals {xq ,xf ,yf ) G choose the orienta- 

tion to coincide with that of the corresponding (xq, yo,xi,yi) £ F(T)). For a 
"vertical" quadrilateral over an edge {x,y) G E{D), assume that x E Vi9), 
y £ V^(S*), and choose the positive orientation corresponding to the cyclic 
order {x,y,y~^ ,x^) of its vertices. Observe that under this convention, two 
opposite "vertical" quadrilaterals are always oriented differently. 



In the case of arbitrary quad-graphs, one has to generalize one more 
ingredient of the Moutard transformation, namely the data (MT^). 

Theorem 5.4 On an arbitrary bipartite quad-graph T), valid initial condi- 
tions for a Moutard transformation of the discrete Cauchy-Riemann equa- 
tions consist of 

(MCRf ) value of f at one point x+(0) G V{'D+); 

(MCR^) values of weights on "vertical" quadrilaterals {x,y,y^,x^) at- 
tached to all edges (x,y) of a Cauchy path in D. 

See Theorem 14.171 for necessary and sufficient conditions for a path to be a 
Cauchy path, i.e., to support initial data for a well-posed Cauchy problem. 
It is natural to assign the weights on the "vertical" quadrilaterals to the 
underlying edges of D. 

Weights u on the faces of V together with the data (MCR^) yield the 
transformed weights i^^ on the faces of 25"*", as well as the weights over all 
edges of E{T)). This can be considered as a Moutard transformation for 





Figure 5.2: Elementary cube of D 
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Cauchy-Riemann equations on D. Finding a solution f : V(T)^) ^ C of the 
transformed equations requires additionally the datum (MCR^). 

Note that the system of weights is highly redundant, due to eq. (|5.2|1 . 
To fix the ideas in writing the equations, we stick to the weights attached to 
the "black" diagonals of the quadrilateral faces of the complex D. On the 
ground floor, these are the edges of the "black" graph S; on the first floor, 
these are the edges of the "black" graph which is a copy of S*; and for the 
"vertical" faces, these are the edges with x G ^(S) and y E V{S*). 

Needless to say that the latter weights can be assigned to the quad-graph 
edges G E{D). So, we write the discrete Cauchy-Riemann equations 

as follows: 

/(yi)-/(yo) = iz^(xo,xi)(/(xi) - /(xo)), (5.4) 

f{x+)-f{xt) = Hy^,yWiyt)-f{y^)), (5-5) 

f{x+)-fiy) = iv{x,y+){f{y+)-f{x)). (5.6) 

Denote, for the sake of brevity, 

v = u{xQ,xi), = u{y^,yf), njk = i^{xj,y'^). 

Regarding the weights z/, ^qq, and //qi as the input of the Moutard transfor- 
mation on an elementary hexahedron of D, its output consists of the weights 
W^, ^10, and /Uii, given by (cf. eq. ()2.46|) ') 

u^v = -/Uii/ioo = -Worn = • (5.7) 

/^oo — A'oi — ^ 

This transformation is well defined for real weights fijk, but it does not 
preserve, in general, positivity of the weights i^. 

To give a different form of this transformation, observe that relation 
/^ii/^oo = A^ioA^oi for each elementary quadrilateral {xo,yo, xi,yi) of D yields 
the existence of the function 9 : V{D) — > C, defined up to a constant factor, 
such that injk = 9{yk)/9{xj). Moreover, choosing 9{xq) real at some point 
xq E ^(9), one sees that 9 takes real values on V{^) and imaginary values 
on y(S*)- An easy computation shows that the last equation in (|5.7|) is 
equivalent to 

6{yi) - ^{yo) = ii^{xo,xi){9{xi) - 9{xo)), 

so that the function 9 is discrete holomorphic with respect to the weights v. 
For the transformed weights z^"^ one finds: 
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Conversely, an arbitrary discrete holomorphic function : V{D) — > C 
defines, via eq. (|5.8|) . a Moutard transformation of the discrete Cauchy- 
Riemann equations. It should be mentioned that the data (MCR2^) can be 
re- formulated in terms of the function 9: 

(MCR^) values of 9 at all vertices along a Cauchy path in 2). 

Remark. Moutard transformation for discrete Cauchy-Riemann equations 
yields, by restriction to the "black" graphs, a sort of Darboux transformation 
of arbitrary discrete Laplacians on 9 into discrete Laplacians on S*. A 
further discussion of this transformation for discrete Laplace operators can 
be found in [NieSDl ID^TTN?^ . 



5.3 Integrable discrete Cauchy-Riemann equations 

We now turn to a fruitful question about "stationary points" of the Moutard 
transformation discussed in the previous section. More precisely, this is 
the question about conditions on the weights : -E'(S) — such that 
there exists a Moutard transformation for which the opposite faces of any 
elementary hexahedron of D (see Fig. 15. 2|) carry identical equations. 

Theorem 5.5 A system of discrete Cauchy-Riemann equations with the 
function v : £'(3) U -E(3*) satisfying i5. gj) admits a Moutard trans- 

formation into itself, if and only if for all xq £ V{S) and for all yo € 1^(3*) 
the following conditions are fulfilled: 

eestar(xo;g) e*Gstar(j/o;g*) 

Proof. Opposite faces of D and carry identical equations, if z^'^z^ = 1 
in eq. 1)5. 7|) . Clearly, this yields also /Un/ioo = /"loA^oi = —1) which means 
that the opposite "vertical" faces also support identical equations (recall 
that opposite "vertical" faces carry different orientation). Moreover, given 
v = i'{xq,Xi) for an elementary quadrilateral {xQ,yQ,xi,yi) of D, we find 
that the input data /xqo, fJ-oi of the Moutard transformation should be related 
as follows: 

i^mm 1 ^ Moo - /^l -A r 1 

- 1 ^ noi = — - = ^00 , 

UnnV +1 \l' I 



Moo - Moi - Mooi^ + 1 



where the standard notation for the action of PGL2{C) on C by Mobius 
transformations is used. This means that all the weights on the vertical 
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faces of a "stationary" Moutard transformation are completely defined by 
just one of them, so that such transformations form a one-parameter family. 
To derive a condition for v for the existence of a "stationary" Moutard 
transformation, consider a flower of quadrilaterals {xo,yk-i, Xk,yk) around 
xq € V{9) (see Fig. 14.10')) . In the natural notations, we find: 

^o,fc-i - ( I -yk\ r 1 

Running around xq should for any ^uoo return its value, which means that 
the matrix product = ^ ^ j^'^^ should be proportional to 

the identity matrix. This matrix product is easily computed: 

^ = ^ ( 11(1 + iuk) + - i^k)) , ^ = ^ ( 11(1 + ii.,) - n(i - i^k)) , 

and the condition i? = is equivalent to the first equality in 1)5. 9() . The sec- 
ond condition in 1)5. 9() is proved similarly, by considering a flower of quadri- 
laterals around yo ^ V{S*)- □ 

Thus, existence of a "stationary" Moutard transformation singles out 
a special class of discrete Cauchy-Riemann equations, which have to be 
considered as 2D systems with the 3D consistency property, see Sect. 14.11 In 
other words, such Cauchy-Riemann equations should be termed integrable. 
The main difference as compared with the examples in Sect. 14.11 is that 
discrete Cauchy-Riemann equations naturally depend on the orientation of 
the elementary quadrilaterals, and that their parameters v are apparently 
assigned not to the edges of the quad-graph, but rather to diagonals of its 
faces. 

The integrability condition ()5.9)) admits a nice geometric interpretation. 
It is convenient (especially for positive real-valued z/) to use the notation 

i/(e)=tan^, 0(e)€(O,7r). (5.10) 

The condition i^(e*) = l/i'(e) is translated into 

^{e*) = 7T-<P{e), (5.11) 

while the condition (|5.9j) says that for all xq € V{3) and for all yo G V{S*) 
there holds: 

n exp(#(e)) = l, n exp(i<^(e*)) = 1. (5.12) 

eestar{xo;g) e*estar{j/o;S*) 
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These conditions should be compared with conditions |KeSj characterizing 
the angles cp ■ E{Q)\JE{Q*) — > (0, vr) of a rhombic embedding of a quad-graph 
D, which consist of H5.11() and 

^ <A(e) = 27r, Yl '^(e*) = 2vr, (5.13) 

e6star(a;o;S) e*6star(j/o;S*) 

for all xq G V{9) and for all yo £ ^(9*)- Thus, the integrability condition 
()5.12|) says that the system of angles (j) : E(3) U E(Q*) — > (0, vr) comes from 
a realization of the quad- graph 2) as a rhombic ramified embedding in C. 
Flowers of such an embedding can wind around its vertices more than once. 

Another formulation of the integrability conditions is given in terms of 
the edges of the rhombic realizations. 

Theorem 5.6 Integrability condition i5.y\) for the weight function : E{Q)\J 
E{9*) M_|_ is equivalent to the following one: there exists a labeling 
6 : E^D) E>^ of directed edges of D such that, in notations of Fig. \5.1\ 

I \ ^ .00-01 iA\ 

u{xo,xi) = — -=t——-. 5.14 

i^{yo,yi) 00 + 01 

Under this condition, the 3D consistency of the discrete Cauchy-Riemann 
equations is assured by the following values of the weights v on the diagonals 
of the vertical faces of D ; 

i/(x,y+) = (5.15) 

where = 0{x,y), and X C is an arbitrary number having the interpreta- 
tion of the label carried by all vertical edges ofD: A = 9{x,x^) = 9{y^y^). 

So, integrable discrete Cauchy-Riemann equations can be given a form with 
parameters attached to directed edges of D: 

f{yi) - f{yo) _ 0i-0o . . 

f{x,)-f{xo) 0i+0o' ^ ' ' 

where 

00 = p{yo) - p{xo) = p{xi) - p{yi), 01 = p{yi) - p{xq) = p{xi) - p{yo), 
and p:y(S)^Cisa rhombic realization of the quad-graph 2). Since 

01 - 00 ^ pjyi) -pjyo) 

01 + 00 Pixi)-p{xo)' 
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we see that for a discrete holomorphic function / : V{2) C, the quotient 
of diagonals of the /-image of any quadrilateral {xo,yo, xi,yi) £ F{T)) is 
equal to the quotient of diagonals of the corresponding rhombus. 

A standard construction of zero curvature representation for 3D con- 
sistent equations, given in Theorem 14.31 leads in the present case to the 
following result. 

Theorem 5.7 The discrete Cauchy-Riemann equations iS.lb]) admit a zero 
curvature representation with spectral parameter dependent 2x2 matrices 
along (x, y) G E{D) given by 

T( ^^ (^^^ -2e{f{x) + f{y))\ 

L{y,x,a;\) = , (5.17) 

y x-e J 

where 9 = p{y) — p{x). 

Linearity of the discrete Cauchy-Riemann equations is reflected in the tri- 
angular structure of the transition matrices. 

Also, all constructions of Sect. 14.41 can be applied to integrable dis- 
crete Cauchy-Riemann equations. In particular, for weights coming from 
a quasicrystallic rhombic embedding of the quad-graph D, with the labels 
Q = {±^1, . . . , ±9d}, discrete holomorphic functions can be extended from 
the corresponding surface Oxi C Z*^ to its hull, preserving discrete holomor- 
phy. Here we have in mind the following natural definition: 

Definition 5.8 A function / : Z'^ — > C is called discrete holomorphic, if 

it satisfies, on each elementary square of JJ^ , the equation 

f{n + ej + ek)- f{n) _ 6^ + 9^ 

f{n + ej)-f{n + ek) 9^-9^' ^ ' ' 

For discrete holomorphic functions in Z'^, the transition matrices along the 
edges (n, n -|- e^,) of Z'^ are given by 

(x + 9, -29k{f{n + ek) + fin))\ 
Lk{n;X) = . (5.19) 

\ X-9k J 

All observations of this section hold also in the case of generic complex 
weights z^, which leads to G C and to parallelogram realizations of T). 
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5.4 Discrete exponential functions 

A discrete exponential function on quad-graphs T) was defined and studied 
in |Mel[ iKej . Its definition, for a given rhombic embedding p : V{T>) C, 
is as follows: fix a point xq € V{T>). For any other point x £ V{T>), choose 
some path {cj}"^]^ C -E(D) connecting xq to x, so that Zj = {xj-i,Xj) and 
Xn = X. Let the slope of the jth edge be 9j = p{xj) — p{xj-i) € S"*^. Then 



e(x;z) = Yl 



Z + i 



j=i -I 

Clearly, this definition depends on the choice of the point xq € V^(I'), but 
not on the path connecting xq to x. A question posed in |Kej is whether dis- 
crete exponential functions form a basis in the space of discrete holomorphic 
functions on 2). 

An extension of the discrete exponential function from VLd to the whole 
of Z*^ is the discrete exponential function, given by the following simple 
formula: 

^("-) - n (^)"'- (^-^o) 

k=l 

For d = 2, this function was considered in [F | IDolj . The discrete Cauchy- 
Riemann equations for the discrete exponential function are easily checked: 
they are equivalent to a simple identity 

fz + 9j z + Ok ^ / (z + ^j z + ek\ _9j+ Ok 



\z — 9j z — 9k J I \z — 9j z — 9^) 9j — 9k 

At a given n E Z'^, the discrete exponential function is rational with respect 
to the parameter z, with poles at the points €i9i, . . . , ed9d, where = 
signnfc. 

Equivalently, one can identify the discrete exponential function by its 
initial values on the axes: 

einek;z) = (^Y. (5.21) 



A still another characterization says that e(-; z) is the Backlund transforma- 
tion of the zero solution of discrete Cauchy-Riemann equations on Z*^, with 
the "vertical" parameter z. 

We now show that the discrete exponential functions form a basis in 
some natural class of functions (growing not faster than exponentially), thus 
answering in affirmative the above mentioned question by Kenyon. 
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Theorem 5.9 Let f be a discrete holomorphic function on V{T)) ~ V{0,X))> 
satisfying 

|/(n)| <exp(C([ni| + ... + |nrf|)), Vn € ^(l^i,), (5.22) 

with some C G M. Extend it to a discrete holomorphic function on "KiV {O,-]^)) . 
There exists a function g defined on the disjoint union of small neighbor- 
hoods around the points ±6^ € C and holomorphic on each one of these 
neighborhoods, such that 

fin)- f{0) = ^. j^g{X)e{n-\)d\ ^ne%{V{^^)), (5.23) 

where T is a collection of 2d small loops, each one running counterclockwise 
around one of the points ±^fc. 

Proof is constructive and consists of three steps. 

(i) Extend / from V{Qx>) to "K^V^Qx)))', inequaUty 1)5. 22(1 propagates in 
the extension process, if the constant C is chosen large enough. 

(ii) Introduce the restrictions fi^^ of / : 'K{V{^'X))) ^ C to the coordinate 
axes: 

Z^'^) =/(nefc), ak{n^)<n<bk{nv)- 

(iii) Set g{X) = X]fc=i(fi'fc(-^) + S'-fc('^))i where the functions 5±fc(A) vanish 
everywhere except in smaU neighborhoods of the points ±0^, respec- 
tively, and are given there by convergent series 

»(A) = ^ (^/!" - /(O) + t (^)"(/.Si - . (5.24) 

and a similar formula for (7_fc(A). Formula (|5.23)) is then easily verified 
by computing the residues at A = ±6k- 

(k) 

It is important to observe that the data /A , necessary for the construction 
of ^(A), are not among the values of / on V(T)) ~ V{flxi) known initially, 
but are encoded in the extension process. 
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5.5 Discrete logarithmic function 

The discrete logarithmic function on T> can be defined as fohows |Kej . Fix 
some point xq £ V{T>), and set 



Here the integration path T is the same as in Theorem 15.91 and fixing xq is 
necessary for the definition of the discrete exponential function on D. To 
make ()5.25|) a valid definition, one has to specify a branch of log(A) in a 
neighborhood of each point ±0^- This choice depends on x, and is done as 
follows. 

Assume, without loss of generality, that the circular order of the points 
ib^fc on the positively oriented unit circle is the following: 6i, . . ., 0d, 
—9i, . . ., —9ii. We set Oi^^d = —9^ for = 1, . . . , d, and then define Or for 
all r £ Ij by 2d-periodicity. For each r G Z, assign to 9r = exp{ijr) G S"^ a 
certain value of argument 7r G M: choose a value 71 of the argument of 61 
arbitrarily, and then extend it according to the rule 



Clearly, there holds jr+d = 7r + tt, and therefore also 7r+2d = 7r + 27r. It 
will be convenient to consider the points Or, supplied with the arguments 
'jr, as belonging to the Riemann surface A of the logarithmic function (a 
branched covering of the complex A-plane) . 

For each m G Z, define the "sector" Um on the embedding plane C of the 
quad-graph T) as the set of all points of V(T)) which can be reached from 
xo along paths with all edges from {9m, ■ ■ ■ , ^m+d-i}- Two sectors Um^ and 
have a non-empty intersection, if and only if \mi — < d. The union 
U = Um=-oo ^rn is a branched covering of the quad-graph D, and serves as 
the definition domain of the discrete logarithmic function. 

The definition (|5.25|) of the latter should be read as follows: for x G 
Um, the poles of e(x; A) are exactly the points 9m, ■ ■ ■ , 9m+d—i G A. The 
integration path F consists of d small loops on A around these points, and 
arg(A) = Silog(A) takes values in a small open neighborhood (in M) of the 
interval 





(5.25) 



7r.+i - 7r. G (0, vr) 



Vr G Z. 



[7m,7m+d-l] 



(5.26) 



of length less than vr. If m increases by 2d, the interval (|5.26|) is shifted 
by 2tt. As a consequence, the function £ is discrete holomorphic, and its 
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restriction to the "black" points V{S) is discrete harmonic everywhere on U 
except at the point xq- 

M{x) = (5.27) 

Thus, the functions in the integral representation H5.23|) of an arbi- 
trary discrete holomorphic function, defined originahy in disjoint neighbor- 
hoods of the points Or, in the case of the discrete logarithmic function are 
actually restrictions of a single analytic function log(A)/(2A) to these neigh- 
borhoods. This allows one to deform the integration path T into a connected 
contour lying on a single leaf of the Riemann surface of the logarithm, and 
then to use standard methods of the complex analysis in order to obtain 
asymptotic expressions for the discrete logarithmic function. In particular, 
one can show |Kej that at the "black" points V{S) there holds: 

£(x) ~ log |a; — xqI, X — > oo. (5.28) 

Properties (|5.27|1 . (|5.28|) characterize the discrete Green's function on S- 
Thus: 

Theorem 5.10 The discrete logarithmic function on T), restricted to the 
set of vertices V{9) of the "black" graph S, coincides with discrete Green's 
function on S. 

Now we extend the discrete logarithmic function to Z"^, which will allow 
us to gain significant additional information about it |BobMeS] . Introduce, 
in addition to the unit vectors S Z'^ (corresponding to 6^ G §^), their 
opposites e^+d = —^k, k G [l,d] (corresponding to 9k+d = —dk)-, and then 
define for all r G Z by 2(i-periodicity. Then 

m+d—l 

5m = Zer C Z'^ (5.29) 

is a d-dimensional octant containing exactly the part of which is the 
P-image of the sector Um C T). Clearly, only 2d different octants appear 
among Sm (out of 2^^ possible d-dimensional octants). Define Sm as the 
octant Sm equipped with the interval (|5.26j) of values for 9 log (6*^.). By 
definition, Smi and Sm2 intersect, if the underlying octants Smi and Sm2 
have a non-empty intersection spanned by the common coordinate semi- 
axes Zcr, and 91og(0r') for these common semi-axes match. It is easy to 
see that Sm-^ and Sm2 intersect, if and only if |mi — < d. The union 
^ = Um=-oo is a branched covering of the set Um=i ^ 
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Definition 5.11 The discrete logarithmic function on S is given by the 
formula 

i{n) = ^ [ i^^l^ e(n; X)dX, Vn G 5, (5.30) 

where for n € Sm the integration path T consists of d loops around 6m, ■ ■ 
6m+d~i on A, and 91og(A) on T is chosen in a small open neighborhood of 
the interval h5.2b}) . 

The discrete logarithmic function on D can be described as the restriction 
of the discrete logarithmic function on 5" to a branched covering of il^) ~ ^• 
This holds for an arbitrary quasicrystallic quad-graph 2) whose set of edge 
slopes coincides with = {±0i, . . . , ±0^}. 

Now we are in a position to give an alternative definition of the discrete 
logarithmic function. Clearly, it is completely characterized by its values 
i{ner), r G [m,m + d—l], on the coordinate semi-axes of an arbitrary octant 
Sm- Let us stress once more that the points ne^ do not lie, in general, on 
the original quad-surface Qd- 

Theorem 5.12 The values £n^ =£{ner), r € [m,m + d—l\, of the discrete 
logarithmic function on Sm C S are given by: 

£(:) = ( + i + + " (5.31) 

log{9r) = i^r, n odd. 

Here the values \og{9r) = i^r <ire chosen in the interval i5.26]) . 

Proof. Comparing formula H5.3U|) with H5.24() . we see that the values in '^ 
can be obtained from the expansion of log(A) in a neighbourhood of X = Or 
into the power series with respect to the powers of (A — 9r)/{X + Or). This 
expansion reads: 

1 /^^ ^ ^ 1 " (-1)" / A - 0, \ " 

log(A)=log(^.) + j:^-^(^) . 

n=l 

Thus, we come to a simple difference equation 

n{ii:l-i'^l,) = l-{-ir, (5.32) 
with the initial conditions 

= e{0) = 0, = ^(er) = log(0.), (5.33) 
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which yields eq. H5.31() . □ 

Observe that values (|5.31() at even (resp. odd) points imitate the be- 
haviour of the real (resp. imaginary) part of the function log(A) along the 
semi-lines arg(A) = arg{6r). This can be easily extended to the whole of S. 
Restricted to black points n £ S (those with ni + . . . + nrf even), the discrete 
logarithmic function models the real part of the logarithm. In particular, it 
is real-valued and does not branch: its values on Sm depend on m (mod 2d) 
only. In other words, it is a well defined function on Sm- On the contrary, 
the discrete logarithmic function restricted to white points n G S (those 
with ni + . . . + Ud odd) takes purely imaginary values, and increases by 27ri, 
as m increases by 2d. Hence, this restricted function models the imaginary 
part of the logarithm. 

It turns out that recurrent relations (|5.32|) are characteristic for an im- 
portant class of solutions of the discrete Cauchy-Riemann equations, namely 
for the isomonodromic ones. In order to introduce this class, recall that 
discrete holomorphic functions in Z'^ possess a zero curvature representa- 
tion with the transition matrices 1)5.19(1 . The moving frame ^'(•, A) : Z'^ — > 
GL2{C)[X] is defined by prescribing some ^'(0;A), and by extending it re- 
currently according to the formula 

^{n + ek;X)=Lk{n;X)^{n;X). (5.34) 

Finally, define the matrices A{-; A) : Z'^ ^ GL2(C)[A] by 

A{n;X) = ^^^^-Hn;X). (5.35) 

These matrices satisfy a recurrent relation, which results by differentiating 
(OH) . 

A{n + efe; A) = ^^^^^^^ '(n; A) + Lfc(n; X)A{n; X)L^\n; A), (5.36) 

and therefore they are defined uniquely after fixing some A{0; A). 

Definition 5.13 A discrete holomorphic function / : Z"' ^ C is called 
isomonodromic, if, for some choice of A{0;X), the matrices A(n; X) are 
meromorphic in X, with poles whose positions and orders do not depend on 
n G Z"^. 

This term originates in the theory of integrable nonlinear differential equa- 
tions, where it is used for solutions with a similar analytic characterization 
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It is clear how to extend Definition 15. 131 to functions on the covering S. 
In the following statement, we restrict ourselves to the octant 5*1 = 
for notational simplicity. 

Theorem 5.14 The discrete logarithmic function is isomonodromic: for a 
proper choice o/^(0; A), the matrices A{n; A) at any point n € have 
simple poles only: 



=1 



with 



A^'Hn) = ^ ^ , (5.38) 

,0 



5^'^(n) = ni\ , (5.39) 







any point n € 5", t/ie following constraint holds: 
d 

ni {i(ri + ei) - i(ri - e^)) = 1 - (_i)«i+-+"d. (5.4I) 

i=\ 

Proof. The proper choice of A(0; A) mentioned in the Theorem, can be read 
off formula ()5.38|) : ^(0; ~ I ■ proof consists of two parts. 

(i) First, one proves the claim for the points of the coordinate semi-axes. 
For any r = 1, . . . , d, construct the matrices A{ner; A) along the r-th 
coordinate semi-axis via formula (|5.36j) with the transition matrices 
1)5. 19(1 . This formula shows that singularities of A{ner; X) are poles 
at A = and at A = ±0r, and that the pole A = remains simple 
for all n > 0. By a direct computation and using mathematical in- 
duction, one shows that it is exactly the recurrent relation (|5.32|) for 
fn^ = /(ne^) that assures the poles A = ±0^ to remain simple for all 
n > 0. Thus, eq. (|5.37j) holds on the r-th coordinate semi-axis, with 
BW(ne,.) = CW(ne,) = for Z 7^ r. 



130 



CHAPTER 5. DISCRETE LINEAR COMPLEX ANALYSIS 



(ii) The second part of the proof is conceptual, and is based upon the mul- 
tidimensional consistency only. Proceed by induction, whose scheme 
follows filling out the hull of the coordinate semi-axes: each new point 
is of the form n + ■ + e^, j ^ k, with three points n, n + ej and 
n + ejt known from the previous steps, where the statements of the 
proposition are assumed to hold. Suppose that (|5.37|) holds at n + ej, 
n+efc. The new matrix A^n+ej + ef;.; A) is obtained by two alternative 
formulas, 

A{n + Ej + Cfc; A) = -jj-^ (n + ej; A) 

+Lfc(n + ej;X)A{n + ej; X)L^^{n + e^; A), (5.42) 

and the one with interchanged roles of k and j. Eq. (|5.42|) shows that 
all poles of A{n+ej + ek; A) remain simple, with the possible exception 
of A = ±0k, whose orders might increase by 1. The same statement 
holds with k replaced by j. Therefore, all poles remain simple, and 
(|5.37|) holds at n + ej + e^. Formulas (|5.38|) - (|5.4flj) and constraint 
(|5.4H) follow by direct computations based on eq. (|5.42j) . □ 



Chapter 6 



Discrete complex analysis. 
Integrable circle patterns 



6.1 Circle patterns 

The idea that circle packings and, more generally, circle patterns serve as 
a discrete counterpart of analytic functions is by now well established, see 
a popular review |Stlj and a monograph |St2j . The origin of this idea is 
connected with Thurston's approach to the Riemann mapping theorem via 
circle packings, see |Thlj . Since then the theory bifurcated to several areas. 

One of them is mainly dealing with approximation problems, and in this 
context it is advantageous to stick from the beginning with some fixed regu- 
lar combinatorics. The most popular are hexagonal packings, for which the 
C°° convergence to the Riemann mapping has been established |RSul IHeSj . 
Similar results are available also for circle patterns with the combinatorics 
of the square grid introduced by Schramm |Schrj . 

Another area concentrates around the uniformization theorem of Koebe- 
Andreev-Thurston, and is dealing with circle packing realizations of cell 
complexes of prescribed combinatorics, rigidity properties, constructing hy- 
perbolic 3-manifolds, etc. [Th2l iMRl iBe^ ITfej . A new variational ap- 
proach to this area is given in |Bob Sp' , where also a review of related re- 
sults can be found. An important differential-geometric application of this 
approach is the construction of discrete minimal surfaces through circle pat- 



terns BobHoSp 



The Schramm circle patterns with the combinatorics of the square grid 
constitute also a prominent example in the area we are mainly interested 
here, which deals with the interrelations of circle patterns with integrable 
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systems. We give here a presentation of several results in this area, based 
on |BohP3l iHobH^ \AgBl\ . 



Definition 6.1 Let S be an arbitrary cell decomposition of an open or closed 
disk inC A map z :V{3) C defines a circle pattern with the combi- 
natorics ofS, if the following condition is satisfied. Let y € i^(9) ~ ^(S*) 



be an arbitrary face of S, and let xi,X2, 



, x„ be its consecutive vertices. 



Then the points z{xi), z{x2), . . . , z{xn) G C lie on a circle, and their circular 
order is just the listed one. We denote this circle by C{y), thus putting it 
into a correspondence with the face y, or, equivalently, with the respective 
vertex of the dual cell decomposition S*. 

As a consequence of this condition, if two faces yo,yi G -^(S) have a com- 
mon edge (xo,xi), then the circles C(yo) and C{yi) intersect in the points 
z{xi),z{x2). In other words, the edges from i?(S) correspond to pairs of 
neighboring (intersecting) circles of the pattern. Similarly, if several faces 
yi,y2, ■ ■ ■ ,ym £ FiS) meet in one point xq G ^(S), then the corresponding 
circles C(yi), C{y2), • • • , C{ym) also have a common intersection point z{xq). 




Figure 6.1: Circle pattern 
Given a circle pattern with the combinatorics of S, we can extend the 
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function z to the vertices of the dual graph, setting 

z{y) = center of the circle C{y), y G F{Q) ~ V{3*)- 

After this extension, the map z is defined on all of VCD) = V{S) U V{3*), 
where D is the double of S. Consider a face of the double. Its z-image is a 
quadrilateral of the kite form, whose vertices correspond to the intersection 
points and the centers of two neighboring circles Cq, Ci of the pattern. De- 
note the radii of Co,Ci by ro,ri, respectively. Let xq,xi correspond to the 
intersection points, and let yo,yi correspond to the centers of the circles. 
Give the circles Co,Ci a positive orientation (induced by the orientation of 
the underlying C), and let cp ^ (0, tt) stand for the intersection angle of these 
oriented circles. This angle (j) is equal to the kite angles at the "black" ver- 
tices z{xo), z{xi), see Fig. 15.21 where the complementary angle cp* = ir — 4> 
is also shown. It will be convenient to assign the intersection angle (p = 0(e) 
to the "black" edge e = {xo,xi) G E{Q), and to assign the complementary 
angle cf)* = (p^e*) to the dual "white" edge e* = {yo,yi) G -^(3*). Thus, the 
function (p : E{3) U E{S*) ^ (0,7r) satisfies eq. (lOTl) . 




Figure 6.2: Two intersecting circles 



The geometry of Fig. IHT^ vields following relations. First of all, the cross- 
ratio of the four points corresponding to the vertices of a quadrilateral face 
of D is expressed through the intersection angle of the circles Cq, Ci: 

q{z{xo), z{yo), z{xi), z{yi)) = exp(2i0*). (6.1) 

Further, running around a "black" vertex of D (a common intersection point 
of several circles of the pattern) , we see that the sum of the consecutive kite 
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angles vanishes (mod 2tt), hence: 

Yl exp(i(/>(e)) = 1, VxoGF(g). (6.2) 

eGstar(a;o;S) 

Finally, let V'oi be the angle of the kite {z{xq), z{yo), z{xi), z{yi)) at the 
"white" vertex z{yo), i.e., the angle between the half-lines from the center 
z{yo) of the circle Co to the intersection points z{xo), z{xi) with its circle 
Ci. It is not difficult to calculate this angle: 

exp zV'oi = ] r—JJ\- 6.3 

Running around the "white" vertex of D, we come to the relation 
ro + r,- exp(i0*) 

n ^^^;,p(L4 -'. v.o.i-e-), (0.4) 

where the product is extended over all edges = {yo,yj) € star(yo; 9*), and 
= (p{ej), while rj are the radii of the circles Cj = C{yj). This formula 
has been used in |BobSp| as the basis of a variational proof of the existence 
and uniqueness of Delaunay circle patterns with prescribed combinatorics 
and intersection angles. 



6.2 Integrable cross-ratio and Hirota systems 

Our main interest is in the circle patterns with prescribed combinatorics and 
with prescribed intersection angles for all pairs of neighboring angles. (This 
is not the only class of patterns deserving a study from the point of view of 
integrability, see, e.g., |BobHo!jul for a different integrable class.) According 
to eq. (|(i.lj) . prescribing all intersection angles amounts to prescribing cross- 
ratios for all quadrilateral faces of the quad-graph T>. Thus, we come to the 
study of cross-ratio equations on arbitrary quad-graphs. 

Let there be given a function Q : £^(9) U -£-(9*) C satisfying the 
condition 

Q{e*) = l/Q{e), Ve G E{S). (6.5) 



Definition 6.2 The cross-ratio system on D corresponding to the func- 
tion Q, consists of the following equations for a function z : V{T)) C, one 
for any quadrilateral face {xQ,yo,xi,yi) ofT>: 



q{z{xo), z{yo), z{xi),z{yi)) = Q{xo,xi) = l/Q{yo,yi). (6.6) 
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An important distinction from the discrete Cauchy-Riemann equations is 
that the cross-ratio equations actually do not feel the orientation of quadri- 
laterals. 

We have already encountered 3D consistent cross-ratio systems on Z"^ 
in Sect. 14.11 (see eq. (|4.ip '). in the version with labelled edges. A natural 
generalization to the case of arbitrary quad-graphs is this: 




Figure 6.3: Quadrilateral, with a labelling of undirected edges 



Definition 6.3 A cross-ratio system is called integrable, if there exists a 
labelling a : E{D) C of undirected edges of D such that the function Q 
admits the following factorization (in notations of Fia. [U~S\) : 

Q{xo,xi) = -—^ = — . (6.7) 

Q{yo,yi) ai 

Clearly, integrable cross-ratio systems are 3D consistent (see Theorem 12. 28() . 
admit Bdcklund transformations, and possess zero curvature representation 
with the transition matrices H4.14() . It is not difficult to give an equivalent 
re- formulation of the integrability condition (|6.7() . 

Theorem 6.4 A cross-ratio system with the function Q : E{S)\-iE{S*) C 
is integrable, if and only if for all xq G V^(9) and for all E V{^*) the 
following conditions are fulfilled: 

n Q(e) = l, n Q{e*) = l. (6.8) 

eGstar(xo;S) e*Gstar(j/o;S*) 

For a labelling of undirected edges a : E{T)) — > C, there can be found a 
labelling 9 : E(T>) ^ C of directed edges such that a = 9"^. The function 
p : V{T)) — > C defined by p{y) — p{x) = 9{x,y) gives, according to eq. (|6.8|) . 
a parallelogram realization (ramified embedding) of the quad-graph 2). The 
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cross-ratio equations are written as 

q{z{xo),z{yo),z{xi),z{yi)) = -| = q{p{xo),p{yo),p{xi),p{yi)), (6.9) 

^1 

in other words, the cross-ratio of the vertices of the /-image of any quadri- 
lateral {xQ,yQ, xi,yi) € F{D) is equal to the cross-ratio of the vertices of 
the corresponding parallelogram. In particular, one always has the trivial 
solution z{x) = p{x) for all x G V{T)). 

A very useful transformation of the cross-ratio system is given by the 
following construction. 

Definition 6.5 The Hirota system for a given labelling of directed edges 
9 : E{T>) C consists of the following equations for the function w : 
V{D) — > C, one for every quadrilateral face {xo,yo,xi,yi) S F{T>): 

0ow{xo)w{yo) + 9iw{yo)w{xi) - 0ow{xi)w{yi) - 9iw{yi)w{xQ) = 0. (6.10) 

Note that this Hirota system is different from eq. 1)4. 3|) of Sect. 14.11 in that 
the present version has parameters defined on directed edges. In terms of the 
parallelogram realization p : V{T>) ^ C of the quad-graph D corresponding 
to the labelling 9, eq. (|6.1flj) reads: 

w{xQ)w{yQ){p{yQ) -p{xq)) + w{yQ)w{xi){p{xi) -p{yo)) + 

w{xi)w{yi){p{yi) - p{xi)) + w(yi)w{xo) {p{xo) - p{yi)) = 0. (6.11) 

Obviously, a transformation w — > cw on V{^) and w — > c'^^w on V{'5*) with 
a constant c G C, called hereafter a black- white scaling, maps solutions of 
the Hirota system into solutions. A relation between the cross-ratio and the 
Hirota system is based on the following observation: 

Theorem 6.6 Let w : V{T)) — > C 6e a solution of the Hirota system. Then 
the relation 

z{y) - z{x) = 9{x, y)w{x)w{y) = w{x)w{y) [p{y) - p{x)) (6.12) 

for all directed edges (x,y) G E{T)) correctly defines a unique (up to an ad- 
ditive constant) function z : V{T)) — > C which is a solution of the cross-ratio 
system i6.y\) . Conversely, for any solution z of the cross-ratio system \d. y\) . 
relation \6.1Bfl defines a function w : V{T)) — > C correctly and uniquely (up 
to a black-white scaling); this function w solves the Hirota system i6.1(J\) . 



6.3. INTEGRABLE CIRCLE PATTERNS 



137 



In particular, the trivial solution z{x) = p{x) of the cross-ratio system 
corresponds to the trivial solution of the Hirota system, w{x) = 1 for all 
X € V(T>). By a direct computation one can establish the following funda- 
mental property. 

Theorem 6.7 The Hirota system llf).l(J\) is 3D consistent. 

As a usual consequence, the Hirota system admits Backlund transformations 
and possesses zero curvature representation with transition matrices along 
the edge (x,y) € E{T)) given by 

1 -Ow{y) \ 

L{y,x,6;X)= ( , where 0=p{y)—p{x). 

-X9/w{x) w{y)/w{x) I 

(6.13) 



6.3 Integrable circle patterns 

Returning to circle patterns, let {z{x) : x G ^(S)} be intersection points 
of the circles of a pattern, and let {z{y) : y G ^(S*)} be their centers. 
Due to eq. ()6.1|) . the function z : V{T)) ^ C satisfies a cross-ratio system 
with Q : E{9) U E{9*) defined as Q(e) = exp(2i(/)(e)). Because of 

eq. 1)6.2(1 . the first one of the integrability conditions ()6.8() is fulfilled for an 
arbitrary circle pattern. Therefore, integrability of the cross-ratio system for 
circle patterns with the prescribed intersection angles (p '■ Ei9*) — > (0, vr) is 
equivalent to: 

n exp(2#(e*)) = 1, yyoeV{9*), (6.14) 

e*estar(j/o;S*) 

This is equivalent to the existence of the edge labelling a : E{T)) — > C such 
that, in notations of Fig. 16.2) 

exp{2i(P*) = ^. (6.15) 
ai 

Moreover, one can assume that the labelling a takes values in S^. 

Our definition of integrable circle patterns will require somewhat more 
than integrability of the corresponding cross-ratio system. 

Definition 6.8 A circle pattern with the prescribed intersection angles (j) : 
E{9*) — > (0,7r) is called integrable, if 



n exp(i0(e*)) = 1, Vyo€l^(g*), (6.16) 

e*Gstar{yo;S*) 
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i.e., if for any circle of the pattern the sum of its intersection angles with 
all neighboring circles vanishes (mod 27r). 

This requirement is equivalent to a somewhat sharper factorization than 
H6.15() . namely, to the existence of a labelling of directed edges 9 : E{D) — > 
such that, in notations of Fig. 16.21 

exp(i0) = ^ <^ exp{i(l)*) = (6.17) 

(Of course, the last condition yields (|6.15j) with a = 6'^.) The parallelogram 
realization p : V{T>) — > C corresponding to the labelling 9 € is actually a 
rhombic one. 

Theorem 6.9 Combinatorial data S CLn-d intersection angles (j) : -E'(S) — > 
(0, vr) belong to an integrable circle pattern, if and only if they admit an 
isoradial realization. In this case, the dual combinatorial data S* and inter- 
section angles (p ■ E{S*) (0, vr) admit a realization as an isoradial circle 
pattern, as well. 

Proof. The rhombic realization p : V(T)) — > C of the quad-graph D cor- 
responds to a circle pattern with the same combinatorics and the same 
intersection angles as the original one and with all radii equal to 1, and, 
simultaneously, to an analogous dual circle pattern. □ 

Consider a rhombic realization p : V{T)) — > C oi T>. Solutions z : 
V{T>) ^ C of the corresponding integrable cross-ratio system which come 
from integrable circle patterns are characterized by the property that the 
z-image of any quadrilateral {xo,yo,xi,yi) from F(T)) is a kite with the 
prescribed angle (j) at the black vertices z{xo), z{xi) (cf. Fig. 16. 2j) . It turns 
out that the description of this class of kite solutions admits a more con- 
venient analytic characterization in terms of the corresponding solutions 
w : V{T>) ^ C of the Hirota system defined by eq. (|6.12|) . 

Theorem 6.10 The solution z of the cross-ratio system corresponds to a 
circle pattern, if and only if the solution w of the Hirota system, correspond- 
ing to z via id.l^) . satisfies the condition 

w{x) G S\ w{y) e M+, Vx e V{9), y € V{5*). (6.18) 

The values w{y) € M-|_ have then the interpretation of the radii of the circles 
C{y), while the (arguments of the) values w{x) € measure the rotation of 
the tangents to the circles intersecting at z{x) with respect to the isoradial 
realization of the pattern. 
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Proof. As easily seen, the kite conditions are equivalent to: 
l^(^°)l-l and ^(^o)^™ 



\w{xi)\ w{yi) 
This yields (|6.18|) . possibly upon a black- white scaling. □ 

The conditions (|6.18j) form an admissible reduction of the Hirota system 
with 6 € S^, in the following sense: if any three of the four points w{xq), 
w{yo), w{xi), w{yi) satisfy the condition (|6.18|) . then so does the fourth one. 
This is immediately seen, if one rewrites the Hirota equation ()6.1U|1 in one 
of the equivalent forms: 

wjxi) ^ 6*1^(^1) - 6ow{yo) ^ w{yi) ^ 9ow{xo) +9iw{xi) 
w{xo) Oiw{yo) - 9ow{yi) w{yo) 9ow{xi) + eiw{xo)' 

As a consequence of this remark, we obtain Backlund transformations for 
integrable circle patterns. 

Theorem 6.11 Let all G S^, and let p : VCD) C be the corresponding 
rhombic realization of D. Let the solution w : V{T>) ^ C of the Hirota 
system correspond to a circle pattern with the combinatorics of ^, i.e., satisfy 
i6.18\) . Consider its Backlund transformation : VCD) C with an 
arbitrary parameter A € S"*^ and with an arbitrary initial value w^{xo) € M+ 
or w~^{yQ) € S"^. Then there holds: 

w+{x)£R+, w+{y)£S\ Vx G F(g), y G F(g*), (6.20) 

so that corresponds to a circle pattern with the combinatorics of 9*; 
which we call a Backlund transform of the original circle pattern. 

We close this section by mentioning several Laplace type equations which 
can be used to describe integrable circle patterns. First of all, the restriction 
of the function z to V{S) (i.e., the intersection points of the circles) satisfies 
the equations 

n 

fr[z{xk) - z{xo) 

Here z{xq) is any intersection point, where n circles C{yi), . . . , C{yn) meet, 
z{xk) are the second intersection points of C{yk) with C{yk+i), and Uk are 
the labels on the edges {xQ,yk) G ECD). Analogously, the restriction of the 
function z to V{S*) (i.e., the centers of the circles) satisfies the equation 

m 

^ z{yj) - z{yo) 



140 



CHAPTER 6. INTEGRABLE CIRCLE PATTERNS 



Here z{yo) is the center of any circle C{y()), which intersects with m circles 
C{yi), . . . ,C{ym) with the centers at the points z{yj); the intersection of 
C{yo) with C{yj) consists of two points and aj are the labels 

on the edges (yoi Xj) G E{D). These both Laplace type equations follow from 
the first claim of Theorem 14.141 applied to the cross-ratio system, which is 
nothing but the case (Ql)5=o of Theorem 14.151 

A similar construction can be applied to the Hirota system, written in 
the three- leg form (|6.19|) . Again, it yields two multiplicative Laplace type 
equations - on 9 and on 9*- It is instructive to look at the equation on 9* 
(for the radii rj = w{yj) of the circles): 



Due to eq. (|6.17l) . this equation can be written in terms of the intersection 
angles 4)j of C{yQ) with C{yj), and takes the form of eq. (|6.4|) . Interest- 
ingly, the latter equation holds for any circle pattern and is not specific for 
integrable ones (as opposed to the similar Laplace type equation on 9)- 

6.4 z"' and logz circle patterns 

Due to the 3D consistency of the cross-ratio and the Hirota systems, we can 
follow the procedure of Sect. 14.41 and to extend solutions to this systems 
from a quasicrystallic quad-graph D, realized as a quad-surface Q-j) C Z'^, to 
the whole of Z*^ (more precisely, to the hull of Ox>)- Then, one can ask about 
isomonodromic solutions. As shown in AgB^ iBobHoj . this leads to discrete 
analogs of the power function. Naturally, these discrete power functions are 
defined on the same branched covering S of the set Um=i C Z"^ as the 
discrete logarithmic function of Sect. 15.51 

Definition 6.12 For a € (0,1), the discrete z^"" is the solution of the 
discrete cross-ratio system on S defined by the values on the coordinate semi- 
axes zii^ = z{ner), r G [m,m -\- d — 1], which solve the recurrent relation 



m 



Ojro - Oj-irj 



n 



= 1. 



(^n+l Zn)(^Zn Zn—\) 



(6.21) 



n 



= azn- 



Zn+1 Zyi—\ 



with the initial conditions. 



= z(0) = 0, ) = z{er) = e 



exp(2alog 9r), 



(6.22) 



where logOj. is chosen in the interval h5.2b}) . 
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By induction, one can derive the following explicit expressions for the solu- 
tions Zn^: 

'S^ul^-~e, 4;'«^ni^-e. (..23, 

k=l k=l 

Observe the asymptotic relation 

= c{a){n0rf^{l + 0{n-^)). (6.24) 

The functions correspond to (intersection points and circle centers 
of) circle patterns. In order to prove this, it is convenient to study the 
corresponding solutions of the Hirota equation. Therefore, we introduce the 
functions w'^"'~^ related to z^" by 

z{n + Ej) — z{n) = 0jw{n)w{n + ej). 



Definition 6.13 For a G (0,1), the discrete w^"^ is the solution of the 
discrete Hirota system on S defined by the values on the coordinate semi- 
axes Wn^ = w{ner), r € [m, m + d — l], which solve the recurrent relation 



n 

Wn+l + Wn-1 

with the initial conditions: 



(a-l)(l-(-l)"). (6.25) 



wl^^ = w{0) = 0, tuf ^ = w{er) = el"-^ = exp((2a - 1) log Or), (6.26) 
where log Or is chosen in the interval \5.2b]) . 
One can easily find a closed expression for Wn^: 

= n ^-^^ , 4ii = e-^. (6.27) 

k=l 



Observe the asymptotics at n — > 00, 

4;) = c(a)n2-i| 



l + 0(n-i)). (6.28) 



The main technical advantage of the w variables is seen from the following 
observation. 
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Theorem 6.14 The function w ~ takes values from R-^- at the white points 
and values from at the black points. Therefore, the function z'^"' defines 
a circle pattern. 

Proof. The claim for w'^"'~^ on the coordinate axes is obvious from the 
explicit formulas (|6.27jl . and can be extended to the whole of S according to 
the remark after Theorem l6.1Ul The statement for z^'^ is now a consequence 
of Theorem I6.1UI □ 

The restriction of z"^"" to various quad-surfaces J^n give the discrete 
analogs of the power function on the corresponding quasicrystallic quad- 
graphs D with the set Q = {±^i, . . . , ±0^} of edge slopes (see Fig. I6.4j) . 



Figure 6.4: Circle patterns z^/^ with the combinatorics of the square grid, 
and z^/'^ with the combinatorics of the regular hexagonal lattice (isotropic 
and non-isotropic) . 



6.4. AND LOG Z CIRCLE PATTERNS 



143 



These pictures lead to the conjecture |BobP3j that the circle patterns z'^^ 
are embedded, i.e., interiors of different kites are disjoint. This conjecture 
has been proven in |Ag| for the case of the square grid combinatorics. The 
fact that the circle patterns z^" are immersed, i.e., the neighboring kites do 
not overlap, was proven in |AgBl| for the square grid and in AgB2 for the 
hexagonal grid combinatorics. One possible approach to the general case 
could be based on applying the well developed techniques of the theory of 
isomonodromic solutions |INj to the asymptotic study of the discrete z'^"', 
because of the following statement. 

Theorem 6.15 The discrete z^" is isomonodromic: for a proper choice of 
A{0; A), the matrices A(n; A) at any point n E have simple poles only: 

A{n- A) = + — i^i , (6.29) 



with 



,„^ / —a/2 —az(n) i 

A^'^\n) = { ' ^ M , (6.30) 

a/2 



n 



ni 



z{n + ei) - z{n - e/) 
z{n + ei) - z{n) {z{n + ei) - z(n))(z(n) - z{n - e/)) 
1 z{n) — z{n — El) 



. (6.31) 



At any point n ^ S, the discrete z^"" satisfies the following constraint: 

fri z{n + e,) - z{n - e,) 

Proof follows the same scheme as the proof of Theorem 15. 141 one first shows 
that the poles of A{ner; A) remain simple, due to the recurrent relations 
()6.21|) . and then shows that the order of poles does not increase at the points 
n away from the coordinate axes, due to the multidimensional consistency. 
Note that this scheme works also for Theorem 16.161 □ 

The transition between z and w variables is a matter of straightforward 
computations. Actually, the next theorem is dealing with the same matrices 
as Theorem 16. 151 but written in different variables. 
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Theorem 6.16 The discrete w'^"'~^ is isomonodromic: for a proper choice 
of A{0;X), the matrices A{n; X) at any point n € C^+Y have simple poles 
only: 

A{n; A) = + J] — ^ , (6.33) 

1=1 ^~ 



with 



A(0)(n)= I ""/^ * 1 , (6.34) 




^(i), . ni w{n + ei) eiw{n + ei)w{n - ei) 



w{n + ei) + w{n - ei) \ ijQ^ w{n - ei) 

(6.35) 

The upper right entry of the matrix A^^^n), denoted by the asterisk in 

i6.34\ ), is given by A^2''(^) = ~ X^f=i ^12 (^)- ^^2/ point n € S, the 
discrete w'^"'''^ satisfies the following constraint: 

^ w{n + ei) + w{n - ei) V 2/v v ; j \ j 

It is interesting do study the limiting behaviour of the function as 
a ^ 0. It is not difficuh to see that for all n 7^ one has z^°(n) 1. 
Denote 

L(n) = lim^-^ . (6.37) 

a^O za 

This function is called the discrete logarithmic function; it should not be con- 
fused with the namesake function £(n) in the linear theory (Sect. 15. From 
eq. ()6.37p the following characterization is found: the discrete logarithmic 
function L is the solution of the discrete cross-ratio system on S defined by 
the values on the coordinate semi-axes L^'^ — L(Tie^i), r (£ [m,m -\- d — 1], 
which solve the recurrent relation 

^ (Ln+i — L„)(L„ — L„_i) _ 1^ 
Ln+i — L„_i 2 

with the initial conditions: 

lJ,"^) = L(0) = 00, hf^ = L(e^) = log^r, (6.39) 



6.5. LINEARIZATION 



145 



where log^,. is chosen in the interval ()5.26|) . Exphcit expressions: 



n-l 



lS = log + E ^ + i ' 41l = log ^'-'"^ 

k=l k=l 



Theorem 6.17 The discrete logarithm is isomonodromic and satisfies, at 
any point n £ S, the following constraint: 

" (L(n + e,)-L(n))(L(n)-L(n-e,)) _ 1 

' L(n + e,)-L(n-e,) "2' ^^'^^^ 

By restriction to quad-surfaces Qx>, we come to the discrete logarithmic 
function on arbitrary quasicrystallic quad- graphs T>. By construction, they 
all correspond to circle patterns. A conjecture that these circle patterns are 
embedded seems plausible (see Fig. IH.5|) . 



Figure 6.5: Discrete logarithm circle patterns with the combinatorics of the 
regular square and hexagonal lattices 



6.5 Linearization 

Let e : E{T)) ^ C be an edge labelhng, and let p : V(Ti) ^ C be the 
corresponding parallelogram realization of T) defined hy p{y)—p{x) = 6{x, y). 
Consider the trivial solutions 



zq{x)=p{x), wq{x) = 1, Mx £V{Ti) 
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of the cross-ratio system ()6.9p and the corresponding Hirota system ()6.11() . 
Suppose that zq : V{T>) C belongs to a differentiable one-parameter 
family of solutions Ze : ^C^) — > C, e € (— eo,eo), of the same cross-ratio 
system, and denote by ■ ^(25) — > C the corresponding solutions of the 
Hirota system. Denote 



Theorem 6.18 Both functions f,g : V{T)) — > C solve discrete Cauchy- 
Riemann equations i5.16]) . 

Proof. By differentiating 1)6. 12() , we obtain a relation between the functions 



9{y)-9{x) = {f{x) + f{y)){p{y)-p{x)), y{x,y) e E{V). (6.43) 



The proof of proposition is based on this relation solely. Indeed, the closeness 
condition for the form on the right-hand side reads: 



{f{xo) + f{yo)) {Pivo) - p{xo)) + (/(yo) + /(xi)) (p(xi) - p{yo)) + 

{f{xi) + /(yi)) {p{yi) - p(xi)) + (/(yi) + /(xq)) {p{xo) - p{yi)) = 0, 



which is equivalent to 1)5. 16(1 for the function /. Similarly, the closeness 
condition for /, that is, 



(/(xo) + /(yo)) - (/(yo) + f{xi)) + (/(xi) + /(yi)) - (/(yi) + /(xq)) = 0, 



gjyo) -gjxo) _ g{xi) -g{yo) ^ g{yi) - gjxi) _ g{xo) - g{yi) ^ ^ 
p(yo)-p(xo) p{xi)-p{yo) p{yi)-p{xi) p{xo)-p{yi) 

Under the condition p{yo) —p{xo) = p{xi) —p{yi), this is equivalent to ()5.16() 
for g. □ 

Remark. This proof shows that, given a discrete holomorphic func- 
tion / : V{T)) — > C, relation (|6.43|) correctly defines a unique, up to an 
additive constant, function g : VCD) C, which is also discrete holomor- 
phic. Conversely, for any g satisfying the discrete Cauchy-Riemann equa- 
tions (|5.16|) . relation (|6.43|) defines a function / correctly and uniquely (up 
to an additive black- white constant); this function / also solves the discrete 





(6.42) 



yields: 
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Cauchy-Riemann equations ()5.16() . Actually, formula H6.43() expresses that 
the discrete holomorphic function / is the discrete derivative of g, so that 
g is obtained from / by discrete integration. This operation was considered 
in tDo1,,Do2. Mel. . 

Summarizing, we have the following statement. 

Theorem 6.19 a) A tangent space to the set of solutions of an integrahle 
cross-ratio system, at a point corresponding to a rhombic embedding of a 
quad-graph, consists of discrete holomorphic functions on this embedding. 
This holds in both descriptions of the above set: in terms of variables z 
satisfying the cross-ratio equations, and in terms of variables w satisfying 
the Hirota equations. The corresponding two descriptions of the tangent 
space are related by taking the discrete derivative (resp. anti- derivative) of 
discrete holomorphic functions. 

b) A tangent space to the set of integrable circle patterns of a given combi- 
natorics, at a point corresponding to an isoradial pattern, consists of discrete 
holomorphic functions on the rhombic embedding of the corresponding quad- 
graph, which take real values at white vertices and purely imaginary values 
at black ones. This holds in the description of circle patterns in terms of 
circle radii and rotation angles at intersection points (Hirota system). 

A spectacular instance of this linearization property is delivered by the 
isomonodromic discrete logarithm studied in Sect. 15.51 and isomonodromic 
z^°' circle patterns of Sect. 16.41 

Theorem 6.20 The tangent vector to the space of integrable circle patterns 
along the curve consisting of patterns w'^"'~^ , at the isoradial point corre- 
sponding to a = 1/2, is the discrete logarithmic function £ of Sect. 15.51 

Proof. We have to prove that the discrete logarithm i and the discrete 
power function w'^'^~^ are related by 



Due to Theorem lf).18( it is enough to prove this for the initial data on the 
coordinate semi-axes. But this follows by differentiating with respect to a 
the initial values (|6.27j) at the point a = 1/2, where all li; = 1: the result 
coincides with 1)5. 31() . □ 
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